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respectively .  19 
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Fig.  5.  Vertical  profile  of  Doppler-derived  u  (solid  curve)  and  v  (dashed  curve)  wind 
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Fig.  6.  Vertical  profile  of  Doppler-derived  u  (solid  curve)  and  v  (dashed  curve)  wind 
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elevation  is  relative  to  sea  level  in  dekameters  and  contour  intervals  are  every  20 
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Fig.  8.  Horizontal  cross-section  of  observed  radar  reflectivities  in  dBZ  and  vector 
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Fig.  9a.  Horizontal  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
as  in  Fig.  8,  except  at  1 134  MST  and  at  3.3  km  MSL . 32 
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(dashed)  contours  show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig. 
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eveiy  3.0  m  s  '  starting  at  1.5  m  s  '  for  CONTROL  simulation  at  6480  s.  Solid 
(dashed)  contours  show  westerly  (easterly)  winds.  Axes  are  ladled  as  in  Fig.  27. 
. 98 
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Fig,  33.  As  in  Fig.  27,  east-west  vertical  cross-section  of  simulated  total  condensate 
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mixing  ratios  for  CONTROL  simulation  at  8640  s.  Contour  label  "  1 "  represents  0. 1 
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Fig.  41.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  as  in  Fig.  27,  except  for  CONTROL  simulation 
at  8820  s.  Contour  label  "1"  represents  0.1  g  kg'‘ .  122 
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Fig.  46a.  Horizontal  cross-section  of  observed  radar  reflectivities  as  in  Fig.  10a,  ex¬ 
cept  at  1134  MST  and  4.3  km  MSL.  Line  AB  shows  location  of  vertical  cross-section 
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Fig.  49.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the 
CONTROL  simulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash 
curve  shows  cloud  water  mixing  ratios,  dashed  curve  shows  rain  water  mixing  ratios, 
double  dot-dash  curve  shows  pristine  ice  crystal  mixing  ratios,  shaded  curve  shows 
graupel  mixing  ratios,  and  dotted  curve  shows  aggregated  snow  flake  mixing  ratios. 
Vertical  axis  is  mixing  ratio  in  grams  per  kilogram.  Horizontal  axis  is  simulation 
time  in  seconds  (bottom)  and  in  hours  and  seconds  MST  (top) . 143 
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every  2.0  g  kg'‘  starting  at  0.1  g  kg  '.  Wind  vector  length  shown  in  lower  right-hand 
comer  of  figure  equals  SO  m  s  '.  Vertical  axis  is  labeled  in  kilometers  above  the 
minimum  terrain  elevation  in  the  simulation  grid  (in  this  case  1.82  km)  as  described 
in  the  text.  Horizontal  axis  is  horizontal  distance  labeled  in  kilometers.  The  cross- 
section  shown  is  located  at  >'  =  0.0  km  and  runs  east-west,  with  the  origin  at  Lang¬ 
muir  Laboratory .  147 

Fig.  5 1 .  East-west  vertical  cross-section  of  simulated  vertical  velocity  contoured  ev¬ 
ery  3.0  m  s  '  starting  at  1.5  m  s  '  for  EXPl  simulation  at  6300  s.  Solid  (dashed)  con¬ 
tours  show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig.  50. 
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Fig.  52.  Analogous  to  Fig.  50,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXPl  simulation  at 
7380  s.  Contour  label  "1"  represents  0.1  g  kg  ' .  150 

Fig.  53.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  as  in  Fig.  50,  except  for  EXPl  simulation  at 
8100  s.  Contour  label  "1"  represents  0.1  g  kg"' .  152 

Fig.  54.  Analogous  to  Fig.  51,  east-west  vertical  cross-section  of  simulated  vertical 
velocity  for  EXPl  simulation  at  8100  s .  153 

Fig.  55.  As  in  Fig.  50,  east-west  vertical  cross-section  of  simulated  total  condensate 
mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXPl  simulation  at  8820  s. 
Contour  label  "1"  represents  0.1  g  kg  ' .  155 
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Fig.  57.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the 
EXPl  simulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash  curve 
shows  cloud  water  mixing  ratios,  dashed  curve  shows  rain  water  mixing  ratios,  double 
dot-dash  curve  shows  pristine  ice  crystal  mixing  ratios,  shaded  curve  shows  graupel 
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Fig.  58.  Time  evolution  of  maximum  and  minimum  simulated  vertical  velocities  for 
the  EXPl  simulation.  Solid  curve  ^ows  maximum  upward  vertical  velocities 
(positive  w).  Dashed  curve  shows  maximum  downward  vertical  velocities  (negative 
vt).  Vertical  axis  is  vertical  velocity  in  m  s  '.  Horizontal  axis  is  simulation  time  in 
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section  shown  is  located  at  >'  =  0.0  km  and  runs  east-west,  with  the  origin  at  Lang¬ 
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Fig.  63.  Analogous  to  Fig.  59,  east-west  vertical  cross-section  of  simulated  total  con¬ 
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Fig.  68.  Time  evolution  of  maximum  and  minimum  simulated  vertical  velocities  for 
the  EXP4  simulation.  Solid  curve  shows  maximum  upward  vertical  velocities 
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Fig.  72.  Analogous  to  Fig.  70,  east-west  vertical  cross-section  of  simulated  total  con¬ 
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Fig,  79.  Analogous  to  Fig.  71,  east-west  vertical  cross-section  of  simulated  vertical 
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Fig.  80.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the 
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Marie  Edwin  Raffensberger,  M.S. 

The  Florida  State  University,  1993 
Major  Professor.  Peter  Sawin  Ray,  Ph.D. 


A  two-dimensional,  non-hydrostatic  cloud  model  was  used  to  isolate  the  processes 
and  initial  conditions  most  important  in  the  initiation  and  development  of  a  small 
mountain  thunderstorm.  Six  numerical  simulations  were  conducted-one  control  and 
five  experiments.  The  control  simulation  was  conducted  with  realistic  initial  condi¬ 
tions  and  physical  processes.  The  simulation's  accuracy  was  evaluated  by  comparison 
to  multiple-Doppler  radar  analyses  of  a  storm  that  occurred  on  31  July,  1984  near 
Langmuir  Laboratory,  New  Mexico  and  to  a  microphysical  retrieval  conducted  for  a 
different,  but  similar  storm.  The  other  five  simulations  were  conducted  to  test  the 
sensitivity  of  the  simulated  storms  to  the  initial  wind  profile,  the  lack  of  solar  heating, 
restriction  to  warm  rain  processes,  and  the  initial  moisture  profile.  Comparison  of  the 
evolution  of  the  control's  simulated  storm  with  observations  and  with  the  microphysi¬ 
cal  retrieval  show  that  the  simulatirm  accurately  captured  the  major  features  observed 
by  the  radars  and  accurately  depicted  the  microphysical  evolution  of  the  storm  cells. 
The  results  of  the  numerical  experiments  showed  the  following  conditions  and 
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processes  to  be  most  important  (in  descending  order)  to  the  development  and  mainte¬ 
nance  of  the  mountain  storm:  solar  heating,  ice  phase  microphysical  processes,  envi¬ 
ronmental  wind  profile,  and  low-level  moisture  profile.  The  environmental  wind 
profile  was  the  single  most  important  factor  in  the  convective  initiation  location. 


NUMERICAL  SIMULATIONS  OF  A  MOUNTAIN  THUNDERSTORM: 
A  COMPARISON  WITH  DOPPLER  RADAR  OBSERVATIONS 


Chapter  1 
Introduction 

Mountain  thunderstorms  have  been  studied  extensively  to  determine  the  processes 
and  conditions  important  to  their  initiation,  development,  and  maintenance.  Mountain 
thunderstorms  play  a  major  role  in  the  economy  of  the  western  United  States  as  they 
produce  the  bulk  of  summertime  precipitatimi,  and  their  lightning  starts  many  of  the 
region's  forest  fires.  The  extensive  research  into  mountain  thunderstorms  has  been 
prompted  by  the  need  to  improve  the  understanding  and  prediction  of  the  storms  in 
order  to  improve  precipitation,  flash  flood,  lightning,  and  civilian  and  military  avi¬ 
ation  forecasts.  In  addition  to  their  economic  importance,  the  frequency  of  occur¬ 
rence  of  mountain  storms  over  the  same  general  location  makes  them  ideal  subjects 
for  intense  scientific  study. 

Over  the  past  40  years,  numerous  descriptive  models  of  mountain  thunderstorms 
were  developed  based  on  observational  studies  which  used  i^otogrammetric  tech¬ 
niques,  instrumented  aircraft,  Doppler  radars,  and  dense  observational  networks. 
Within  the  past  30  years,  two-  and  three-dimensional  numerical  models  have  been  de¬ 
veloped  which  allow  researchers  to  simulate  mountain  thunderstorms.  These  models 
permit  "experimentation"  with  the  storms  by  allowing  the  researcher  to  vary  impor¬ 
tant  physical  processes  and  initial  conditions  in  a  way  not  permitted  in  nature. 
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The  current  research  examines  the  initiation,  development,  and  maintenance  of  a 
small  mountain  thunderstorm  using  numerical  simulations  and  observaticNially-based 
information.  The  storm  considered  in  this  study  occurred  on  31  July  1984  during  the 
Thunderstorm  Research  International  Project  (TRIP),  a  field  experiment  conducted 
near  the  Langmuir  Laboratory  in  the  Ms^dalena  Mountains  of  central  New  Mexico 
during  the  summer  of  1984.  Four  Doppler  radars,  two  instrumented  research  aircraft, 
rawinsondes,  and  a  surface  network  for  measuring  electrical  fields  were  employed  to 
collect  storm  data  during  TRIP. 

This  study  employs  the  two-dimensional,  non-hydrostatic  cloud  model  configura¬ 
tion  of  the  Colorado  State  University  Regional  Atmospheric  Modeling  System 
(RAMS)  (Tripoli  and  Cotton,  1982,  1986;  Cotton  et  al.,  1982,  1986;  Tripoli,  1986) 
to  simulate  the  31  July  1984  storm.  The  study  also  uses  multiple  Doppler  analyses 
performed  following  the  technique  of  Ray  et  al.  (1980)  and  described  in  detail  by 
Lang  (1991),  as  well  as  the  three-dimensional  microphysical  retrievals  performed  by 
Lang  ( 1991)  for  a  similar  storm  that  occurred  on  3  August  1984. 

This  research  has  two  primary  goals:  1)  to  assess  the  validity  of  a  RAMS  simu¬ 
lation  compared  with  observations  from  31  July  1984,  and  2)  to  use  RAMS  to  conduct 
numerical  simulation  experiments  to  test  the  sensitivity  of  the  simulated  storms  to  dif¬ 
ferent  physical  processes  and  initial  conditions.  A  total  of  six  simulations  were  con¬ 
ducted.  The  first  simulation  was  the  control  and  included  all  appropriate  physical 
parameterizations  and  processes  and  was  initialized  with  observationally-based  base 
state  atmospheric  profiles.  The  results  of  the  control  simulation  were  compared  with 
Doppler  analyses  and  with  Lang's  (1991)  retrieval  to  assess  the  validity  of  the  simu¬ 
lation  and  its  usefulness  in  studying  the  storm.  The  other  five  simulations  were  con¬ 
ducted  to  test  the  sensitivity  of  the  simulated  storms  to  the  initial  wind  profile,  the 
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lack  of  solar  heating,  restriction  to  warm  rain  processes,  and  the  initial  moisture  pro¬ 
file.  The  results  of  the  numerical  experiments  were  compared  with  those  of  the  con¬ 
trol  to  determine  the  impact  of  the  change  imposed  for  each  experiment. 

The  thesis  is  organized  in  the  following  manner.  Chapter  2  provides  an  historical 
review  of  both  observational  and  numerical  mountain  thunderstorm  studies.  Chapter 
3  describes  the  observations  of  the  31  July  1984  storm  and  its  environment.  This 
chapter  includes  discussions  of  the  data  collection  and  analysis  procedures  and  a  de¬ 
tailed  description  of  the  evolution  of  the  storm  as  observed  by  the  Doppler  radars. 
Chapter  4  describes  the  numerical  model  as  it  was  configured  for  this  study.  Chapter 
5  explains  the  experimental  design  chosen  for  this  study.  This  chapter  includes  in¬ 
formation  on  the  construction  of  the  base-state  profiles  used  to  initialize  the  simu¬ 
lations.  Chapter  6  details  the  results  of  the  control  simulation  and  compares  those 
results  with  the  Doppler  analyses,  synthesized  wind  fields,  and  Lang’s  (1991)  micro¬ 
physical  retrieval  for  the  3  August  1984  case.  Chapter  7  describes  the  results  of  the 
experiments  conducted  to  test  the  sensitivity  of  the  simulated  storms  to  the  initial  wind 
profile.  Chapter  8  presents  the  results  of  the  experiments  conducted  to  test  the  sensi¬ 
tivity  of  the  simulated  storms  to  diabatic  heating  effects;  namely  elimination  of  solar 
heating  and  restriction  to  warm  rain  processes.  Chapter  9  describes  the  results  of  the 
experiment  conducted  to  test  the  sensitivity  of  the  simulated  storms  to  the  initial  mois¬ 
ture  profile.  Finally,  Chapter  10  summarizes  with  the  conclusions. 


Chapter  2 

Historical  Review  Of  Mountain  Thunderstorm  Studies 
2.1  Observational  Studies 

Mountain  thunderstorms  are  believed  to  result  from  thermally  driven  mountain- 
valley  circulations  (Defant,  19S1).  Solar  heating  of  the  mountain  slopes  causes  the 
adjacent  air  to  become  warmer  than  ^ee  air  at  the  same  elevation.  The  consequent 
buoyant  rise  of  the  warmer  air  to  the  condensation  level  produces  clouds  (Raymond 
and  Wilkening,  1980).  At  this  point  the  effects  of  latent-heat  release  mask  the  origi¬ 
nal  circulation  and  may  soon  become  the  dominant  mechanism  for  cloud  and  storm 
maintenance. 

Banta  and  Schaaf  (1987)  presented  several  other  mechanisms  thought  to  be  respon¬ 
sible  for  the  production  of  thunderstorms  in  the  mountains,  the  first  of  which  was  oro¬ 
graphic  lifting.  Orographic  lifting  involves  the  mechanical  lifting  of  air  as  it 
approaches  a  barrier.  Wake  effects  were  the  second  mechanism.  They  include  wake 
turbulence,  gravity  waves  and  confluence  generated  by  flow  around  a  mountain.  Oth¬ 
er  mechanisms  described  include:  leeside  convergence  (discussed  below),  flow  into 
convergent  valleys,  and  other  mechanisms  not  dependent  on  topography  (e.g.,  fronts, 
local  gradients  of  moisture,  etc.) 

Environmental  conditions  may  play  an  important  role  in  determining  which  of  the 
above  processes,  if  any,  dominate  the  development  of  mountain  convection.  Among 
the  environmental  conditions  to  be  considered  are  ambient  wind,  temperature,  and 
humidity  profiles  and  surface  characteristics  such  as  vegetation,  soil  type,  and  soil 
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moisture.  Many  observational  and  numerical  studies  have  been  conducted  to  docu¬ 
ment  the  nature  of  mountain  convection  and  to  determine  the  processes  and  conditions 
important  to  its  initiation  and  maintenance.  In  order  to  concentrate  on  the  thermally 
driven  circulations  leading  to  convection,  most  of  these  studies  focused  on  conditions 
of  light  winds  and  strong  sc^  heating.  A  review  of  some  of  these  studies  is  given 
below  beginning  with  observational  studies  and  followed  by  numerical  studies. 

Several  of  the  early  observational  studies  were  not  motivated  purely  by  an  interest 
in  mountain  convection,  but  rather  were  designed  to  settle  a  long-standing  debate  over 
the  actual  nature  of  convection.  The  "bubble  theory"  (Scorer  and  Ludlam,  1953; 
Scorer,  1958)  held  that  convection  consisted  of  a  series  of  short  impulses,  or  bubbles, 
of  buoyant  air  rising  from  the  surface  with  mixing  occurring  mainly  at  the  tops  of  the 
bubbles.  The  "jet  theory",  on  the  other  hand,  held  that  convection  could  be  thought 
of  as  a  jet,  or  plume,  of  buoyant  air  with  a  continuous  source  of  buoyant  air  entering 
through  the  bottom  and  entrainment  occurring  at  the  sides.  The  observation  that  cu¬ 
mulus  clouds  have  flat  bases  was  difficult  to  explain  with  the  bubble  theory.  To  solve 
this  problem,  Turner  (1962)  proposed  that  convection  be  viewed  instead  as  a  "starting 
plume"  which  consisted  of  a  plume  of  buoyant  air  capped  by  a  rising  thermal.  By 
linking  the  jet  and  bubble  theories,  this  theory  allowed  for  a  continuous  source  of 
buoyant  air  through  the  bottom,  which  could  explain  the  existence  of  flat  cloud  bases. 
Because  mountain  convection  was  known  to  occur  frequently  over  the  same  location, 
it  soon  became  a  "natural  laboratory"  for  testing  these  theories.  The  information 
from  these  studies  showed  that  the  topic  of  mountain  convection  itself  was  one  worthy 
of  further  research. 

Braham  and  Draginis  (1960)  used  an  instrumented  aircraft  to  measure  temperature 
and  humidity  over  the  Santa  Catalina  Mountains  near  Tucson.  They  found  that  by 
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mid-moming  a  "convective  core”  measuring  about  1-4  km  in  diameter  had  formed 
over,  or  slightly  downwind  of  the  peaks.  In  contrast  to  the  relatively  dry  and  moder¬ 
ately  stable  environment,  the  core  was  characterized  by  an  excess  of  water  vapor  mix¬ 
ing  ratio  and  a  nearly  dry  adiabatic  lapse  rate.  Excesses  of  temperature  and  virtual 
temperature  also  were  noted  in  the  lower  portion  of  the  core.  Although  no  wind  mea¬ 
surements  were  collected,  negative  mixing  ratio  anomalies  and  positive  potential  tem¬ 
perature  anomalies  surrounding  the  core  at  higher  elevations  implied  that  downdrafts 
as  well  as  updrafts  were  integral  parts  of  the  observed  convective  system. 

Anderson  (1960)  studied  congestus  clouds  in  Arizona  with  photogrammetric  tech¬ 
niques.  The  study  showed  that  as  the  clouds  grew,  they  pulsated  with  a  frequency  of 
about  10  minutes.  Anderson  attributed  these  pulsations  to  the  buoyancy-restoring 
force  of  statically  stable  air  and  proposed  a  model  in  which  pulsating  flow  was  super¬ 
imposed  on  a  slowly  increasing,  large-scale  updraft. 

Glass  and  Carlson  (1963)  used  photogrammetric  techniques  to  study  small  cumu¬ 
lus  clouds  over  the  San  Francisco  Peaks  near  Flagstaff.  They  found  that  many  of  the 
clouds  in  their  study  seemed  to  fit  the  bubble  theory  of  convection.  The  clouds  in  this 
study  typically  measured  1  km  or  less  in  diameter.  Todd  (1964)  combined  aircraft 
observations  with  time-lapse  cloud  photography  to  trace  the  continuity  of  individual 
updrafts  over  the  San  Francisco  Peaks.  His  study  was  different  from  those  above  in 
that  it  followed  the  convective  development  through  to  the  precipitation  stage.  He 
found  that  many  updrafts  lasted  from  7  to  18  minutes  and  that  individual  updrafts  had 
maximum  velocities  two  to  three  times  greater  than  the  maximum  rate  of  rise  of  the 
convective  tower  tops.  Todd  used  these  observations  to  suggest  that  the  clouds  he  stu¬ 
died  followed  the  "starting  plume"  theory  of  Turner  (1962).  The  non-precipitating 
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clouds  in  this  study  also  typically  measured  less  than  1  km  in  diameter,  and  the  pre¬ 
cipitating  clouds  measured  from  1  to  3  km  in  diameter. 

In  another  photogrammetric  study,  Orville  (196Sb)  studied  convection  over  the 
Santa  Catalinas  when  the  low-level  wind  speeds  were  between  S  and  10  m  s'^  He 
found  that  the  location  of  convective  initiation  strongly  depended  on  the  direction  of 
the  environmental  wind.  Orville  speculated  that  this  dependency  resulted  from  the 
translation  of  the  updrafts  by  the  environmental  wind. 

In  their  study  of  dry  mountain  convection,  Raymond  and  Wilkening  (1980)  used 
an  instrumented  aircraft  to  collect  temperature,  humidity  and  wind  data  over  the  San 
Mateo  Mountains  near  Socorro,  New  Mexico.  By  carefully  analyzing  and  filtering  the 
data,  the  authors  isolated  two  important  scales  of  motion.  Superimposed  on  the  ambi¬ 
ent  flow  was  a  heat-island  circulation,  convergent  in  the  lower  levels  and  divergent 
aloft,  which  measured  approximately  20  km  in  diameter.  The  strength  of  this  mesos- 
cale  flow  averaged  over  its  area  ranged  from  about  one  meter  per  second  in  the  verti¬ 
cal  wind  to  a  few  meters  per  second  in  the  horizontal.  The  authors  compared  this 
feature  to  the  "convective  core"  discussed  earlier  by  Braham  and  Draginis  (1960). 
Embedded  within  the  convective  core  were  numerous  convective  eddies  with  horizon¬ 
tal  wavelengths  of  3-4  km  and  vertical  velocities  measuring  about  5  m  s  ^  The  eddies 
were  responsible  for  most  of  the  heat  and  moisture  fluxes  over  the  mountain.  The 
authors  felt  that  the  eddies  and  the  mesoscale  circulation  coexisted  in  a  "symbiotic 
relationship":  the  heat  carried  aloft  by  the  eddies  drove  the  mesoscale  circulation, 
and  the  mesoscale  circulation  carried  the  heat  out  of  the  region,  thus  preserving  a  fa¬ 
vorable  lapse  rate  for  the  eddies. 

In  a  follow-up  study  of  moist  convection  Raymond  and  Wilkening  (1982)  showed 
that  moisture  negligibly  affected  the  low-level  convergence  during  the  early  stages  of 
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convection.  They  found  that  the  dry-convective  pattern  of  convergence  below  the 
peaks  and  divergence  above  persisted  through  the  cumulus  stages  of  develc^ment  and 
into  the  early  cumulonimbus  stage.  Tliey  found  only  a  single  storm  that  developed  its 
own  low-level  flow  pattern.  The  authors  presented  a  schematic  diagram  of  the  esti¬ 
mated  mass  budget  of  a  typical  mountain  thunderstorm.  It  depicted  a  convergent  flow 
below  the  peaks,  which  supported  both  the  mesoscale  upcurrent  and  a  divergent  flow 
which  extended  from  above  the  peaks  to  about  half  the  height  of  the  cloud. 

In  a  paper  describing  the  evolution  of  the  daytime  boundary-layer  over  mountain¬ 
ous  terrain,  Banta  (1984)  documented  the  existence  of  a  "leeside  convergence  zone". 
He  showed  that  as  the  mountain  surface  heated  during  the  morning,  a  shallow  upslope 
flow  developed  near  the  slope  beneath  the  nocturnal  inversion  cold  pool.  At  eleva¬ 
tions  above  the  cold  pool  the  winds  were  convectively  mixed.  On  the  lee  side  of  the 
mountain  these  two  flows  opposed  each  other  and  thus  created  a  convergence  zone. 
Banta  felt  that  this  leeside  convergence  zone  might  explain  the  slight  displacement 
downwind  from  the  peaks  of  the  convective  elements  observed  by  Braham  and  Dragi- 
nis  (1960),  Glass  and  Carlson  (1963),  and  Orville  (1965b). 

Recently  several  investigators  (Barker  and  Banta,  1985;  Banta  and  Schaaf,  1987; 
Schaaf  et  al.,  1988)  have  used  geostationary  satellite  imagery  to  trace  mountain  con¬ 
vection  back  to  its  initiation  location.  The  authors  found  that  certain  topographical 
featured  act  as  preferred  storm  initiation  locations  or  "genesis  zones".  When  they 
stratified  the  initiation  locations  by  the  direction  of  the  ridge  top  wind,  the  authors 
found  that  the  wind  direction  strongly  affected  the  ability  of  a  particular  topographic 
feature  to  act  as  an  initiation  point.  All  three  of  these  studies  showed  the  same  ten¬ 
dency  toward  leeside  development  that  appeared  in  earlier  studies. 
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Rayin<xid  and  Wilkening  (1985)  used  budget  residual  techniques  to  study  thunder¬ 
storms  and  cumulus  congestus  clouds  near  Socorro,  New  Mexico.  As  part  of  the 
study  the  authors  compiled  a  composite  thermodynamic  sounding  for  thunderstorm 
and  congestus  cases  and  found  the  congestus  sounding  to  be  much  drier  at  mid-levels 
than  the  thunderstorm  sounding.  The  authors  also  found  that  the  peak  vertical  mass 
flux  for  thunderstorms  was  only  slightly  greater  than  that  for  ccmgestus  clouds.  They 
felt  this  ruled  out  a  significant  feedback  between  latent  heat  release  and  the  generation 
of  low-level  convergence.  This  implied  that  the  low-level  forcing  for  both  cases  was 
about  the  same.  They  also  found  that  for  thunderstorms  the  vertical  mass  flux  near 
700  mb  was  only  about  half  that  found  in  congestus  clouds.  This  difference  was  at¬ 
tributed  to  the  thunderstorm  downdraft,  which  was  estimated  to  be  about  one-half  the 
updraft  strength. 

Toth  and  Johnson  (1985)  used  wind  data  collected  from  a  dense  surface  mesonet- 
work  to  study  surface  flow  characteristics  over  the  Front  Range  in  northeast  Colora¬ 
do.  The  observations  showed  the  evolution  of  upslope  flow  during  the  morning  hours 
with  a  convergent  region  just  to  the  lee  of  the  ridges.  In  late  afternoon,  a  downslope 
flow  began  in  the  high  elevations.  The  confluence  zone  between  the  downslope  and 
the  upslope  regimes  gradually  moved  eastward  onto  the  plains.  The  authors  found 
that  this  propagating  confluence  line  coincided  with  a  maximum  radar  echo  frequency 
in  the  late  afternoon  over  the  western  plains. 

Ziegler  et  al.  (1986)  used  Doppler  radars  in  their  study  of  a  thunderstorm  near 
Socorro.  They  found  vertical  extensions  in  the  radar  reflectivity  field  that  occurred 
with  a  period  of  about  12  min.  These  upward  extensions  were  suggestive  of  pulses  in 
the  vertical  velocity  field  similar  to  those  observed  by  Anderson  (1960).  The  maxi¬ 
mum  updrafts  inferred  from  the  Doppler  velocities  were  between  15  and  20  m  s'*. 
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Additionally,  the  inferred  flow  field  reported  in  this  study  supported  the  convergent 
below/divergent  aloft  flow  field  reported  by  Raymond  and  Wilkening  (1980). 

McCutchan  and  Fox  (1986)  used  automatic  surface  weather  stations  to  collect 
wind,  temperature,  and  humidity  data  at  various  aspects  and  elevations  on  an  isolated 
conical  mountain  in  New  Mexico.  The  station  positions  allowed  the  authors  to  use 
analysis  of  variance  techniques  to  study  the  effects  of  elevation  and  aspect  on  the  mea¬ 
sured  variables.  The  data  were  divided  into  two  groups  of  wind  speeds  at  the  peak: 
light  (less  than  5  m  s'‘)  and  strong  (greater  than  or  equal  to  5  m  s  ‘).  Analysis  of 
variance  tests  showed  that  for  daytime  cases  in  both  groups:  1 )  the  effect  of  aspect  on 
the  u  and  v  wind  components  was  significant,  and  2)  the  effects  of  both  elevation  and 
aspect  were  significant  on  the  potential  temperature.  Of  the  four  case  studies  pres¬ 
ented,  the  two  daytime  cases  are  of  interest  here.  The  first  case  illustrated  a  typical 
"light  wind"  situation.  The  observations  were  marked  by  a  pronounced  upslope  flow 
at  all  stations.  The  lapse  rate  indicated  neutral  to  slightly  stable  air  over  the  southern 
slopes  and  significant  instability  over  the  northern  slopes.  Temperature  observations 
very  near  the  surface  on  the  southwest  slope  indicated  a  very  strong  positive  heat  flux 
on  that  day.  The  second  case  illustrated  a  typical  "strong  wind"  situation.  The  strong 
ambient  flow  dominated  the  (4>servations  at  all  stations.  The  authors  found  they  were 
able  to  recover  an  upslope  flow  similar  to  the  light  wind  case  by  subtracting  a  fraction 
of  the  mountain-top  wind  from  the  wind  velocities  everywhere.  This  implies  that  the 
strong  ambient  winds  simply  masked  the  weaker  thermally-driven  upslope  flow.  The 
lapse  rates  in  this  case  indicated  nearly  neutral  conditions  everywhere.  The  near- 
surface  temperatures  on  the  southwest  slope  indicated  a  much  smaller  heat  flux  in  this 
case  than  in  the  light  wind  case.  The  authors  speculated  that  the  strong  winds  pre- 
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vented  the  buildup  of  sufficient  near-surface  temperature  gradients  to  produce  the 
positive  heat  flux  magnitude  found  in  the  light  wind  case. 

2.2  Numerical  ModeUi^  Stutfies 

Many  numerical  models  have  been  used  to  simulate  mountain  convectitm  and  to 
improve  our  understanding  of  some  of  the  above  observations.  Early  numerical  stu¬ 
dies  were  almost  exclusively  two-dimensional  and  contained  only  crude,  if  any,  para- 
meterizations  of  subgrid-scale  processes,  such  as  turbulence  and  cloud  micn^ysics. 
Later  models  expanded  to  three  dimensions  and  included  more  sophisticated  treat¬ 
ments  of  subgrid-scale  processes.  A  brief  summary  of  some  of  these  studies  is  given 
below. 

Orville  (1964)  developed  a  two-dimensional  model  based  on  the  equations  of  Ogu- 
ra  (1962)  to  study  mountain  upslope  winds.  In  this  study  Orville  simulated  the  devel¬ 
opment  of  an  upslope  circulation,  the  center  of  which  moved  from  the  lower  slope  of 
the  mountain  to  well  above  the  peak.  Several  elements  of  the  simulated  wind  field 
matched  those  described  by  Defant  (1951).  The  thermal  fields  from  this  simulation 
also  showed  some  similarities  to  those  described  by  Braham  and  Draginis  (1960). 

In  a  later  paper,  Orville  (1965a)  reported  the  results  of  simulations  that  included 
the  effects  of  moisture  in  his  two-dimensional  model.  Comparison  of  these  results 
with  those  from  the  "dry"  simulations  (Orville,  1964)  showed  that  moisture  effects 
caused  the  upslope  flow  to  develq)  10  to  20  percent  filter  than  in  the  dry  model. 
Again,  the  model  showed  some  qualitative  agreement  with  the  observations  of  Braham 
and  Draginis  (1960).  Other  simulations  in  this  study  were  compared  to  the  photo- 
grammetric  observations  reported  by  Orville  (l%5b).  The  simulations  produced 
clouds  in  a  reasonable  time  scale  and  at  reasonable  elevations.  Orville  pointed  out, 
however,  that  neither  the  fields  of  moisture  and  temperature  nor  the  time  scale  of 
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cloud  growth  were  realistic.  One  simulation  in  this  study  produced  a  pulsation  remi¬ 
niscent  of  those  observed  by  Todd  (1964)  and  Anderson  (1960). 

Orville  (1968)  extended  his  previous  studies  to  include  the  effects  of  ambient 
winds  (M)  the  initi^on  and  development  of  simulated  nnountain  cumulus  clouds.  In 
general,  he  found  that  vertical  wind  shear  enhanced  cloud  initiation  but  damped  fur¬ 
ther  cloud  development.  In  addititm  the  simulations  showed  that  the  clouds  formed 
downwind  of  the  peak  as  they  had  in  most  of  the  observational  studies.  Orville  also 
discussed  the  simulation  of  a  mountain  wave  and  its  associated  cloud  development. 
He  found  that  in  this  situation,  the  damming  effect  of  the  mountain  and  heating  of  the 
upwind  slope  acted  together  to  produce  a  wave  in  the  upwind  airflow.  When  small 
perturbations  of  moisture  and  potential  temperature  reached  this  wave,  they  were  am¬ 
plified  and  advected  over  the  top  of  the  mountain  into  air  moistened  by  the  large 
downwind  circulation.  At  this  point,  the  amplified  perturbations  initiated  clouds.  Or¬ 
ville  felt  this  method  could  explain  the  pulsations  observed  in  stationary  clouds  that 
form  in  the  downwind  circulation. 

Liu  and  Orville  (1969)  added  the  effects  of  precipitation  and  cloud  shadows  to  Or¬ 
ville's  model.  The  precipitation  processes  were  parameterized  following  Srivastava 
(1967)  and  Kessler  (1969),  and  the  cloud  shadow  effects  were  simulated  by  decreas¬ 
ing  the  temperature  and  water  vapor  changes  at  the  ground  points  directly  below  the 
cloud.  Comparison  of  precipitating  and  non-precipitating  clouds  showed  the  effect  of 
the  precipitation  processes  to  be  negligible  during  early  development.  During  later 
stages  of  development,  however,  evaporation  below  the  cloud  in  the  precipitating  case 
generated  a  downdraft  that  tended  to  shorten  the  cloud's  lifetime.  The  shadow  effects 
caused  both  the  clouds  to  move  out  of  the  grid  faster  and  subsequent  clouds  to  be 
smaller. 
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Gal-Chen  and  Somerville  (197Sa)  developed  a  terrain-following  coordinate  system 
for  use  in  simulations  with  irregular  lower  boundary  crxiditions.  This  coordinate 
transformation  greatly  simplified  the  lower  boundary  conditions  but  only  slightly 
complicated  the  basic  model  equations.  The  complications  take  the  form  of  terrain- 
dependent  transformation  coefficients,  which  must  be  included  in  some  terms  of  the 
model  equations.  Details  of  this  coordinate  transformation  will  be  discussed  further 
in  Chapter  4.  Gal-Chen  and  Somerville  (197Sb)  used  a  two-dimensional  model  based 
on  the  Navier-Stokes  equations  to  test  their  new  coordinate  transformation.  They 
simulated  upslope  mountain  winds  and  found  their  results  to  be  similar  to  Orville's 
(1964)  observations. 

In  the  late  1970' s  numerous  three-dimensional,  non-hydrostatic  tiKxlels  were  de¬ 
veloped  (see  for  example  Tapp  and  White,  1976;  Clark,  1977;  Schlesinger,  1975, 
1978;  Cotton  and  Tripoli,  1978;  Klemp  and  Wilhelmson,  1978a, b).  Since  that  time 
these  models  have  been  used  mainly  to  study  severe  storms.  In  these  studies  the  con¬ 
vection  was  usually  started  by  the  imposition  of  a  temperature  and/or  a  moisture  per¬ 
turbation  on  an  initially  horizontally  homogenous  environment.  Only  a  few 
investigators  have  used  three-dimensional  models  to  study  mountain  convection. 
These  studies  will  be  summarized  next. 

Clark  and  Gall  (1982)  used  a  three-dimensional  model  to  simulate  the  airflow  over 
mountainous  terrain.  One  part  of  their  study  attempted  to  simulate  the  case  of  dry 
convection  over  Mount  Withington  in  New  Mexico  presented  by  Raymond  and  Wil- 
kening  (1980).  A  prominent  feature  of  the  simulations  was  the  development  of  strong 
longitudinal  rolls  after  sunrise.  This  longitudinal  development  did  not  compare  well 
with  the  more  cellular  development  found  in  the  observations.  The  authors  attributed 
this  difference  to  the  crude  surface  layer  parameterization  used  in  the  model,  to  the 


14 


difficulty  in  choosing  low  level  winds,  and  to  poor  horizontal  grid  resolution.  The 
simulated  flow  field  also  failed  to  show  the  low-level  convergence  and  higher  level 
divergence  described  by  Raymond  and  Wilkening. 

Smolarldewicz  and  Clark  (1985)  developed  a  surface  boundary  layer  model  which 
used  single-level  surface  mesonetwork  data  to  calculate  surface  fluxes  of  momentum, 
heat,  and  moisture.  These  fluxes  were  then  used  in  a  three-dimensional  model  to 
simulate  the  evolution  of  a  field  of  cumulus  clouds  over  complex  terrain  in  western 
Montana.  The  simulated  cloud  field  qualitatively  agreed  quite  well  with  photographs 
and  other  observations  of  the  actual  clouds.  Sensitivity  experiments  showed  that  dur¬ 
ing  the  early  stages  of  cloud  development,  the  dynamical  effects  of  flow  over  the  ter¬ 
rain  and  the  thermodynamical  effects  of  soil  type  and  vegetation  played  equal  roles  in 
the  rate  of  development  of  the  cloud  field.  During  the  later  stages  of  development, 
the  thermodynamical  effects  primarily  influenced  the  local  properties  of  the  cloud 
field.  The  authors  found  that  the  gross  characteristics  of  the  cloud  field  were  mostly 
determined  by  the  topography  rather  than  the  thermodynamics  of  the  surface  layer. 

Banta  (1983,  1986)  used  a  two-dimensional  version  of  the  Cotton-Tripoli  model  to 
study  the  boundary  layer  evolution  over  mountainous  terrain.  His  numerical  experi¬ 
ments  showed  the  importance  of  the  inversion  layer  and  the  ridgetop  winds  to  the 
formation  and  duration  of  the  upslope  flow.  He  found  that  the  duration  of  the  upslope 
flow  was  inversely  proportional  to  the  strength  of  the  ridgetop  winds.  This  finding 
implied  that  the  leeside  convergence  effect  discussed  earlier  is  less  effective  when  the 
upper-level  winds  are  strong. 


Chapters 

Observations  Of  The  31  July  1984  Storm  And  Its  Environment 
This  chapter  describes  the  observed  initiation  and  evolution  of  a  mountain  thiui- 
derstorm  over  the  Magdalena  Mountains  near  Socorro,  New  Mexico  on  31  July  1984. 
Section  3.1  describes  the  general  nature  of  central  New  Mexico  mountain  storms. 
Section  3.2  describes  the  field  experiment  facilities,  and  the  data  collection  and  analy¬ 
sis  methods  used  to  acquire  and  analyze  the  data  used  in  this  research.  The  local  ter¬ 
rain  in  the  field  experiment  area  also  is  described  in  this  section.  Finally,  section  3.3 
details  the  31  July  1984  storm  history  as  observed  by  the  Doppler  radars. 

3.1  New  Mexico  Mountain  Thunderstorms 

The  small,  isolated  summer  thunderstorms  that  frequently  develop  over  the  Mag¬ 
dalena  Mountains  and  other  nearby  mountain  ranges  of  New  Mexico  have  been  stu¬ 
died  extensively  for  at  least  the  past  15  years  (see  for  example  Winn  et  al.,  1974; 
Raymond  and  Wilkening,  1980,  1982,  1985;  Dye  etal.,  1987;  and  Lang,  1991).  The 
storms  usually  occur  under  conditions  of  weak  environmental  winds  and  wind  shear 
(Raymond  and  Wilkening,  1985),  and  so  generally  are  thought  to  develop  primarily  as 
a  result  of  the  thermally  driven  mountain-valley  circulations  described  in  Chapter  2 
(Raymond  and  Wilkening,  1980,  1982;  and  Dye  era/.,  1987).  Orographic  lifting  and 
flow  into  convergent  valleys  resulting  from  the  weak  environmental  winds  may  play  a 
secondary  role  in  storm  development  and  a  more  important  role  in  determining  the 
location  of  development  (see,  for  example  Barker  and  Banta,  1985;  Banta  and  Schaaf, 
1987;  and  Schaaf  et  al. ,  1988).  Initial  cloud  development  is  often  restricted  weak 
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stable  layers  which  are  eroded  slowly  by  the  convective  action  until  small  thunder¬ 
storms  develop  with  tops  ranging  from  9  to  14  km  MSL  (Dye  et  al.,  1987),  As  a  re¬ 
sult  of  the  weak  environmental  winds,  the  storm  movements  are  small,  making  the 
Magdalena  Mountain  range  and  surrounding  ranges  ideal  locations  for  thunderstorm 
studies. 

3.2  Data  Collection  And  Analysis 

As  part  of  the  Thunderstorm  Research  International  Project  (TRIP),  a  field  experi¬ 
ment  was  conducted  near  the  Langmuir  Laboratory  (elevation  3255  m  MSL)  in  the 
Magdalena  Mountains  of  central  New  Mexico  during  the  summer  of  1984.  During 
the  field  experiment,  22  storms  occurred  on  19  different  days  from  14  July  to  24  Au¬ 
gust  (Dye  et  al.,  1987).  This  study  uses  primarily  the  data  collected  on  31  July  1984. 
The  numerical  simulations  conducted  for  this  study  were  initialized  with  rawinsonde 
and  Doppler  radar-derived  wind  data  from  the  31  July  case.  The  simulations  were 
then  compared  with  the  31  July  Doppler  radar  observations  and  with  Lang's  (1991) 
study  of  the  microphysical  evolution  of  the  3  August  storm. 

3.2.1  Facilities 

Four  Doppler  radars,  two  instrumented  research  aircraft,  rawinsondes,  and  a  sur¬ 
face  network  for  measuring  electrical  fields  were  employed  to  collect  storm  data  over 
the  Mi^dalena  Moimtains  during  the  TRIP  experiment.  Rawinsonde  observations 
were  collected  from  a  site  (elevation  3223  m  MSL)  in  the  mountains  near  the  Lang¬ 
muir  Laboratory  and  from  the  Socorro  Airport  (elevation  1456  m  MSL)  in  the  valley 
during  the  mornings  and  early  afternoons  throughout  the  experiment.  The  radar  net¬ 
work  (Fig.  1)  consisted  of  two  National  Center  for  Atmospheric  Research  (NCAR) 
5-cm  radars,  CP-3  and  CP-4,  and  two  National  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA)  3-cm  radars,  NOAA-C  and  NOAA-D.  Lang  (1991)  and  Dye  et  al. 
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LANGMUIR  TERRAIN 
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Fig.  1.  Radar  network  and  topographical  view  of  terrain  features  near  Langmuir  Labo¬ 
ratory.  Terrain  elevation  is  relative  to  mean  sea  level  in  kilometers  and  contour  inter¬ 
vals  are  every  0.2  km.  The  horizontal  and  vertical  axes  are  labeled  in  kilometers.  The 
origin  (0,0)  is  located  at  Langmuir  Laboratoty.  A  sample  radar  analysis  domain  is 
shown  for  reference. 
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(1987)  described  in  detail  the  location  and  characteristics  of  the  Doppler  radars.  De¬ 
pending  on  storm  development,  the  radars  collected  data  from  mid-morning  into  the 
afternoon  throughout  the  experiment. 

3.2.2  Rawinsonde  Observations 

Two  thermodynamic  soundings  were  collected  from  rawinsondes  launched  on  the 
morning  of  31  July  1984.  The  first  sonde  was  launched  at  0719  MST  from  the  Socor¬ 
ro  airport  located  in  the  valley  several  miles  east  of  the  Magdalena  Mountains.  The 
second  was  launched  at  0945  MST  from  the  balloon  hangar  near  the  Langmuir  Labo¬ 
ratory.  This  sounding  coincided  with  the  onset  of  early  convective  development  as 
observed  by  the  NCAR  CP-4  radar. 

Figures  2  and  3  show  the  airport  and  Langmuir  thermodynamic  soundings,  respec¬ 
tively.  Both  soundings  exhibit  essentially  the  same  characteristics.  The  soundings 
show  a  conditionally  unstable  atmosphere  with  a  mixed  layer  extending  from  the  sur¬ 
face  to  SSO  mb  (500  mb  in  the  Langmuir  sounding).  The  mixed  layer  was  relatively 
moist  throughout  with  a  very  moist  layer  30-40  mb  deep  near  its  top.  The  mixed  lay¬ 
er  and  the  associated  moist  layer  aloft  were  capped  by  an  inversion  approximately 
50-100  mb  deep.  Below  300  mb,  the  Langmuir  Laboratory  sounding  tended  to  be 
more  moist  and  slightly  cooler  than  the  airport  sounding.  The  thermodynamic  struc¬ 
ture  observed  in  both  31  July  1984  soundings  was  very  similar  to  that  reported  by 
Lang  (1991)  for  the  3  August  1984  case  and  by  Raymond  and  Wilkening  (1985)  for 
their  "composite  thunderstorm  sounding”.  The  observed  thermodynamic  profiles 
were  combined  and  smoothed  as  described  in  Section  5.2.1  to  initialize  the  numerical 
simulations. 

The  wind  data  (not  shown)  from  both  the  airport  and  Langmuir  rawinsondes  are 
very  noisy  and  dissimilar  due  to  problems  with  the  rawinsonde  tracking  equipment 
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SOCORRO  0719  MST 
31  JUL  84 


TEMPERATURE  (OQ 


Fig.  2.  Skew  T,  log  P  plot  of  the  0719  MST  31  July  1984  Socorro  Airport  thermody¬ 
namic  sounding.  Solid  skewed  lines  are  temperatures  in  degrees  Celsius;  dashed 
skewed  lines  are  mixing  ratios  in  grams  per  kilogram;  curved  dashed  lines  are  dty 
adiahats  in  degrees  Celsius;  solid  horizontal  lines  are  pressures  in  millibars.  Heavy 
solid  and  dashed  lines  represent  sensible  and  dewpoint  temperature  profiles,  respective¬ 
ly- 


PRESSURE  (mb) 


20 


LANQMUIR  0945  M8T 
31  JUL  84 


Fig.  3.  As  in  Fig.  2,  0945  MST  31  July  1984  Langmuir  Laboratory  thermodynamic- 
sounding. 
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(Winn,  personal  communication,  1987;  Lang,  1991).  The  tracking  problem  persisted 
throughout  the  field  experiment. 

To  overcome  the  problems  with  tlM  rawinsonde  wind  data,  alternative  wind  pro¬ 
files  for  heights  above  the  mountain  top  were  constructed  from  the  synthesized 
pier  radar  wind  data  (see  Section  3.2.4).  The  only  wind  information  available  (other 
than  rawinsonde)  for  below  the  height  of  the  mountain  top  was  an  estimated  profile 
provided  by  a  TRIP  scientist  who  participated  in  the  observational  phase  of  the  ex¬ 
periment  (Ziegler,  personal  communication,  1986).  The  profile  was  estimated  by  us¬ 
ing  Doppler  radar  wind  data  above  the  mountain  top  and  extending  the  profile  down 
to  the  valley  floor  using  an  Ekman  spiral  approximation.  Figure  4  shows  this  profile. 
These  alternative  Doppler-derived  and  estimated  wind  profiles  were  combined  and 
smoothed  as  described  in  Section  S.2.2  to  initialize  the  numerical  simulations. 

3.2.3  Doppler  Radar  Observations  And  Analyses 

All  four  Doppler  radars  in  the  network  surrounding  Langmuir  Laboratory  were 
operational  on  31  July  1984.  The  radars  collected  radial  velocity  and  reflectivity  data 
intermittently  between  0946  and  1058  MST,  and  collected  data  regularly  between 
1058  and  1219  MST  with  one  break  between  1119  and  1134  MST.  The  radial  veloc¬ 
ity  data  were  used  to  synthesize  three-dimensional  wind  fields  at  the  same,  analysis 
times  as  the  reflectivity  analyses.  The  wind  and  reflectivity  analyses  then  were  used 
to  investigate  the  31  July  1984  storm  development  and  to  evaluate  the  numerical 
simulations  of  the  storm  development. 

The  multiple-Ekjppler  wind  analyses  used  in  this  study  were  performed  by  Lang 
(personal  communication,  1992)  and  follow  exactly  his  implementation  of  Ray  etal.'s 
(1980)  iterative  method  (Lang,  1991)  for  the  3  August  case.  The  method  determines 
the  wind  components  at  each  analysis  grid  point  by  solving  the  set  of  equations  con- 
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ZIEGLER  ESTIMATED  WIND  PROFILE 


WIND  SPEED  (M/S) 


Fig.  4.  Vertical  profiles  of  estimated  u  (solid  curve)  and  v  (dashed  curve)  wind  compo¬ 
nents  from  Ziegler,  personal  communication  (1986).  The  horizontal  axis  is  wind  speed 
in  meters  per  second,  and  pressure  in  millibars  is  along  the  vertical  axis. 


23 


sisting  of  the  overdetermined  solutions  for  the  u  and  v  wind  components  (derived  by 
using  the  radial  velocity  relationship  for  each  radar),  Joss  and  Waldvogel's  (1970) 
particle  terminal  velocity  relationship,  and  the  anelastic  continuity  equation.  The  it¬ 
eration  first  solves  the  u  and  v  equations  using  a  first  guess  for  w,  which  is  the  pre¬ 
viously  calculated  w  for  the  layer  immediately  above  the  analysis  grid  point.  At  the 
top  of  the  analysis  grid,  this  value  is  a  specified  kinematic  uj^r  boundary  condition 
for  w.  The  continuity  equation  is  next  integrated  downward  to  give  a  value  for  w  at 
the  bottom  of  the  analysis  grid  box.  A  mean  grid  box  w  is  calculated  using  the  new 
values  of  w  at  the  top  and  bottom  of  the  grid  box,  and  u  and  v  are  recalculated.  The 
iteration  continues  until  the  solution  converges.  As  in  Lang's  (1991)  3  August  analy¬ 
sis,  a  vertical  velocity  of  2.0  m  s'‘  was  imposed  as  the  upper  boundary  condition  over 
the  entire  analysis  domain.  Unlike  the  3  August  analysis  which  only  used  radial  ve¬ 
locity  data  from  three  radars,  the  wind  analyses  for  the  31  July  case  incorporated  ra¬ 
dial  velocity  data  from  all  four  radars.  Lang  (1991)  showed  Doppler-derived  vertical 
velocities  for  the  3  August  case  compared  favorably  with  observed  vertical  velocities 
collected  from  an  instrumented  aircraft. 

The  Doppler  radar  reflectivity  and  wind  analyses  were  performed  at  3-minute  in¬ 
tervals  starting  at  0946  MST.  As  in  Lang's  (1991)  3  August  analysis,  the  31  July 
wind  analysis  was  performed  on  a  three-dimensional  grid  with  30  gridpoints  on  each 
side  and  20  points  in  the  vertical.  The  grid  spacing  was  one-half  kilometer  in  all 
directions.  During  the  period  of  greatest  storm  development,  the  grid  was  centered 
on  Langmuir  Laboratory.  At  several  times  early  in  the  analysis  period,  the  analysis 
grid  was  centered  south  of  Langmuir  Laboratory  to  capture  early  convective  develop¬ 
ment  in  that  area.  Throughout  the  entire  analysis  period,  the  lowest  analysis  level  was 
3.3  km  MSL,  putting  it  very  near  the  elevation  of  the  mountain  tops  near  Langmuir. 


24 


This  study  follows  Lang  (1991)  and  takes  the  analyzed  reflectivities  directly  from 
the  NCAR  CP-4  observations.  As  with  the  3  August  1984  case  (Lang,  1991),  the 
CP-4  radar,  located  on  the  plain  approximately  22  km  to  the  northwest  of  Langmuir, 
afforded  the  best  view  of  the  31  July  1984  storm  develt^ment.  Additionally,  as  noted 
by  Lang  (1991),  the  S-cm  wavelength  CP-4  was  less  likely  to  be  attenuated  than  the 
3-cm  NOAA  radars. 

3.2.4  Doppler-Derived  Environnieiital  Wind  Profiles 

Section  3.2.2  described  the  problems  associated  with  both  rawinsonde  wind  pro¬ 
files  collected  on  31  July  1984.  To  avoid  these  problems  and  still  achieve  an 
observation-based  estimate  of  the  pre-storm  environmental  wind  field,  a  scheme  was 
developed  to  estimate  the  wind  profile  using  Doppler  radar  data. 

The  Doppler-derived  wind  profiles  were  constructed  with  a  scheme  which  uses  the 
observed  Doppler  radar  reflectivity  and  synthesized  wind  data  described  in  Section 
3.2.3.  An  inherent  problem  with  using  Doppler  wind  data  to  estimate  pre-storm  envi¬ 
ronmental  winds  is  that  the  data  are  available  only  after  convection  develops  enough 
to  produce  a  radar  echo.  Thus,  the  Doppler-derived  wind  profile  may  be  strongly  in¬ 
fluenced  by  convective  circulations  and  may  not  necessarily  represent  the  pre-storm 
wind  profile.  However,  by  assuming  the  wind  field  near  the  storm  boundary  is  not 
significantly  affected  by  the  storm  circulations,  the  pre-storm  wind  profile  may  be 
reasonably  estimated  by  using  only  the  wind  inforroati<m  from  the  outer  boundary  of 
the  radar  echo.  The  scheme  used  in  this  study  follows  this  approach. 

As  described  in  Section  3.2.3,  Doppler  radar  reflectivity  and  synthesized  wind 
data  were  available  from  0946  through  1219  MST  31  July  1984.  The  data  were  avail¬ 
able  above  the  mountain  top  (approximately  3.3  km)  on  horizontal  cross-sections 
spaced  O.S  km  ^xut  vertically.  Very  little  wind  data  could  be  synthesized  for  the 
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times  between  1019  and  1 1 16  MST.  For  the  purpose  of  creating  environmental  wind 
profiles,  the  radar  reflectivity  and  synthesize  wind  data  were  divided  into  two  sec¬ 
tions;  the  first  section  included  data  collected  between  0946  and  1019  MST,  and  the 
second  section  included  data  collected  between  1116  and  1216  MST.  For  simplicity, 
the  "cloud  boundary"  was  define  as  any  region  with  a  radar  reflectivity  of  less  than 
IS  dBZ.  For  each  radar  analysis  time,  the  mean  wind  within  the  clcHid  boundary  was 
compute  at  each  horizontal  level.  Figures  5  and  6  show  the  resulting  Doppler- 
derive  wind  profiles  after  interpoIatitM)  to  pressure  levels  for  the  data  peries  begin¬ 
ning  at  0946  MST  and  1116  MST,  respectively.  These  profiles  were  combine  as 
describe  in  Section  5.2.2  to  construct  the  base-state  wuid  profiles  use  to  initialize 
the  numerical  simulations. 

3.2.5  Local  Terrain 

Langmuir  Laboratory  is  locate  on  South  Baldy  Peak  at  an  elevation  of  3255  m  in 
the  center  of  the  Magdalena  Mountains.  Two  ridges  extend  south  and  southeast 
(denote  R1  and  R2  respectively)  ^m  South  Baldy  Peak.  A  valley  (denote  VI) 
separates  R1  and  R2.  Two  minor  ridges  extend  eastward  from  R2,  one  to  the  south¬ 
east  (denote  R3)  and  the  other  to  the  northeast  (denote  R4).  A  small  valley 
(denote  V2)  separates  R2  and  R3.  A  minor  valley  (denote  V3)  separates  R3  from 
R4.  Finally,  a  deep  valley  (denote  V4)  separates  R4  ftom  the  main  ridge  that  con¬ 
tains  South  Baldy  Peak.  Figure  7  shows  the  above  features  on  a  topographic  map  of 
the  terrain  surrounding  Langmuir  Laboratory. 

3.3  31  July  1984  Storm  History:  Radar  Observations 

This  section  provides  a  description  of  the  31  July  1984  storm  development  using 
the  Doppler  radar  reflectivity  and  synthesize  wind  analys.'^s  describe  in  previous 
sections.  The  detaiie  analyses  which  follow  track  the  storm  development  from  weak 
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EARLY  DOPPLER  WIND  PROFILE 


Fig.  5.  Vertical  profile  of  Doppler-derived  u  (solid  curve)  and  v  (dashed  curve)  wind 
components  calculated  as  described  in  text  using  Doppler-derived  wind  data  from  the 
period  between  0946  and  1019  MST.  Axes  labeled  as  in  Fig.  4. 
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LATE  DOPPLER  WIND  PROFILE 


Fig.  6.  Vertical  profile  of  Doppler-derived  u  (solid  curve)  and  v  (dashed  curve)  wind 
components  calculated  as  described  in  text  using  Doppler-derived  wind  data  from  the 
period  between  1116  and  1216  MST.  Axes  labeled  as  in  Fig.  4. 
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UNSMOOTHED  MODEL  TERRAIN 


Fig.  7.  Topographical  view  of  terrain  features  near  Langmuir  Laboratory.  Terrain 
elevation  is  relative  to  sea  level  in  dekameters  and  contour  intervals  are  every  21)  deka- 
meters.  The  horizontal  and  vertical  axes  are  labeled  in  kilometers.  The  origin  (0.0)  is 
located  at  Langmuir  Laboratory.  Valleys  and  ridges  are  labeled  as  described  in  text. 
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initial  activity  in  the  mid-ntoming  hours  to  full  development  shortly  before  noon  to 
gradual  decay  shortly  after  noon.  All  (x,,v)  coordinate  references  are  in  kilometers 
and  based  on  a  coordinate  system  with  its  origin  at  Langmuir  Laboratory. 

3.3.1  Pre-Storm  Developmeiit 

The  NCAR  CP-4  radar  first  detected  an  echo  at  0946  MST  over  the  eastern  slope 
of  the  main  ridge  approximately  l.S  km  north  of  Langmuir  Laboratory.  The  first 
echo  extended  vertically  to  only  6.3  km  MSL.  The  radar  detected  several  shallow, 
weak,  short-lived  showers  between  0946  and  1110  MST,  During  this  period  one 
shallow  shower  persisted  north  of  Langmuir  Laboratory.  The  shower's  radar  returns 
rarely  exceeded  40  dBZ  and  never  extended  above  6.3  km  MSL.  Although  these  ear¬ 
ly  showers  never  developed  into  thunderstorms,  they  served  to  moisten  the  low-level 
environment  for  later  storm  development. 

3.3.2  Storm  Initiatioii  And  Early  Development 

By  1116  MST  the  CP-4  radar  detected  an  organized  cell  (denoted  hereafter  as  Cl) 
(-1 .0,  -4.5)  with  returns  of  approximately  20  dBZ  over  and  east  of  R1  south  of  Lang¬ 
muir  Laboratory.  Cl  was  the  first  in  a  series  of  cells  that  developed  and  combined  to 
form  a  multi-cellular  storm.  For  the  purposes  of  this  study  a  cell  is  defined,  as  in 
Lang's  (1991)  analysis,  as  any  relative  reflectivity  and/or  updraft  maximum  that  ex¬ 
hibits  vertical  continuity.  The  radar  also  detected  weak  returns  over  R2  southeast  of 
the  Laboratory.  The  radar  returns  at  elevation  3.3  km  MSL  located  along  the  main 
ridge  just  to  the  north  of  Langmuir  Laboratory  appear  to  be  ground  clutter  and  per¬ 
sisted  through  most  of  the  remaining  analysis  times.  By  1119  MST,  Cl's  radar  re¬ 
turns  (-1.0,  -4.5)  reached  30  dBZ,  and  20  dBZ  echoes  extended  to  7.3  km  MSL.  The 
activity  over  R2  southeast  of  Langmuir  continued  to  develop  slowly  (Fig.  8). 


30 


07/31/814 

1119  MST  3.3  KM  MSL 

tl  RADARS  — *1S.0H/S 


HORIZONTAL  FLOW  FIELD/REFLECTIVITT  (DBZ) 


Fig.  H.  Horizontal  cross-section  of  observed  radar  reflectivities  in  dBZ  and  vector  velo¬ 
cities  at  an  altitude  of  3.3  km  relative  to  mean  sea  level  at  1119  MST  31  July  1984.  Re¬ 
flectivities  (solid  contours)  are  contoured  every  10  dBZ.  Vector  velocities  are 
ground-relative  with  a  vector  length  of  one  grid  intern  al  (0.5  km)  equal  to  15  m  s  '. 
Terrain  elevations  (dashed  contours)  are  relative  to  sea  level  in  meters,  and  contour 
intervals  are  every  200  meters.  The  origin  is  at  Langmuir  Laboratory,  and  the  axes  are 
labeled  in  kilometers. 
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At  the  next  available  analysis  time,  1134  MST,  the  radar  detected  two  well- 
defined  cells  over  the  ridges:  Cl  (0.0,  -5.0)  over  R1  south  of  Langmuir;  and  a  new 
cell  (denoted  hereafter  as  C2)  (2.5,  -1.0)  over  R2  approximately  3  km  east-southeast 
of  Langmuir.  Cl's  radar  reflectivities  exceeded  40  dBZ  and  C2's  exceeded  45  dBZ. 
Cl's  20  dBZ  returns  extended  to  9.3  km  MSL,  while  C2's  extended  to  10.0  km.  The 
horizontal  wind  analysis  at  3.3  km  MSL  showed  significant  upslope  flow  on  the  west¬ 
ern  side  of  R1  along  with  upslope  flow  through  VI,  The  two  flows  converged  very 
near  the  location  of  Cl.  Additionally,  upslope  flow  on  the  eastern  side  of  R2  con¬ 
verged  with  the  VI  flow  in  the  vicinity  of  C2.  Above  6.3  km  MSL,  the  horizontal 
wind  analysis  also  showed  strong  divergent  flow  over  Cl  and  weak  divergence  over 
C2  (Fig.  9).  Cl  contained  a  1.5  km-wide  updraft  located  on  the  southeastern  edge  of 
the  cell  and  an  approximately  1.0  km-wide  downdraft  located  just  to  the  north  of  the 
updraft.  The  updraft  exceeded  9.0  m  s  *  between  3.3  and  6.3  km  MSL  with  a  maxi¬ 
mum  updraft  of  12.0  m  s’*  located  at  5.3  km  MSL.  The  updraft  diminished  rapidly 
above  7.3  km  MSL.  The  downdraft  averaged  approximately  3.0  m  s  *  between  3.3 
and  4.3  km  MSL.  C2  also  contained  a  1.5  km-wide  updraft  located  on  the  southeast¬ 
ern  edge  of  the  cell  and  an  approximately  1.0  km-wide  downdraft  located  north  of  the 
updraft.  The  updraft  exceeded  9.0  m  s’*  between  5.3  and  7.3  km  MSL  and  remained 
identifiable  up  to  10.3  km  MSL.  Like  Cl's  downdraft,  C2's  extended  between  3.3 
and  4.3  km  MSL. 

3.3.3  Mature  Storm  Development 

By  1137  MST,  both  Cl's  and  C2's  maximum  reflectivities  exceeded  50  dBZ  at 

3.3  km  MSL.  The  CP-4  radar  again  detected  20  dBZ  returns  up  to  9.3  km  MSL  from 
Cl  (0.0,  -4.5)  and  to  10.0  km  MSL  from  C2  (2.5,  -1.5).  The  horizontal  wind  analy¬ 
sis  at  3.3  km  MSL  showed  little  change  from  the  1134  MST  analysis.  The  strong  di- 
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Fig.  9a.  Horizontal  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
as  in  Fig.  S.  except  at  1134  MST  and  at  3.3  km  MSL. 
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Fig.  9b.  As  in  Fig.  8,  horizontal  cross-section  of  observed  radar  reJJf'ctivities  and  vec¬ 
tor  velocities  for  1 134  MST  and  at  7.3  km  MSL. 
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vergence  above  6.3  km  MSL  continued  over  Cl  while  virtually  non-divergent  flow 
existed  over  C2.  Cl's  updraft  expanded  to  a  width  of  approximately  3.0  km.  The 
enlarged  updraft  contained  two  areas  of  maximum  vertical  velocities  at  low  levels 
which  merged  into  a  single  updraft  at  6.3  km.  Cl's  updraft  exceeded  9.0  m  s'*  be¬ 
tween  4.3  and  7.3  km  MSL  with  a  maximum  of  12.0  m  s'*  between  5.3  and  6.3  km 
MSL.  Cl's  downdraft  became  extremely  weak  (less  than  3.0  m  s'*).  C2's  updraft 
began  to  weaken  considerably  from  the  bottom  with  vertical  velocities  exceeding  9.0 
m  s  *  only  between  7.3  and  8.3  km  MSL.  C2's  downdraft  became  much  stronger 
than  at  1134  MST,  exceeding  3.0  m  s'*  between  3.3  and  6.3  km  MSL  with  a  maxi¬ 
mum  of  6.0  m  s  *  at  3.3  km  MSL.  Strong  divergent  flow  near  7.0  km  MSL  over  Cl 
contributed  to  the  weakening  of  C2's  updraft. 

The  trend  of  Cl  strengthening  and  C2  weakening  continued  at  1140  MST.  The 
radar  continued  to  detect  maximum  returns  of  50  dBZ  from  both  cells,  but  Cl's  stron¬ 
gest  returns  (-0.25,  -4.75)  extended  higher  (to  4.3  km  MSL)  than  at  the  previous 
time.  Cl's  double-maxima  updraft  continued  with  the  northern  branch  strengthening 
to  15.0  m  s'*  between  5.3  and  6.3  km  MSL  while  the  southern  branch  remained  at 
about  9.0  m  s'*.  The  two  branches  merged  near  8.3  km  MSL.  Meanwhile,  C2's 
(2.25,  -0.75)  updraft  weakened  to  about  6.0  m  s'*.  Both  cells  contained  weak  down- 
drafts  at  this  analysis  time.  Cl's  downdraft  remained  less  than  3.0  m  s'*,  while  C2's 
decreased  to  about  3.0  m  s'*  between  3.3  and  4.3  km  MSL.  The  vertical  cross- 
section  shown  in  Figure  10  again  shows  strong  divergent  flow  near  7.0  km  MSL  over 
Cl  undercutting  the  remnants  of  C2's  updraft. 

By  1143  MST,  both  cells  had  weakened  from  their  previous  strengths.  Cl's 
maximum  reflectivity  (-0.25,  -4.0)  decreased  to  45  dBZ  and  C2's  (2.25,  -1.75)  de¬ 
creased  to  40  dBZ.  The  horizontal  flow  at  3.3  km  MSL  remained  virtually  un- 
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Fig.  lOa.  Horizontal  cross-section  of  observed  radar  reflectivities  in  dBZ  at  an  altitude 
of  4.3  km  relative  to  mean  sea  level  at  1140  MST  31  July  1984.  Reflectivities  are  con¬ 
toured  every  5  dBZ.  The  origin  Ls  at  Langmuir  Laboratory.  The  axes  are  labeled  in 
kilometers.  Line  AB  shows  the  location  of  the  vertical  cross-section  in  Fig.  I  Oh. 
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Fig.  I  Oh.  Vertical  cross-section  of  observed  radar  r^ectivities  in  dBZ  and  vector  velo¬ 
cities  at  1140  MST  31  July  1984.  Reflectivities  are  contoured  every  5  dBZ.  Vector  ve¬ 
locities  are  ground-relative  with  the  vector  length  shown  in  upper  right-hand  corner  of 
figure  equal  to  20  m  .s  '.  Cross-section  location  AB  is  shown  in  Fig.  10a. 
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changed  from  previous  times.  The  horizontal  wind  analyses  between  4.3  and  S.3  km 
MSL,  however,  showed  that  the  circulation  induced  by  Cl  significantly  diminished 
the  southwesterly  flow  into  C2  found  in  earlier  analyses  at  these  levels.  Cl's  updraft 
remained  split  with  the  northern  branch  dominant  with  vertical  velocities  exceeding 
12.0  m  s'^  at  7.3  km  MSL.  The  southern  branch  contained  vertical  velocities  as 
strong  as  9.0  m  s  '.  Cl  again  developed  an  identifiable  downdraft  of  about  3.0  m  s  ' 
between  3.3  and  4.3  km  MSL.  C2's  updraft  became  very  weak  (approximately  3.0 
m  s  ')  and  narrow  (approximately  0.5  km  wide).  C2's  downdraft  remained  at  about 
3.0  m  s  '  at  3.3  km  MSL. 

At  1146  MST  both  cells  began  to  show  some  signs  of  regeneration.  The  reflecti¬ 
vities  from  both  cells  remained  about  the  same  as  at  1143  MST.  The  horizontal  wind 
analysis  shows  continued  strong  low-level  convergence  near  the  location  of  Cl  (0.0, 
-3.5)  and  increased  low-level  convergence  in  the  vicinity  of  C2  (2.25,  -1.25).  Strong 
upper-level  divergence  continued  over  Cl  (Fig.  11).  The  split  updraft  continued  in 
Cl  with  the  northern  branch  remaining  dominant.  A  new  updraft  formed  at  this  time 
on  the  eastern  edge  of  Cl  (over  the  western  slopes  of  R2).  This  updraft  was  approxi¬ 
mately  1.0  km  wide  and  reached  a  maximum  strength  of  9.0  m  s  '  before  merging 
with  the  other  two  branches  of  Cl's  updraft  at  8.3  km.  C2's  updraft  remained  very 
narrow,  but  increased  in  strength  to  over  9.0  m  s  '  between  6.3  and  8.3  km  MSL  with 
a  maximum  of  12.0  m  s  '  at  7.3  km  MSL.  Cl's  downdraft  remained  3.0  m  s  ',  but 
increased  in  vertical  extent  to  5.3  km  MSL.  C2's  downdraft  remained  unchanged 
from  the  previous  analysis  time  (Fig.  12). 

The  regeneration  trend  of  Cl  (0.0,  -4.25)  continued  at  1149  MST  as  its  maximum 
reflectivity  increased  to  50  dBZ.  C2's  (3.0,  -1.75)  maximum  reflectivity  remained 
steady  at  40  dBZ,  but  extended  in  height  to  7.3  km  MSL  from  the  previously  ob- 
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Fig.  I  la.  Horizontal  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
as  in  Fig.  ft,  except  at  1146  MST  and  at  3.3  km  MSL. 
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Fig.  lib.  As  in  Fig.  8,  horizontal  cross-section  of  observed  radar  reflectivities  and 
tor  velocities  for  1146  MST  and  at  9.3  km  MSL. 


40 


07/31/814 

1146  MST  3.3  KM  MSL 

4  RADARS 


W  COMPONENT  CM/S) 


Fig.  1 2a.  Horizontal  cross-section  of  observed  vertical  velocities  at  an  altitude  of  3.3 
km  relative  to  mean  sea  level  at  1146  MST  31  July  19S4.  Velocities  are  contoured  ev¬ 
ery  3.0  m  s  '  .starting  at  0  m  s  '.  Solid  contours  denote  upward  vertical  velocities 
(positive  w/.  and  dashed  contours  denote  downward  vertical  velocities  (negative  w). 
Origin  and  axes  are  as  in  Fig.  H. 
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Fig.  12h.  Horizontal  cross-section  of  observed  vertical  velocities  at  1146  MST  31  July 
19H4  as  in  Fig.  P'  'cccept  at  8.3  km  MSL. 
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served  3.3  km  MSL.  Low-level  convergence  continued  in  the  vicinity  of  Cl  and  C2, 
while  strong  divei^ence  continued  over  Cl.  The  reflectivity  analysis  also  shows  the 
development  of  a  new  cell  (-0.5,  -7.S)  (denoted  hereafter  as  C3)  as  an  elongation  in 
the  reflectivity  contours  along  Rl's  ridge  axis  (Fig.  13).  Cl's  updraft  remained  split 
at  this  time.  The  new  updraft  noted  on  the  eastern  edge  of  cell  Cl  at  1146  MST 
dominated  Cl  with  a  maximum  vertical  velocity  of  18.0  m  s'‘  at  6.3  km  MSL  (1.25, 
-4.0).  Cl's  southern  updraft  remained  relatively  strong  with  a  maximum  updraft  of 
12.0  m  s'*  (0.0,  -5.5).  The  two  updrafts  merged  with  each  other  and  with  the  rem¬ 
nants  of  the  older  third  updraft  at  8.3  km  MSL.  C2's  updraft  (3.25,  -2.5)  remained 
relatively  weak,  barely  reaching  9.0  m  s'*  at  8.3  km  MSL.  C3  contained  no  distin¬ 
guishable  updrafts  or  downdrafts  at  this  time.  The  downdrafts  fit>m  both  Cl  and  C2 
combined  to  form  a  broad  area  of  downward  motion  covering  approximately  12  km" 
below  4.3  km  MSL  northwest  of  the  updrafts.  The  downdrafts  reached  a  maximum 
strength  of  6.0  m  s'*  (Fig.  14). 

By  1 152  MST,  Cl  (-0.25,  -4.0)  had  drifted  slowly  northward.  C2  lost  its  identity 
after  slowly  weakening  and  being  absorbed  by  the  northward-moving  Cl.  Cl  re¬ 
mained  strong  with  its  maximum  reflectivity  exceeding  50  dBZ.  C3  (-0.25,  -6.5) 
showed  further  signs  of  development  with  its  maximum  reflectivity  exceeding  30  dBZ 
between  3.3  and  8.3  km  MSL.  The  horizcmtal  wind  analysis  continued  to  show 
strong  convergence  over  both  R1  and  R2,  especially  at  4.3  km  MSL.  Strong  diver¬ 
gence  existed  above  both  Cl  and  C3  at  and  above  8.3  km  MSL.  Cl  continued  with  a 
very  complex  updraft  structure  below  8.3  km  MSL.  Several  small  updrafts  formed  at 
this  time  below  8.3  km  MSL,  especially  along  the  southeastern  edge  of  the  cell. 
These  were  likely  caused  by  the  convergence  of  the  downdraft  air  as  it  spilled  down 
the  valleys  and  the  southeasterly  flow  up  the  valleys  ahead  of  the  storm.  The  updrafts 


43 


07/31/84 

1149  MST  7.3  KM  MSL 

4  RflOflRS  — ►  15.0M/S 


HORIZONTAL  FLOW  FIELD/REFLECTIVITT  (DBZ) 


Fig.  13.  As  in  Fig.  8.  horizontal  cross-section  of  observed  radar  reflectivities  and  vec¬ 
tor  velocities  for  1149  MST  and  at  7.3  km  MSL. 
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Fig.  1 4a.  Horizontal  cross-section  of  observed  vertical  velocities  as  in  Fig.  1 2a,  except 
at  1149  MST  and  at  3.3  km  MSL  Labels  indicate  locations  of  updrafts  associated  with 
cells  Cl  and  C2. 
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Fig.  I4h.  A.i  in  Fig.  1 2a.  horizontal  cross-section  of  observed  vertical  velocities  for 
1149  MST  and  at  6.3  km  MSL.  Labels  indicate  locations  of  updrafts  associated  with 
cells  Cl  and  C2. 
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combined  to  form  the  main  Cl  updraft  with  a  maximum  strength  of  12.0  m  s'.  C3's 
narrow  updraft  (approximately  0.5  km  wide)  strengthened  markedly  from  the  pre¬ 
vious  analysis  to  a  maximum  strength  of  15.0  m  s'.  Cl's  downdraft  continued  to 
cover  a  broad  area  with  a  maximum  strength  of  6.0  m  s  '.  The  vertical  cross-sections 
shown  in  Figure  15  show  that  Cl  had  a  significant  backward  tilt.  As  a  result,  the  pre¬ 
cipitation  and  the  precipitation-induced  downdraft  occurred  north  and  west  of  the  up¬ 
draft  so  that  they  did  not  interfere  with  it.  Similar  storm  structures  and  circulations 
have  been  observed  in  some  long-lived  severe  storms  (see  e.g.,  Browning  et  al., 
1976). 

Poor  data  quality  at  1155  MST  prevented  an  analysis  for  that  time.  At  1158 
MST,  Cl  (0.5,  -3.5)  reached  its  peak  intensity,  and  C3  (-0.5,  -6.5)  increased  in 
strength  slightly.  Cl's  area  of  50  dBZ  reflectivities  grew  substantially  from  1152 
MST,  and  C3's  reflectivities  increased  to  40  dBZ.  A  new  cell  (5.0,  -3.5)  (denoted 
hereafter  as  C4)  developed  at  this  time  over  R3  (Fig.  16).  The  maximum  reflectivity 
from  C4  reached  40  dBZ  at  7.3  km  MSL.  Strong  low-level  convergence  continued 
over  both  R1  and  R2.  The  horizontal  wind  analysis  showed  strong  divergence  over 
Cl  between  10.3  and  12.3  km  MSL  and  over  C3  and  C4  at  8.3  km  MSL.  Cl's  mul¬ 
tiple  updrafts  continued,  but  weakened  considerably  at  low  levels.  The  analysis 
showed  updrafts  greater  than  9.0  m  s'*  only  above  8.3  km  MSL  with  maximum  up¬ 
drafts  of  12.0  m  s'.  The  updrafts  merged  at  11.3  km  MSL.  The  downdraft  asso¬ 
ciated  with  Cl  became  very  broad  and  strong  by  this  time,  exceeding  12.0  m  s  '  at 
3.3  km  MSL.  The  downdraft  impacted  the  surface  very  near  the  top  of  VI.  The 
analysis  also  showed  a  strong  updraft  (maximum  strength  18.0  m  s  ')  associated  with 
C3  located  directly  over  Rl.  The  vertical  cross-sections  shown  in  Figure  17  show  the 
strong  interrelationships  between  the  circulations  of  all  three  cells.  For  example,  the 
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Fig.  1 5a.  Horizontal  cross-section  of  observed  radar  r^lectivities  as  in  Fig.  lOa.  ex¬ 
cept  at  1152  MST  and  4.3  km  MSL.  Lines  AB  and  CD  show  locations  of  vertical  cross- 
sections  in  Figs.  15b  and  15c.  respectively. 
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Fig.  15b.  A.s  in  Fig.  /  Oh,  vertical  cross-section  of  observed  radar  reflectivities  and  vec¬ 
tor  velocities  at  1152  MST.  Cross-section  location  AB  is  shown  in  Fig.  15a. 
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Fig.  15c.  As  in  Fig.  lOb,  vertical  cross-section  of  observed  radar  reflectivities  and  vec¬ 
tor  velocities  at  1152  MST.  Cross-section  location  CD  Ls  shown  in  Fig.  15a. 


07/31/84 

1158  MST  8.3  KM  MSL 

4  RADARS  —►  15.0M/S 


HORIZONTAL  FLOW  FIELD/REFLECTIVITT  [DBZl 


Fig.  /  6.  Horizontal  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
as  in  Fig.  8.  except  at  1158  MST  and  at  8.3  km  MSL. 
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Fig.  1 7a.  As  in  Fig.  lOa,  horizontal  cross-section  of  observed  radar  reflectivities  at 
/ 1 58  MST  and  4.3  km  MSL.  Lines  AB  and  CD  show  locations  of  vertical  cro.ss-sections 
in  Figs.  I7h  and  1 7c.  respectively. 
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Fig.  1 7h.  Vertical  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
in  Fig.  I  Oh.  except  at  1158  MST.  Cross-section  location  AB  is  shown  in  Fig.  1 7a. 
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Fig.  1 7c.  Vertical  cross-section  of  observed  radar  reflectivities  and  vector  velocities  as 
in  Fig.  I  Oh.  except  at  1 158  MST.  Cross-section  location  CD  is  shown  in  Fig.  1 7a. 
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outflow  at  the  top  of  C3,  the  newer  cell,  feeds  directly  into  the  much  higher  updraft 
in  the  older  cell,  Cl.  Additionally,  the  flow  through  C4  feeds  directly  into  Cl's  up¬ 
draft. 

3.3.4  Storm  Decay 

The  1201  MST  analysis  suffered  from  a  lack  of  data  over  most  of  the  southern 
half  of  the  analysis  grid.  The  available  data  showed  Cl  (O.S,  -3.0)  began  a  slight 
weakening  trend  at  this  time.  The  maximum  reflectivities  observed  from  Cl  de¬ 
creased  to  40  dBZ.  C3  (-0.5,  -6.5)  moved  slightly  northward  and  weakened  slightly, 
but  radar  still  detected  45  dBZ  returns.  The  actual  amount  of  weakening  is  question¬ 
able  due  to  the  lack  of  data.  C4  (5.0,  -3.0)  continued  to  strengthen  with  40  dBZ  re¬ 
turns  detected  between  6.0  and  8.3  km  MSL.  The  horizontal  wind  analysis  showed 
strong  mid-level  divergence  southeast  of  the  storm  area  over  the  eastern  slope  of  Rl. 
Additionally,  the  analysis  showed  strong  upper-level  divergence  over  both  Cl  and  C3. 
The  updrafts  in  both  Cl  and  C3  weatened  at  this  time:  Cl  to  a  maximum  strength  of 
9,0  m  s'‘,  C3  to  a  maximum  of  12.0  m  s'‘.  C4's  updraft  intensified  rapidly  to  12.0 
ms'.  Cl's  downdraft  remained  broad  and  strong  with  a  maximum  strength  of  12.0 
m  s  '  at  3.3  km  MSL.  The  missing  data  prevented  analysis  of  C3's  downdrafts.  C4's 
rapid  intensification  led  to  the  development  of  a  narrow,  but  strong  downdraft  (9.0 
m  s  '). 

The  reflectivity  analysis  at  1204  MST  showed  only  minor  changes  in  Cl  (0.5, 
-2.5)  and  C3  (0.0,  -5.5)  from  the  1201  analysis.  By  1204,  C4  (3.5,  -2.5)  had  drifted 
northwestward  to  the  location  where  C2  had  previously  dissipated,  and  its  area  of  40 
dBZ  returns  extended  downward  to  3.3  km  MSL.  The  radar  detected  30  dBZ  returns 
between  4.3  and  7.3  km  MSL  from  a  new  cell  (-4.0,  -5.0)  (denoted  hereafter  as  C5) 
just  west  of  Rl  (Fig.  18).  The  horizontal  wind  analysis  showed  a  weakening  of  the 
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Fig.  IS.  As  in  Fig.  8,  horizontal  cross-section  of  observed  radar  reflectivities  and  vec 
tor  velocities  for  1204  MST  and  at  6.3  km  MSL. 
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channelling  effect  in  VI  with  a  slight  increase  in  the  convergent  channelling  in  V2. 
Moderate  low-level  convergence  existed  in  the  region  between  Cl  and  C4  directly 
over  the  southern  tip  of  R2.  The  mid-level  divergence  continued  over  the  eastern 
slopes  of  R1  and  a  new  divergence  area  formed  over  the  western  slt^s  of  R1  in  the 
vicinity  of  C5.  The  strongest  upper-level  divergence  existed  over  C4.  The  vertical 
velocity  analysis  for  1204  MST  showed  a  slight  increase  in  updraft  strength  in  Cl, 
C3,  and  C4.  The  analysis  showed  no  distinct  updraft  regions  associated  with  C5. 
The  downdrafts  in  Cl  and  C4  weakened  to  6.0  m  s  '  and  covered  a  slightly  smaller 
area  than  earlier.  The  analysis  showed  a  well-established  (6.0  m  s  ')  downdraft  asso¬ 
ciated  with  C3.  The  vertical  cross-section  (Fig.  19)  shows  the  strong  development  of 
C4  and  that  its  circulation  continued  to  feed  into  Cl’s  circulation. 

The  reflectivity  analysis  at  1207  MST  showed  very  few  changes  from  the  1204 
and  1201  MST  analyses.  By  this  time  both  Cl  (0.5,  -2.5)  and  C3  (0.0,  -5.5)  had 
drifted  slowly  northward.  Radar  detected  several  new  echoes  developing  over  the 
eastern  slopes  of  R3,  These  appear  to  have  developed  as  a  result  of  the  strong  ups- 
lope  flow  into  C4  (3.5,  -2.5).  The  horizontal  wind  analysis  showed  continued  strong 
low-  and  mid-level  convergence  in  the  vicinity  of  Cl.  The  analysis  showed  strong 
upper-level  divergence  over  both  C3  and  C4.  Cl’s  updraft  remained  at  about  the 
same  strength  as  at  1201,  while  C3's  weakened  slightly.  C4’s  updraft  weakened  con¬ 
siderably  to  a  strength  of  6.0  m  s'*.  C5’s  (-3.75,  -4.75)  updraft  became  more  orga¬ 
nized,  reaching  a  strength  of  6.0  m  s'.  All  four  active  cells  showed  organized 
downdrafts  at  this  time  with  Cl's  and  C3's  increasing  in  strength  again  to  over  9.0 
ms'.  The  vertical  velocity  analyses  showed  that  by  this  time  the  updrafts  existed 
mainly  at  mid-  and  upper-levels. 


Fig.  1 9a.  A.s  in  Fig.  lOa.  horizontal  cross-section  of  observed  radar  reflectivities  at 
1204  MST  and  4.3  km  MSL.  Line  CD  shows  location  of  vertical  cro!:s-section  in  Fig. 
1 9h.  Line  AB  was  not  used. 
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Fig.  /  9h.  V ertical  cross-section  of  observed  radar  reflectivities  and  vector  velocities 
in  Fig.  /  Oh,  except  at  1204  MST.  Cross-section  location  CD  is  shown  in  Fig.  1 9a. 
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By  1210  MST,  the  defining  reflectivity  contours  of  Cl  (0.5,  -2.0)  and  C3  (0.5, 
-4.5)  had  merged.  Radar  returns  from  C4  (3.0,  -2.5)  weakened  slightly  while  those 
from  C5  (-3.5,  -4.0)  increased  in  strength  to  40  dBZ  at  3.3  km  MSL.  The  weak  de¬ 
velopment  over  R3  continued  while  several  new  cells  formed  in  the  extreme  north¬ 
eastern  comer  of  the  analysis  area.  The  horizontal  wind  analysis  showed  strong 
low-level  convergence  over  the  western  slopes  of  R2  and  over  the  northern  end  of  V2. 
Strong  upper-level  divergence  also  existed  over  these  same  regions.  Cl's  updraft  re¬ 
mained  as  its  onl  distinguishing  feature  at  1210  MST.  At  this  time,  it  still  main¬ 
tained  a  strength  of  9.0  m  s  ‘.  C3's  updraft  remained  about  the  same  strength  as  at 
1207  MST,  but  C4's  updraft  strengthened  again  to  9.0  m  s',  and  C5's  updraft 
reached  9.0  m  s  '  for  the  first  time.  At  this  analysis  time.  Cl's  downdraft  had  de¬ 
creased  significantly  to  less  than  3.0  m  s'.  C3,  C4,  and  C5  all  contained  maximum 
downdrafts  of  6.0  m  s  '. 

By  1213  MST,  C3  (0.5,  -3.5)  had  overtaken  all  remnants  of  Cl  and  most  of  the 
remnants  of  C4  (2.5,  -2.0).  C3's  maximum  reflectivity  increased  to  over  50  dBZ  at 
3.3  km  MSL  at  approximately  the  midpoint  between  the  previous  locations  of  C3, 
Cl,  and  C4.  C4  appeared  in  the  reflectivity  analysis  as  only  a  bulge  in  the  40  dBZ 
contour  at  3.3  km  MSL.  C5's  (-3.0,  -4.0)  maximum  reflectivity  remained  about  40 
dBZ,  however  the  maximum  expanded  in  vertical  extent  to  5.3  km  MSL.  A  new  cell 
(2.0,  -7.0)  (denoted  hereafter  as  C6)  developed  along  the  southeastern  edge  of  C3  al¬ 
most  directly  over  the  center  of  VI.  This  cell  was  best  defined  in  the  reflectivity 
analyses  at  6.3  and  7.3  km  MSL  (Fig.  20).  The  horizontal  wind  analysis  showed 
weak  low-level  flow  over  the  northern  half  of  the  analysis  domain,  but  a  relatively 
strong  southwesterly  flow  over  the  southern  half.  The  analysis  showed  very  little 
low-level  convergence.  Some  upper-level  divergence  still  existed  over  C3  and  C6. 


60 


07/31/84 

1213  MST  7.3  KM  MSL 

4  RADARS  — ►  15.0M/S 


HORIZONTAL  FLOW  FIELD/REFLECTIVITT  CDBZ) 


Fig.  20.  Horizontal  cro.ss-section  of  observed  radar  reflectivities  and  vector  velocities 
as  in  Fig.  H,  except  at  1213  MST  and  at  7.3  km  MSL. 
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C3’s  updraft  continued  with  a  maximum  strength  of  9.0  m  s  '  between  5.3  and  9.3  km 
MSL  while  C4's  updraft  became  indistinguishable  at  this  time.  C5's  updraft  contin¬ 
ued  to  strengthen,  reaching  a  maximum  strength  of  12.0  m  s'*.  By  1213  MST,  C6’s 
updraft  had  already  reached  9.0  m  s  '  between  5.3  and  6.3  km  MSL.  The  area  cov¬ 
ered  by  C3's  downdraft  continued  to  shrink,  although  the  downdraft's  maximum 
strength  continued  at  6.0  m  s  '.  The  weak  downdraft  (3.0  m  s  ')  remained  as  the  only 
structure  that  distinguished  C4  from  C3  at  1213  MST.  C5's  downdraft  increased  in 
strength  to  6.0  m  s''. 

At  1216  MST,  C3  (0.5,  -4.0)  had  absorbed  the  final  remnants  of  C4.  C3  contin¬ 
ued  with  a  maximum  reflectivity  of  40  dBZ.  C5  (-3.0,  -4.0)  maintained  its  previous 
strength,  while  C6  (2.0,  -7.0)  strengthened  rapidly,  with  its  40  dBZ  maximum  reflec¬ 
tivities  extending  down  to  the  lowest  analysis  level  of  3.3  km  MSL.  A  weak  new  cell 
(5.5,  -2.0)  (denoted  hereafter  as  C7)  developed  at  this  time  over  R3.  The  horizontal 
wind  analysis  at  this  time  again  showed  very  weak  flow  over  the  northern  half  of  the 
analysis  domain  with  a  relatively  strong  southerly  to  southwesterly  flow  over  the 
southern  half.  Some  weak  low-level  convergence  existed  just  east  of  C6.  The  verti¬ 
cal  velocity  analysis  showed  that  by  this  time  most  of  the  vertical  motion  had  weak¬ 
ened  considerably.  Only  C6's  updraft  exceeded  9.0  m  s  '.  Downdrafts  of  6.0  m  s  ' 
existed  only  in  C3  and  C5  by  this  time. 

3.3.5  Summary 

The  previous  subsections  provide  a  detailed  description  of  the  initiation  and  evolu¬ 
tion  of  the  31  July  1984  storms  as  observed  by  the  Doppler  radars.  These  observa¬ 
tions  provide  some  insight  into  the  general  nature  of  the  storms. 

Early  cell  development  occurred  directly  over  or  slightly  east  of  the  ridge  lines. 
These  observations  are  consistent  with  the  theory  that  the  New  Mexico  storms  evolve 
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mainly  as  a  result  of  thermally  driven  mountain-valley  circulations  and  are  focused  by 
local  topographic  features  (see,  for  example,  Barker  and  Banta,  1985;  Banta  and 
Schaaf,  1987;  Schaaf  era/.,  1988).  The  weak  low-level  southeasterly  flow  into  the 
convergent  valleys,  VI  and  V2,  produced  areas  of  upslope  and  convergent  flow  at  the 
northwestern  (upslope)  end  of  those  valleys.  The  upslope  flow  and  convergence  fo¬ 
cused  the  areas  of  convective  initiation  over  the  ridges  nearest  the  ends  of  VI  and  V2 
as  described  by  Banta  and  Schaaf  (1987). 

After  the  first  cells  had  developed,  the  mountain-valley  circulation  and  the  conver¬ 
gent  flows  continued  to  help  focus  the  location  of  new  cell  initiation.  However,  cir¬ 
culations  from  the  existing  cells  also  began  to  affect  the  location  of  cell  initiation. 
For  example,  at  1152  MST,  several  short-lived  cells  developed  along  the  outflow 
boundaries  from  Cl  and  C3.  In  some  cases  the  low-level  flow  induced  by  the  exist¬ 
ing  cells  enhanced  or  created  upslope  and/or  convergent  flows  which  in  turn  helped 
focus  new  cell  development. 

The  storm  observed  near  Langmuir  Laboratory  on  31  July  1984  generally  had  a 
lifetime  of  approximately  1.0  to  1.5  hr.  The  storm  contained  several  cells  with  life¬ 
times  of  from  20  to  60  min.  Many  of  the  cells  exhibited  periodic  intensification  and 
weakening  with  a  period  of  6  to  12  min.  For  example,  both  Cl  and  C2  underwent  at 
least  two  cycles  of  intensification/ weakening  during  the  analysis  period  (Cl  from 
1134  to  1143  MST  and  1146  to  i210  MST;  C2  from  before  1134  to  1143  MST  and 
1146  to  1152  MST).  This  periodic  activity  resembles  closely  the  updraft  pulsations 
suggested  by  Ziegler  et  al.  (1986).  Radar  returns  from  the  strongest  cells  exceeded  50 
dBZ  several  times,  but  maximum  returns  from  most  cells  ranged  from  30  to  40  dBZ. 
Most  of  the  cells  never  extended  higher  than  10.0  km  MSL,  however,  the  stronger 
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cells.  Cl  and  C4,  frequently  extended  above  10.0  km  and  in  some  cases  to  over  12.0 
km  MSL. 

Most  cells  contained  single  updrafts  approximately  2.0  km  wide  with  maximum 
vertical  velocities  of  12.0  m  s'*  or  less.  Cl  frequently  contained  multiple  low-level 
updrafts  which  merged  into  a  single  updraft  at  mid-levels.  The  strongest  observed 
updrafts  of  18.0  m  s'*  occurred  in  Cl  and  C3.  In  the  early  stages  of  cell  develop¬ 
ment,  the  downdrafts  were  only  about  1.5  km  wide  with  maximum  vertical  velocities 
of  3.0  m  s'*  or  less.  As  the  cells  matured,  however,  the  downdrafts  expanded  in  areal 
coverage  and  increased  in  strength  with  maximum  vertical  velocities  frequently  ex¬ 
ceeding  6.0  m  s  *  and  with  an  absolute  maximum  of  12.0  m  s'*.  In  at  least  one  cell, 
Cl,  during  the  cell's  mature  stage,  the  interaction  between  the  updraft  and  the  down- 
draft  appeared  to  be  symbiotic  in  a  manner  similar  to  that  often  observed  in  long-lived 
severe  storms.  As  discussed  in  a  previous  section,  the  updraft  tilted  northwestward 
(or  backward  relative  to  the  low-level  storm  inflow).  The  precipitation  and  its  asso¬ 
ciated  downdraft  were  located  behind  and  below  the  updraft  (relative  to  the  low-level 
inflow)  and  thus  did  not  interfere  with  the  inflow  or  the  updraft.  Several  other  cells 
contained  tilted  updrafts,  but  to  a  lesser  extent  than  Cl.  As  the  cells  entered  the 
decay  stage,  the  updrafts  became  nearly  vertical  so  that  the  precipitation  and  the 
downdraft  interfered  considerably  with  the  updraft. 

The  Doppler  radar  observations  and  the  associated  wind  analyses  reveal  a  strong 
interrelationship  between  all  the  cells  in  the  analysis  domain.  In  the  mature  stage,  the 
outflow  from  downdrafts  appeared  in  some  cases  to  enhance  low-level  convergence 
that  either  focused  the  development  of  new  cells  or  helped  sustain  existing  cells.  In 
several  instances  the  updrafts  of  two  or  more  cells  merged  into  one  updraft  at  upper- 
levels.  Often  the  updraft  of  a  newer  cell  developed  in  the  low-  to  mid-levels  and  fed 
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into  a  higher-level  updraft  associated  with  an  older  cell.  Additionally,  the  divergent 
flow  located  near  the  top  of  a  strong  cell  often  cut  off  the  updraft  of  a  weaker  neigh¬ 
boring  cell.  The  periodic  intensification  noted  above  frequently  appeared  to  alternate 
among  the  cells.  For  example,  during  Cl's  first  period  of  intensification,  C2  was 
weakening.  This  certainly  did  not  occur  without  exception,  but  the  trend  was  readily 
apparent.  Finally,  the  observations  showed  several  cases  where  older  or  weaker  cells 
were  overtaken  and  absorbed  by  newer  and/or  stronger  cells. 


Chapter  4 

The  Numerical  Model 

This  chapter  describes  the  numerical  model  chosen  for  this  study.  Section  4.1  de¬ 
scribes  the  criteria  used  to  select  the  numerical  model.  Section  4.2  discusses  the  mod¬ 
el's  ability  to  simulate  mountain  thimderstorms.  In  section  4.3  the  model  selected  for 
this  study  is  briefly  described.  Sections  4.4  and  4.5  describe  the  model's  vertical  and 
horizontal  grid  structures,  respectively,  and  section  4.6  describes  the  topography  used 
in  the  numerical  simulations.  Finally,  the  model  initialization  procedure  is  discussed 

in  section  4.7. 

4.1  Model  Selection 

Before  conducting  any  numerical  study,  a  numerical  model  must  be  selected. 
Choosing  a  model  is  extremely  important  and  should  be  deliberate  and  well-planned. 
Essential  criteria  must  be  defined  and  used  in  the  selection  process.  Flexibility,  appli¬ 
cability,  and  availability  must  also  be  considered.  If  existing  models  do  not  meet  the 
needs  of  the  study,  then  a  new  model  must  be  developed. 

To  study  the  type  of  mountain  thunderstorms  of  interest  here,  it  was  necessary  to 
choose  a  model  with  fully  compressible,  non-hydrostatic  equations  of  motion;  the  op¬ 
tion  of  conducting  two-  or  three-dimensional  simulations;  a  terrain-following  coordi¬ 
nate  system;  a  surface  layer  parameterization  and  a  soil  model  which  interact  with 
realistic  long-  and  short-wave  radiation  parameterizations;  and  a  cloud  microphysical 
parameterization  which  includes  the  ice  phase.  Additionally,  the  model  needed  to  be 
compatible  with  locally  available  computer  software  and  hardware.  The  Colorado 
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State  University  (CSU)  Regional  Atmospheric  Modeling  System  (RAMS)  was  chosen 
for  this  study.  RAMS  was  chosen  over  other  potential  models  because  it  met  all  of 
the  above  criteria,  because  of  its  flexibility,  its  applicability  to  modeling  mountain 
circulations  and  storms,  and  its  availability. 

4.2  Model  Ability  To  Simulate  Mountain  Storms 

Before  undertaking  any  numerical  study,  it  is  essential  to  assess  the  model's  ability 
to  simulate  the  basic  processes  responsible  for  the  phenomenon  under  study.  This 
study  is  concerned  with  the  initiation  and  development  of  mountain  thunderstorms  un¬ 
der  conditions  of  strong  solar  heating  and  weak  ambient  flow.  As  discussed  in  Chap¬ 
ter  3,  the  3  processes  are  expected  to  include  primarily  thermodynamic  processes  and 
possibly  dynamic  processes  such  as  orographic  lifting  and  valley  channeling.  Only 
after  verifying  the  model's  skill  with  these  processes  can  we  use  it  with  confidence  to 
study  the  more  complex  phenomenon  of  mountain  thunderstorms. 

The  simple,  three-dimensional,  no-wind  simulation  used  by  Klemp  and  Wilhelm- 
son  (1978a)  to  test  their  model  also  provided  an  excellent  test  for  this  model's  skill 
with  thermodynamic  processes.  In  this  test,  simulated  convection  is  initiated  by  in¬ 
serting  a  warm,  moist,  symmetric  perturbation  in  the  lower  part  of  an  otherwise  ho¬ 
mogeneous  atmosphere.  Because  this  type  of  perturbation  should  yield  an 
axisymmetric  storm,  this  type  of  test  can  immediately  reveal  many  coding  errors. 
This  was  an  important  benefit  in  this  case  because  RAMS  was  revised  and  recoded 
shortly  before  the  start  of  this  study.  The  results  of  this  test  did,  in  fact,  reveal  sever¬ 
al  coding  errors.  After  these  errors  were  corrected,  the  test  results  agreed  very  well 
with  those  reported  by  Klemp  and  Wilhelmson.  Some  minor  differences  were  noted 
between  the  two  tests.  These  could  be  attributed  to  the  different  physical  parameter- 
izations  used  by  the  two  models.  Tripoli  (1986)  tested  the  model's  ability  to  simulate 
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simple  gravity  waves  and  reported  good  agreement  between  his  simulations  and  those 
of  Durran  (1981). 

In  addition  to  these  simple  tests,  the  model  has  been  used  for  simulations  of  more 
complex  processes  important  to  mountain  convection.  Banta  (1983,  1986),  Bader  and 
McKee  (1983,  1985),  and  Bader  etal.  (1987)  used  the  cloud  model  portion  of  RAMS 
to  study  mountain  upslope  flow  and  boundary  layer  evolution.  Cotton  and  Tripoli 
(1978),  Tripoli  and  Cotton  (1980,  1982,  1986),  and  Cotton  etal.  (1982)  also  used  the 
cloud  model  portion  of  RAMS  successfully  in  studies  of  cumulus  and  thunderstorm 
development.  The  simple  Klemp-Wilhelmson  tests  described  above,  along  with  the 
successful  complex  simulations  carried  out  by  other  investigators,  indicate  that  the 
model  is  well-suited  to  simulate  the  initiation  and  development  of  mountain  thunder¬ 
storms  of  interest  in  this  study. 

4.3  Model  Description 

As  its  name  implies,  the  CSU  RAMS  is  a  modeling  system,  not  a  single  model. 
RAMS  is  a  combination  of  the  Tripoli  and  Cotton  (1982)  non-hydrostatic  cloud  model 
and  of  two  hydrostatic  mesoscale  models  (Tremback  etal.,  1985;  Mahrer  and  Pielke, 
1977).  Because  of  this  combination,  RAMS  is  very  flexible  and  allows  the  user  to 
select  from  a  wide  variety  of  physical  processes  and  numerical  techniques. 

This  study  employed  the  non-hydiostatic  cloud  model  portion  of  RAMS.  This 
section  briefly  describes  the  characteristics  of  the  non-hydrostatic  model  and  the  vari¬ 
ous  options  employed  in  this  study.  The  description  here  closely  follows  that  given  by 
Tripoli  (1986). 

The  model  formulation  is  a  revised  and  slightly  modified  form  of  the  CSU 
cloud/mesoscale  model  described  by  Tripoli  and  Cotton  (1982)  and  Cotton  et  al. 
(1982,  1986).  The  model  employs  a  non-hydrostatic,  quasi-Boussinesq,  fiilly  com- 


68 


pressible  set  of  primitive  equations  cast  in  two  or  three  dimensions  (Tripoli  and  Cot¬ 
ton,  1986;  Tripoli,  1986).  The  current  model  formulation  follows  Klemp  and 
Wilhelmson  (1978a)  and  predicts  the  Exner  function  (n)  and  diagnoses  density 
(Tripoli,  1986;  Bader  eta/.,  1987). 

The  model  uses  time-dependent  partial  differential  equations  to  predict  the  three 
wind  components  (u,  v,  w),  the  Exner  function  (n)  ,  the  liquid-ice  potential  tempera¬ 
ture  (6,/)  (Tripoli  and  Cotton,  1981),  total  water  mixing  ratio  (r, ),  mixing  ratios  of 
rain  water  (r, ),  pristine  ice  crystals  (r, ),  graupel  or  hail  (r^ ),  and  aggregated  snow 
flakes  (r^ ).  The  model  also  predicts  the  total  concentration  of  pristine  ice  crystals 
(A^^^,).  This  gives  a  total  of  eleven  prognostic  variables.  The  model  uses  a  diagnostic 
procedure  described  by  Tripoli  and  Cotton  (1982)  to  diagnose  potential  temperature 
(6),  temperature  (f ),  cloud  water  mixing  ratio  (r^ ),  water  vapor  mixing  ratio  (/\, ), 
pressure  (p),  and  density  (p). 

The  model  includes  parameterizations  for  several  important  physical  processes. 
Microphysical  processes  are  treated  by  a  generalized  microphysical  parameterization 
scheme  described  by  Tripoli  et  al.  (1988).  Within  this  generalized  framework  it  is 
possible  to  vary  parameters  to  formulate  the  desired  microphysical  parameterization. 
The  microphysical  parameterizations  chosen  for  this  study  are  similar  to  those  of 
Tripoli  and  Cotton  (1980)  for  warm  rain  processes  and  Cotton  etal.  (1982,  1986)  for 
ice  phase  processes.  Long-  and  short-wave  radiative  tendencies  on  0^  are  computed 
following  Chen  and  Cotton  (1983).  For  the  simulations  conducted  in  this  study,  the 
radiative  tendencies  were  updated  at  120  s  intervals  starting  at  /  =  900  s.  Surface 
layer  fluxes  of  heat,  moisture,  and  momentum  are  parameterized  following  the  meth¬ 
ods  of  Louis  (1979)  and  Manton  and  Cotton  (1977).  A  surface  roughness  of  1.0  cm 
was  prescribed  for  the  surface  layer  parameterization.  Soil  temperature  and  soil  mois- 
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ture  are  predicted  by  an  underlying  soil  model  described  by  McC umber  and  Pielke 
(1981)  and  modified  by  Tremback  and  Kessler  (198S).  The  soil  model  inclutted  five 
levels  located  at  the  surface  and  1,4,  16,  and  32  cm  beneath  the  surface.  Chapter  5 
provides  further  details  on  the  soil  parameters  prescribed  for  the  simulations.  Turbu¬ 
lent  mixing  above  the  surface  layer  is  modeled  using  the  eddy  viscosity  approach  de¬ 
scribed  by  Lilly  (1962).  Mixing  is  enhanced  in  unstable  regions  as  suggested  by  Hill 
(1974)  and  implemented  by  Tripoli  (1986). 

The  predictive  equations  are  solved  on  a  staggered  grid  described  by  Tripoli  and 
Cotton  (1982)  with  a  terrain-following  vertical  coordinate  system  similar  to  that  of 
Gal-Chen  and  Sommerville  (1975a,  1975b).  The  time-differencing  scheme  used  here 
was  described  by  Tripoli  and  Cotton  (1982)  and  modified  by  Tripoli  (1986).  The 
time-split  leapfrog  scheme  developed  by  Klemp  and  Wilhelmson  (1978a)  provides  the 
basic  time-differencing  framework.  In  this  scheme  the  terms  responsible  for  sound 
wave  propagation  are  integrated  on  a  small  time  step  while  the  non-acoustic  terms  are 
held  fixed  for  several  small  time  steps  and  updated  on  the  larger  leapfrog  time  step. 
For  the  small  time  step,  the  model  employs  a  forward-backward  differencing  scheme 
for  the  horizontal  terms  and  an  implicit  Crank-Nicholson  scheme  described  by  Durran 
(1981)  and  Tripoli  and  Cotton  (1986)  for  the  vertical  terms.  The  physical  p)arameter- 
izations  are  all  computed  with  a  forward  time-differencing  scheme.  The  turbulent 
mixing  parameterization  uses  a  combined  forward  and  Crank-Nicholson  scheme  sug¬ 
gested  by  Paegle  et  al.  (1976).  For  this  study,  the  acoustic  terms  were  integrated  us¬ 
ing  a  2.5  s  small  timestep,  while  the  non-acoustic  terms  were  integrated  with  a  5.0  s 
large  timestep.  All  advective  terms  are  computed  with  a  fourth-order  finite- 
differencing  method. 
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Because  RAMS  is  a  limited  area  model,  ^propriate  boundary  conditions  must  be 
specified  in  order  to  determine  unique  solutions  to  the  time  dependent  equations.  For 
the  top  boundary  condition  a  rigid  upper  lid  similar  to  that  employed  by  Klemp  and 
Wilhelmson  (1978a)  was  used.  The  bottom  boundary  condition  assumed  zero  mean 
vertical  fluxes  at  the  surface  with  surface  layer  turbulent  fluxes  of  momentum,  heat, 
and  water  vapor  diagnosed  by  the  surface  layer  parameterization.  The  open  lateral 
boundary  conditions  described  by  Klemp  and  Wilhelmson  (1978a)  were  employed 
with  a  specified  intrinsic  gravity  wave  phase  velocity  of  30  m  s  '. 

4.4  Horizontal  Grid  Structure 

The  horizontal  grid  structure  for  this  study  was  chosen  to  maximize  the  horizontal 
resolution  over  a  model  domain  large  enough  to  accommodate  the  simulation  of 
mountain  convection  while  minimizing  artificial  model  boundary  effects  on  the  simu¬ 
lated  storms.  Computational  limitations  dictated  the  choice  of  a  medium-sized  do¬ 
main  to  allow  for  a  desired  high-resolution  grid  structure. 

As  discussed  in  Chapter  3,  the  storms  of  interest  here  developed  over  the  Magda¬ 
lena  Mountains,  a  range  that  extends  approximately  28  km  in  the  north-south  direction 
and  approximately  20  km  in  the  east-west  direction.  Most  of  the  convection  on  31 
July  1984  developed  and  remained  near  the  center  of  the  range.  A  33  km  model  do¬ 
main  centered  on  the  Langmuir  Laboratory  was  chosen  for  the  simulations  in  this 
study.  This  domain  was  large  enough  to  allow  ample  room  for  storm  development 
near  the  domain  center  and  also  provided  for  "buffer"  regions  around  the  domain 
boundaries  to  minimize  the  boundary  '*ffects  on  the  domain's  interior. 

A  500  m  horizontal  .,;solution  was  chosen  for  this  study.  This  resolution  was  cho¬ 
sen  based  on  three  considerations.  First,  an  accurate  numerical  simulation  of  the 
types  of  storms  of  interest  here  required  a  high-resolution  grid  to  capture  the  cloud- 
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and  sub-cloud-scale  processes  known  to  be  important  in  mountain  convection.  Se¬ 
cond,  the  extremely  rough  terrain  in  the  model  domain  (see  Section  4.6)  necessitated 
a  high-resolution  grid  to  capture  the  important  terrain  features  and  resultant  effects  on 
storm  development.  Finally,  the  desire  for  the  highest  resolution  possible  was  neces¬ 
sarily  tempered  by  the  availability  of  computer  resources.  The  500  m  grid  structure 
provided  a  suitable  compromise  among  the  considerations  above. 

Several  other  studies  have  employed  similar  grid  structures  and  domain  sizes  to 
simulate  thunderstorm  development  in  general  and  mountain  convection  in  particular. 
For  example,  Klemp  and  Wilhelmson  (1978a)  employed  a  24  km  domain  with  1000 
m  grid  spacing  to  simulate  a  generic  severe  storm;  Tripoli  and  Cotton  (1980)  used  a 
35.25  km  domain  with  750  m  grid  spacing  to  simulate  a  Florida  thunderstorm;  and 
Clark  and  Gall  (1982)  and  Smolarkiewicz  and  Clark  (1985)  used  a  51  km  domain 
with  1000  m  grid  spacing  to  simulate  New  Mexico  and  Montana  thunderstorms,  re¬ 
spectively. 

4.5  Vertical  Grid  Structure 

The  vertical  grid  structure  chosen  for  this  study  reflects  an  attempt  to  adequately 
resolve  most  of  the  important  features  of  thunderstorm  initiation  and  development 
with  a  minimum  of  comput  er  resources.  It  was  computationally  prohibitive  to  use  an 
extremely  high-resolution  vertical  grid  structure  required  to  resolve  such  important 
processes  as  turbulent  fluxes,  gust  fronts,  etc.  Instead,  a  vertical  grid  structure  was 
chosen  which  resolves  the  processes  thought  to  be  most  important  and  at  least  allows  a 
gross  parameterization  of  the  unresolvable  processes. 

Banta  (1984).  Braham  and  Draginis  (1960),  Orville  (1965a),  and  others  have 
shown  that  shallow  mountain-valley  circulations  play  a  critical  role  in  he  develop¬ 
ment  of  mountain  storms.  These  flows  typically  are  only  a  few  hundred  meters  deep. 


72 


Properly  resolving  such  circulations  would  require  a  vertical  grid  spacing  of  100-200 
m.  Tripoli  (1986)  showed  that  a  grid  spacing  of  500-1000  m  is  required  to  adequately 
resolve  gravity  waves  under  conditions  similar  to  those  ooserved  over  the  Magdalena 
Mountains.  Resolving  turbulent  effects  would  require  a  grid  spacing  on  the  order  of 
tens  of  meters. 

It  is  reasonable  to  expect  the  actual  and  simulated  thunderstorms  to  occasionally 
penetrate  the  tropopause.  To  adequately  simulate  this  phenomenon,  the  model's  verti¬ 
cal  domain  should  extend  at  least  several  kilometers  into  the  stratosphere.  The  tropo¬ 
pause  on  31  July  1984  was  located  at  approximately  12.4  km.  The  desired  vertical 
resolution  of  100  m  would  require  a  computationally  prohibitive  160  vertical  grid 
points.  In  order  to  facilitate  the  desired  horizontal  resolution  of  500  m,  a  vertical 
resolution  of  500  m  was  chosen.  Above  the  tropopause,  the  grid  is  stretched  to  a 
maximum  spacing  of  2000  m  at  the  top  to  provide  the  desired  extension  into  the 
stratosphere.  The  complete  vertical  grid  structure  (for  the  vertical  velocity  portion  of 
the  staggered  grid)  is  shown  in  Table  1 .  Note  that  the  z  and  Az  values  given  in  Table 
1  represent  the  vertical  grid  structure  over  flat  terrain  with  a  terrain  elevation  of  zero. 
In  areas  of  mountainous  terrain  the  actual  z  and  Az  values  vary  considerably  from 
those  shown  in  Table  1  as  a  result  of  the  vertical  coordinate  transformation  developed 
by  Gal-Chen  and  Sommerville  (1975a,  1975b). 

Several  previous  studies  have  employed  similar  grid  structures.  For  example, 
Tripoli  and  Cotton  (1986)  used  a  500  m  grid  to  model  an  intense,  quasi-steady  left- 
moving  thunderstorm  over  mountainous  terrain;  Smolarkiewicz  and  Clark  (1985)  used 
a  200  m  grid  to  model  Montana  thunderstorm  development;  Tripoli  (1986)  used  a 
variable  grid  spacing  from  250  m  to  750  m  to  model  mountain-induced  mesoscale 
convective  systems;  and  Banta  (1986)  used  a  100  m  grid  spacing  to  model  the  evolu- 
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Table  1:  Model  vertical  grid  structure  for 
vertical  velocity  grid  points. 


HEIGHT 

(Meters) 

AZ 

(Meters) 

0 

500 

500 

1000 

500 

1500 

500 

2000 

500 

2500 

500 

3000 

500 

3500 

500 

4000 

500 

4500 

500 

5000 

500 

5500 

500 

6000 

500 

6500 

500 

7000 

500 

7500 

500 

8000 

500 

8500 

500 

9000 

500 

9500 

500 

10000 

500 

10500 

500 

11000 

500 

11500 

500 

12000 

500 

12500 

500 

13000 

500 

14000 

1000 

NOTE:  The  model  vertical  grid  uses  a  terrain¬ 
following  coordinate  system  similar  to  that  described 
by  Gal-Chen  and  Sommerville  (1975a.  1975b).  The 
grid  .structure  shown  above  is  valid  only  in  areas 
where  the  terrain  elevation  is  zero  and  the  terrain  is 
flat.  Elsewhere,  the  transformed  values  are  less 
than  those  shown  above. 


74 


tion  of  the  boundary  layer  in  mountainous  terrain. 

4.6  Model  Topography 

The  topography  employed  in  the  model  for  this  study  was  derived  from  actual  3 
arc-seicond  terrain  data  acquired  from  the  National  Cartt^raphic  Information  Center. 
The  model  topography  was  constructed  in  the  following  manner.  First,  the  appropri¬ 
ate  data  were  extracted  from  the  original  terrain  data  set  and  interpolated  to  a  33  X  33 
km  grid  (centered  on  Langmuir  Lab)  with  0.5  km  grid  spacing.  The  data  were  then 
smoothed  three  times  using  a  simple  five-point  smoother.  It  was  necessary  to  smooth 
the  terrain  field  in  order  to  minimize  the  occurrence  of  numerical  instabilities  induced 
by  extremely  rough  terrain.  Around  the  boundaries  of  the  grid,  the  topography  was 
flattened  to  minimize  terrain-induced  effects  at  the  model's  boundaries.  Next,  to  sat¬ 
isfy  an  internal  model  requirement  that  the  lowest  terrain  value  be  zero,  the  lowest 
terrain  elevation  value  in  the  smoothed  field  was  subtracted  from  all  terrain  values  in 
the  grid.  The  resulting  terrain  grid  was  used  in  three-dimensional  test  simulations. 
Figure  21  shows  the  "3-D"  model  topography. 

The  model  topography  used  in  the  two-dimensional  simulations  for  this  study  was 
constructed  by  extracting  the  appropriate  slab  of  terrain  values  (y  =  0.0  km)  from  the 
33  X  33  km  gridded  data  discussed  above.  Again  the  lowest  terrain  elevation  value 
was  subtracted  from  this  subset  of  data  to  meet  the  model  requirement.  Figure  22 
shows  the  final  model  topography  used  in  the  two-dimensional  simulations. 

4.7  Initialization 

RAMS  is  initialized  with  base-state  atmospheric  temperature,  moisture,  and  wind 
profiles  and  base-state  soil  moisture  and  temperature  profiles.  The  numerical  simu¬ 
lations  in  this  study  were  all  initialized  with  base-state  atmospheric  profiles  (described 
in  Chapter  5)  representative  of  the  state  of  the  atmosphere  over  the  Magdalena  Moun- 
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MODEL  TERRAIN 


Fig.  21.  "Three-dimensionar  model  terrain  elevation  data.  Terrain  elevations  are 
relative  to  lowest  terrain  value  in  grid  domain  after  smoothing  as  described  in  text. 
Terrain  elevations  are  labeled  in  dekameters  and  contoured  every  20  dekameters.  The 
horizontal  and  vertical  axes  are  labeled  in  kilometers.  The  origin  (0.0)  is  located  at 
Langmuir  Laboratory. 


z 


TWO-DIMENSIONAL  MODEL  TERRAIN 


X  (KM) 


Fig.  22.  "Two-dimensionar  model  terrain  elevation  data.  Terrain  elevation.'!  are  rela¬ 
tive  to  lowest  terrain  elevation  value  in  slab  at  y  =  0.0  km  as  described  in  text.  Both 
axes  are  labeled  in  kilometers,  with  the  origin  at  Langmuir  Laboratoiy. 
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tains  at  0945  MST,  31  July  1984.  It  was  assumed  that  at  initialization  time  the  atmo¬ 
sphere  was  horizontally  homogeneous  throughout  the  model  domain.  While  some 
horizontal  variability  certainly  existed  and  probably  played  a  role  in  the  thunderstorm 
development,  it  was  neglected  in  this  study  since  it  was  not  possible  to  determine  the 
exact  nature  of  the  variability  from  the  available  (^servational  data.  The  base-state 
soil  moisture  and  temperature  profiles  (described  in  Chapter  5)  were  likewise  assumed 
to  be  horizontally  homogeneous  for  the  numerical  simulations. 

For  this  study,  no  special  initialization  procedures  were  employed.  The  model 
was  initialized  approximately  60  min  prior  to  the  expected  onset  of  convection  to  al¬ 
low  time  for  model  "spin-up."  The  employment  of  a  simplified  base-state  wind  pro¬ 
file,  described  later  in  Chapter  5,  combined  with  the  assumption  of  a  nearly 
dynamically  balanced  base-state  atmosphere  permitted  the  short  spin-up  time.  After 
an  initial  IS  min  period  designed  to  allow  dynamical  balancing  without  external  forc¬ 
ing  mechanisms,  the  short-wave  radiation  package  was  activated  to  produce  the  sur¬ 
face  heating  crucial  to  the  development  of  the  thunderstorms. 


Chapter  5 

Experimental  Design 

As  discussed  in  Chapter  1 ,  the  purpose  of  this  study  is  to  determine  the  important 
processes  involved  in  the  initiation,  development,  and  maintenance  of  mountain  thun¬ 
derstorms  by  comparing  the  results  of  numerical  simulation  experiments  with  Doppler 
radar  observations.  This  chapter  describes  the  experimental  design  chosen  to  achieve 
this  goal.  A  total  of  six  numerical  simulations  (one  control  and  five  experiments) 
were  conducted  with  the  RAMS  to  test  the  effect  of  differing  initial  conditions  and 
forcing  mechanisms.  Section  S.l  completes  the  description  of  the  basic  numerical 
model  configuration  chosen  for  the  control  simulation.  Section  S.2  describes  the  con¬ 
struction  of  the  base-state  atmospheric  environmental  profiles  used  to  initialize  the 
control  simulation  and  most  of  the  experiments.  Finally,  section  5.3  outlines  the  five 
numerical  simulation  experiments. 

5.1  Control  Simulation  Configuration 

This  section  completes  the  description  begun  in  Chapter  4  of  the  model  configura¬ 
tion  employed  in  the  control  simulation.  This  configuration  formed  the  basis  for  all 
the  numerical  experiments;  changes  to  the  configuration  for  those  experiments  are 
briefly  described  in  section  5.3  and  in  more  detail  in  later  chapters. 

The  simulation  initialization  time  and  location  were  carefully  chosen  to  pn^rly 
depict  the  radiative  heating  effects  on  the  surface.  The  control  simulation  was  ini¬ 
tialized  at  0945  MST  and  integrated  forward  in  time  9000  s  (2.5  hr).  The  model  do¬ 
main's  western  edge  was  placed  at  33.975®  N,  107.36“  W. 
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The  soil  model  assumes  a  uniform  soil  temperature  profile  throughout  the  depth  of 
the  model.  The  soil  temperature  is  taken  to  be  the  same  as  the  base-state  temperature 
at  the  lowest  level  in  the  atmospheric  profile.  The  base-state  soil  moisture  profile 
constructed  for  the  control  simulation  assumed  an  increasing  moisture  content  with 
increasing  depth  beneath  the  surface.  The  initial  soil  moisture  profile  was  specified  as 
a  percentage  of  the  soil  moisture  content  at  saturation  (expressed  in  terms  of  volume 
of  water  per  volume  of  soil)  (Tremback  and  Kessler,  1985).  Table  2  provides  the 
base-state  soil  moisture  profile  used  to  initialize  the  control  simulation.  A  sandy  clay 
loam  soil  type  was  chosen  for  use  in  the  soil  model.  This  choice  is  consistent  with  the 
soil  type  typically  found  in  the  mountains  of  the  southwestern  United  States. 


Table  2:  Soil  model  base-state  soil 
moisture  profile.  Soil  moisture  given  as 
percentage  of  saturated  soil  moisture. 


SOIL  MODEL 
LEVEL  (CM) 

SOIL  MOISTURE 
(%  OF  SATURATION) 

0 

36 

1 

41 

4 

48 

16 

56 

32 

60 

5.2  Atmospheric  Base-State  Profiles 

RAMS  requires  base-state  temperature,  moisture,  and  wind  profiles  for  initializa¬ 
tion.  This  section  describes  the  construction  of  the  atmospheric  base-state  profiles 


80 


used  to  initialize  most  of  this  study’s  numerical  experiments.  Modified  base-state 
profiles  used  to  initialize  some  of  the  simulation  experiments  will  be  discussed  as  ap¬ 
propriate  in  later  chapters. 

The  base-state  profiles  were  derived  from  the  observations  collected  over  and  near 
the  Magdalena  Mountains  on  the  morning  of  31  July  1984  (see  Chapter  3).  Data  gaps 
and  unreliable  wind  information  from  the  morning  soundings  prevented  their  use 
alone  to  initialize  the  model.  Instead,  the  two  rawinsonde  temperature  and  moisture 
soundings  were  combined  to  construct  the  base-state  temperature  and  moisture  pro¬ 
files.  The  alternative  wind  profiles  derived  from  synthesized  Doppler  radar  wind 
fields  (see  section  3.2.4)  were  used  to  construct  the  base-state  wind  profile. 

5.2.1  Temperature  And  Moisture  Profiles 

Figures  2  and  3  show  the  morning  radiosonde  temperature  and  moisture  soundings 
available  for  31  July  1984  (referred  to  as  the  airport  and  Langmuir  soundings  respec¬ 
tively).  The  Langmuir  sounding  was  collected  at  nearly  the  simulation  initialization 
time  but  does  not  extend  down  to  the  valley  floor.  The  airport  sounding  extends  to 
the  valley  floor  but  was  collected  about  1.5  hr  prior  to  initialization  time.  Therefore, 
neither  set  of  observed  temperature  and  moisture  data  could  be  used  by  itself  to  ini¬ 
tialize  the  model.  In  constructing  the  base-state  temperature  and  moisture  profiles, 
the  observational  data  was  preserved  to  die  greatest  extent  possible. 

The  airport  sounding  was  chosen  as  the  basis  for  the  base-state  temperature  pro¬ 
file.  The  later  sounding  (Langmuir)  was  not  used  explicitly  because  it  contained  un¬ 
usually  cold  temperatures  between  550  and  400  mb.  At  other  levels,  the  two 
soundings  were  nearly  identical.  The  airport  sounding  was  used  above  780  mb  with 
only  minor  smoothing.  Below  780  mb  the  observed  temperatures  were  increased 
slightly  to  account  for  the  effect  of  morning  heating  between  the  early  morning  ob- 
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servation  time  and  the  model  initialization  time.  The  bottom  of  the  profile  was  cho¬ 
sen  to  coincide  with  the  lowest  terrain  elevation  in  the  model  domain.  Figure  23 
shows  the  final  base-state  temperature  profile. 

Both  soundings  exhibited  basically  the  same  vertical  moisture  structure,  but  with 
widely  varying  differences  in  magnitude.  Above  660  mb  the  final  base-state  profile 
represents  a  smoothed  average  of  the  two  soundings.  Below  660  mb  the  mixing  ratio 
was  increased  by  about  l.S  g  kg  '  over  that  in  the  airport  sounding  to  account  for  the 
low-level  moistening  typically  observed  in  the  area  in  the  early  morning.  This  moist¬ 
ening,  described  by  Raymond  and  Wilkening  (1982),  is  brought  on  by  a  weak 
thermally-driven  upslope  flow  which  advects  moisture  from  the  lower  elevations  of 
the  Rio  Grande  Valley.  As  with  the  temperature  profile,  the  bottom  of  the  profile 
was  chosen  to  coincide  with  the  lowest  terrain  elevation  in  the  model  domain.  The 
final  moisture  profile  (Fig.  23)  is  similar  to  the  composite  thunderstorm  profile  de¬ 
scribed  by  Raymond  and  Wilkening  (1985). 

5.2.2  Wind  Profiles 

Figures  5  and  6  show  the  two  Doppler-derived  wind  profiles  available  for  31  July 
1984  (referred  to  as  Profiles  1  and  2,  respectively).  Profile  1  was  constructed  with 
data  collected  nearest  the  initialization  time  but  extends  to  only  approximately  400 
mb.  Profile  2  was  constructed  with  data  collected  well  after  the  initialization  time  but 
extends  to  at  least  200  mb.  Neither  profile  provides  wind  data  below  the  mountain 
top  (approximately  3280  m).  Figure  4  shows  the  Ziegler  (personal  communication, 
1986)  estimated  wind  profile  (referred  to  as  Profile  3)  that  extends  from  the  valley 
floor  to  2000  m  above  the  mountain-top  level  (approximately  5280  m  MSL).  Clearly 
neither  set  of  Doppler-derived  wind  profiles  nor  the  estimated  wind  profile  could  be 
used  alone  to  initialize  the  model. 


PRESSURE  (mb) 
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Fig.  23.  Skew  T,  log  P  plot  of  base-state  thermodynamic  profile  used  to  initialize  the 
numerical  model.  Heavy  solid  and  dashed  lines  represent  sensible  and  dewpoint  tem¬ 
perature  profiles,  respectively.  Skew  T,  log  P  diagram  labeled  as  in  Fig.  2. 


83 


This  subsection  describes  the  construction  of  a  preliminary  composite  wind  profile 
from  the  three  profiles  described  above.  It  then  describes  the  modifications  made  to 
the  preliminary  profile  to  arrive  at  the  base-state  wind  profile  used  to  initialize  the 
two-dimensional  (2-D)  numerical  simulations. 

5.2.2. 1  Preliiiiiiiary  Compo^te  Wind  Profile 

Since  Profile  1  contained  data  collected  nearest  the  initialization  time  and  before 
significant  storm  development,  it  was  used  wherever  possible  to  construct  the  prelimi¬ 
nary  wind  profile.  Profile  1  was  combined  with  Profile  2  (above  440  mb)  and  with 
Profile  3  (below  640  mb)  to  produce  the  preliminary  profile.  Figure  24  shows  the 
three  profiles  and  the  smoothed  preliminary  composite  wind  profile.  The  composite 
wind  profile  is  similar  to  the  composite  thunderstorm  profile  described  by  Raymond 
and  Wilkening  (1985). 

5.2.2.2  Base-State  Wind  Profiles 

The  preliminary  composite  wind  profile  served  as  the  basis  for  constructing  two 
base-state  wind  profiles:  one  a  full  profile  containing  both  u  and  v  wind  components 
(referred  to  as  the  3-D  profile),  and  the  other  a  profile  containing  only  the  u  compo¬ 
nent  (referred  to  as  the  2-D  profile).  Since  the  2-D  base-state  profile  was  based  di¬ 
rectly  on  the  3-D  profile,  the  3-D  profile  will  be  described  first. 

Above  480  mb  the  3-D  profile  is  identical  to  the  preliminary  composite  profile. 
Below  740  mb  (approximate  height  of  the  model  mountain  top)  the  wind  speed  was 
set  to  zero.  This  was  done  for  two  reasons.  First,  the  wind  profile  available  for  be¬ 
low  the  mountain  top  was  only  estimated— other  qualitative  observations  that  day  indi¬ 
cated  the  winds  were  light  and  variable  below  the  mountain  (Winn,  personal 
communication,  1987).  Second,  setting  the  wind  speed  to  zero  below  the  mountain 
top  significantly  simplified  the  model's  spin-up  phase.  Between  740  mb  and  480  mb 
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CX)MPOSITE  U  WIND  PROFILE 


WIND  SPEED  (M/S) 


Fig.  24a.  Composite  vertical  profile  of  u  wind  component  used  to  construct  profile  to 
initialize  the  numerical  model.  Heavy  solid  line  represents  preliminary  composite  pro¬ 
file  of  u  wind  component.  Observed  and  estimated  profiles  (shown  in  Figs.  4.  5.  and  6) 
are  shown  for  reference.  Axes  labeled  as  in  Fig.  4. 
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CX)MPOSITE  V  WIND  PROFILE 


WIND  SPEED  (M/S) 


Fig.  24h.  Composite  vertical  profile  of  v  wind  component  used  to  construct  profile  to 
initialize  the  numerical  model.  Heavy  solid  line  represents  preliminary  composite  pro¬ 
file  ofv  wind  component.  Observed  and  estimated  profiles  (shown  in  Figs.  4.  5.  and  6) 
are  shown  for  reference.  Axes  labeled  as  in  Fig.  4. 
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the  wind  speeds  were  gradually  increased  from  zero  at  740  mb  to  the  full  preliminary 
composite  profile  speed  at  480  mb.  Figure  25  shows  the  preliminary  profile  and  the 
final  3-D  base-state  profile. 

The  construction  of  the  2-D  base-state  profile  required  further  manipulation  to  sat¬ 
isfy  the  RAMS  requirements  for  2-D  simulations.  The  version  of  RAMS  employed 
for  this  study  allows  2-D  simulations  to  be  run  only  in  an  x-z  slab.  Thus  only  the  u 
component  of  the  wind  is  allowed  in  such  simulations.  The  2-D  profile  was  con¬ 
structed  from  the  3-D  base-state  profile  by  using  the  magnitude  of  the  total  wind  at 
each  pressure  level  as  the  -u  component  of  the  wind.  The  shape  of  the  mountain  and 
its  orientation  with  respect  to  the  predominant  wind  flow  permitted  this  approach. 
The  actual  wind  was  predominantly  southeasterly  (particularly  at  low-levels)  (Figs.  5, 
6,  and  25)  and  perpendicular  to  the  mountain  slopes.  Similar  wind  vs.  slope  relation¬ 
ships  can  be  produced  within  the  2-D  limitations  by  using  the  magnitude  of  the  actual 
wind  as  the  -u  component.  The  actual  wind  orientation  is  even  less  important  in  these 
simulations  since  the  magnitude  was  set  to  zero  below  the  mountain  top.  Figure  26 
shows  the  2-D  base-state  wind  profile. 

5.3  Numerical  Ebtperiments 

This  section  describes  the  five  numerical  simulation  experiments  conducted  to  test 
the  effect  of  differing  initial  conditions  and  forcing  mechanisms  on  the  simulated 
storms.  For  each  experiment,  only  one  such  condition  or  mechanism  was  varied  from 
the  control's  configuration.  The  simulation  experiments  will  be  referred  to  as  EXPl 
through  EXP5  for  the  rest  of  this  discussion.  Table  3  provides  a  brief  overview  of  all 
the  simulations  conducted  for  this  study.  The  term  "standard”  as  used  in  the  table  re¬ 
fers  to  the  profiles  or  configuration  used  for  the  control  simulation  and  described  in 
Chapter  4  and  in  Sections  5.1  and  5.2.  Details  of  base-state  profiles  differing  from 
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3-D  BASE-STATE  U  WIND  PROFILE 


WIND  SPEED  (M/S) 


Fig.  25a.  Vertical  profile  of  the  u  wind  component  of  the  "three-dimensional"  base- 
state  profile  (heavy  solid  line  with  ticks)  used  to  build  final  profile  for  input  to  the  mod¬ 
el.  The  preliminary  composite  profile  .shown  in  Fig.  24a  is  .show  i  for  reference.  Axes 
as  labeled  in  Fig.  4. 
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3-D  BASE-STATE  V  WIND  PROFILE 


WIND  SPEED  (M/S) 


Fig.  25h.  Vertical  profile  of  the  v  wind  component  of  the  "three-dimensional"  base- 
state  profile  (heavy  .solid  line  with  ticks)  used  to  build  final  profile  for  input  to  the  mod¬ 
el.  The  preliminary  composite  profile  shown  in  Fig.  24b  is  .shown  for  reference.  Axes 
as  labeled  in  Fig.  4. 


WIND  SPEED  (M/S) 


Fig.  26.  Vertical  wind  profile  used  to  initialize  the  numerical  model.  Wind  speed  in 
this  profile  equals  the  magnitude  of  the  total  wind  in  the  profiles  shown  in  Fig.  25. 
Axes  labeled  as  in  Fig  4. 
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Table  3:  Numerical  simulation  experiment  overview.  "Standard"  refers  to  condition 
or  process  employed  in  the  control  sinmilation. 


SIMULATION 

NAME 

WIND  PROFILE 

THERMODYNAMIC 

PROFILE 

MICROPHYSICS 

SHORT-WAVE 

RADIATION 

CONTROL 

STANDARD 

STANDARD 

STANDARD 

STANDARD 

EXPl 

CALM 

STANDARD 

STANDARD 

STANDARD 

EXP2 

50%  SPEED 
PROFILE 

STANDARD 

STANDARD 

STANDARD 

EXP3 

STANDARD 

STANDARD 

STANDARD 

NONE 

EXP4 

STANDARD 

STANDARD 

WARM  RAIN 
MICROPHYSICS 

STANDARD 

EXPS 

STANDARD 

MOISTENED 
LOWER  LEVELS 

STANDARD 

STANDARD 

those  of  the  control  simulation  will  be  described  as  appropriate  in  the  following  chap¬ 
ters, 

EXPl  and  EXP2  evaluated  the  sensitivity  of  the  simulated  storms  to  the  initial 
base-state  wind  profiles.  The  intent  behind  this  set  of  experiments  was  to  isolate  the 
contribution  of  the  environmental  wind  field  to  the  development  of  the  storms.  EXPl 
was  initialized  with  calm  winds  throughout  the  depth  of  the  atmosphere,  while 
EXP2's  initial  wind  profile  contained  wind  speeds  half  the  magnitude  of  the  control's. 

EXP3  and  EXP4  tested  the  sensitivity  of  the  simulations  to  diabatic  heating  ef¬ 
fects.  The  purpose  of  EXP3  was  to  evaluate  the  importance  of  solar  heating  in  storm 
development.  EXP3's  simulation  was  conducted  with  the  RAMS  short-wave  radiation 
package  deactivated;  in  effect,  "turning  off  the  sun”.  It  is  important  to  note  that  deac¬ 
tivation  of  the  short-wave  radiation  package  strongly  impacts  the  behavior  of  the  un¬ 
derlying  soil  model.  As  in  nature,  removing  the  model's  solar  radiation  also  cuts  off 
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surface  heating  and  any  thermal  forcing  resulting  from  such  heating.  The  purpose  of 
EXP4  was  to  determine  the  importance  of  the  i(%-phase  cloud  micio|Aysical  processes 
and  their  associated  latent  heating  effects  in  storm  develq>ment.  EXP4  was  con¬ 
ducted  with  all  the  ice-phase  microphysical  processes  deactivated  and  only  the  warm 
rain  microphysical  parameterizadons  activated. 

The  final  experiment,  EXP5,  was  conducted  to  evaluate  the  effect  on  the  simu¬ 
lated  storm  of  a  moistened  base-state  moisture  protile.  The  base-state  moisture  profile 
was  moistened  by  up  to  0.5  g  kg  *  below  720  mb.  Section  9.1  describes  the  profile 


used  for  EXP5. 


Chapter  6 

Analy^  And  Verification  Of  The  Control  Simulation 
This  chapter  describes  the  results  of  the  control  simulatitm  (referred  to  as  the  con¬ 
trol).  Section  6.1  provides  a  brief  analysis  of  the  simulated  storm's  evolution,  and 
Section  6.2  compares  the  simulated  storm  with  the  31  July  1984  Doppler  radar  ob¬ 
servations  described  in  Chapter  3.  Section  6.3  compares  the  microphysical  evolution 
of  the  simulated  storm  with  the  results  of  Lang's  (1991)  microphysical  retrieval  study 
of  a  similar  storm  that  occurred  on  3  August  1984. 

Throughout  the  simulation  descriptions  which  follow  in  this  and  later  chapters,  the 
following  conventions  will  be  employ^.  First,  all  time  references  will  be  simulation 
times  in  seconds  unless  otherwise  noted.  Table  4  provides  a  conversion  between 
simulation  times  in  seconds  and  hours  and  actual  times  with  the  0945  MST  model  ini¬ 
tialization  time  employed  for  these  simulations. 

All  references  lO  heights  will  be  based  on  the  model's  grid  coordinates  unless 
otherwise  noted.  As  discussed  in  Chapter  4,  the  model  terrain  elevations  were  ad¬ 
justed  so  that  the  lowest  terrain  elevation  in  the  model  grid  had  a  value  of  zero. 
Thus,  the  model's  lower  boundary  is  located  at  approximately  1820  m  MSL.  All  oth¬ 
er  model  grid  heights  are  reduced  from  actual  MSL  heights  by  this  same  amount.  Ta¬ 
ble  S  provides  a  conversion  between  model  heights  and  actual  MSL  heights. 

Finally,  for  the  purposes  of  the  following  descriptions,  a  cell  is  defined  as  any 
relative  condensate  and/or  updraft  maximum  that  exhibits  vertical  continuity.  This 
definition  follows  closely  that  used  by  Lang  (1991)  and  employed  in  Chapter  3. 
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Table  4:  Simulation  time  to  actual  time  conversion  chart. 


SIMULATION 

TLVKCi) 

SIMULATION 
TIME  (HR) 

ACTUAL  TIME 
(MST) 

SIMULATION 

T1ME(S) 

SIMULATHm 
TIME  (HR) 

ACTUAL  TINK 
(MST) 

0 

0.00 

9:45 

4680 

1.30 

11:03 

180 

0.05 

9:48 

4860 

1.35 

11:06 

360 

0.10 

9:51 

5040 

1.40 

11:09 

540 

0.15 

9:54 

5220 

1.45 

11:12 

720 

0.20 

9:67 

5400 

1.50 

11:15 

900 

0.25 

10:00 

6580 

1.55 

11:18 

1080 

0.30 

10:03 

5760 

1.60 

11:21 

1260 

0.35 

10:06 

5940 

1.65 

11:24 

1440 

0.40 

10:09 

6120 

1.70 

11:27 

1620 

0.46 

10:12 

6300 

1.75 

11:30 

1800 

0.60 

10:16 

6480 

1.80 

11:33 

1980 

0.55 

10:18 

6660 

1.85 

11:36 

2160 

0.60 

10:21 

6840 

1.90 

11:39 

2340 

0.65 

10:24 

7020 

1.95 

11:42 

2520 

0.70 

10:27 

7200 

2.00 

11:45 

2700 

0.75 

10:30 

7380 

2.05 

11:48 

2880 

0.80 

10:33 

7560 

2.10 

11:51 

3060 

0.85 

10:36 

7740 

2.15 

11:54 

3240 

0.90 

10:39 

7920 

2.20 

11:57 

3420 

0.95 

10:42 

8100 

2.25 

12K)0 

3600 

1.00 

10:46 

8280 

2.30 

12:03 

3780 

1.06 

10:48 

8460 

2.35 

12.06 

3960 

1.10 

10:51 

8640 

2.40 

12:09 

4140 

1.16 

10:54 

8820 

2.46 

12:12 

4320 

1.20 

10:57 

9000 

2.50 

12:15 

4500 

1.25 

11K)0 
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Table  5:  Model  height  to  actual  height  conversion. 


MODEL  HEIGHT 
(KM) 

ACTUAL  HEIGHT 
(KMMSL) 

0.00 

1.82 

0.25 

2.07 

0.75 

2.57 

1.25 

3.07 

1.75 

3.57 

2.25 

4.07 

2.75 

4.57 

3.25 

5.07 

3.75 

5.57 

4.25 

6.07 

4.75 

6.57 

5.25 

7.07 

5.75 

7.57 

6.25 

8.07 

6.75 

8.57 

7.25 

9.07 

7.75 

9.57 

8.25 

10.07 

8.75 

10.57 

9.25 

11.07 

9.75 

11.57 

10.25 

12.07 

10.75 

12.57 

11.25 

13.07 

11.75 

13.57 

12.25 

14.07 
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6.1  Analysis  of  Simulated  Storm  Evohition 

Three  primary  cells  developed  during  the  control  simulaticm.  The  origins  of  two 
of  the  cells  can  be  traced  back  to  the  earliest  cloud  and  updraft  development  in  the 
simulation.  The  third  cell  developed  primarily  as  a  result  of  the  upslope  wind  flow 
induced  by  circulation  of  the  first  two  cells. 

The  first  condensate  developed  at  5760  s  over  the  western  slope  of  the  mountain  at 
X  =  -2.5  km  and  z  =  2.0  km.  For  the  next  9  min,  the  cloud  grew  slightly  and  drifted 
westward.  By  6480  s  the  cloud  had  dissipated  and  was  replaced  by  a  new  cloud 
which  formed  at  nearly  the  same  position  as  the  first.  Figures  27  and  28  show  analy¬ 
ses  of  total  condensate  mixing  ratio  and  vertical  velocity,  respectively,  for  6480  s. 
An  analysis  of  the  wind  field  for  this  time  period  showed  that  the  first  cloud  drifted 
away  from  the  weak  updraft  which  contributed  to  its  initial  development.  The  updraft 
remained  fixed  at  approximately  its  original  position  and  led  to  the  subsequent  initia¬ 
tion  of  the  second  cloud  which  quickly  developed  into  a  cell.  During  this  same  peri¬ 
od,  a  low-level  westerly  flow  of  approximately  1.5  m  s  '  developed  below  2.0  km 
west  of  the  mountain  ridge.  Further  investigation  and  the  results  of  EXP3  (discussed 
in  Chapter  8)  show  that  this  westerly  flow  appeared  to  develop  as  a  result  of  a  dy¬ 
namically  induced  mountain-wave  circulation  on  the  lee  side  of  the  mountain.  The 
interaction  of  this  westerly  flow  with  the  easterly  flow  over  the  mountain  created  an 
area  of  low-level  convergence  at  af^roximately  x  =  -2.5  km  which  supported  the 
fixed  updraft  described  above.  Figure  29  shows  an  example  of  this  convergence. 

By  6660  s  the  first  cell  (denoted  hereafter  as  CCl)  was  located  at  x  =  -2.5  km  and 
z  =  2.5  and  was  approximately  1.0  km  deep.  A  new  cell  (denoted  hereafter  as  CC2) 
had  developed  just  over  the  ridge  at  x  ==  -1.0  km  and  z  =  2.0  km.  These  cells 
continued  to  develop  slowly  and  drift  slowly  westward  until  72(X)  s.  During  this  peri- 
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Fig.  27.  East- west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors  for  CONTROL  simulation  at  6480  .s.  Contours 
are  every  2.0  g  kg  '  starting  at  O.I  g  kg  '.  Wind  vector  length  shown  in  lower  right- 
hand  corner  of  figure  equals  50  m  s''.  Vertical  axes  is  labeled  in  kilometers  above  the 
minimum  terrain  elevation  in  the  simulation  grid  (in  this  case  1.82  km)  as  described  in 
the  text.  Horizontal  axis  is  horizontal  distance  labeled  in  kilometers.  The  cross-section 
shown  is  located  at  y  =  0.0  km  and  runs  east-west,  with  the  origin  at  Langmuir  Labora¬ 
tory. 
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Fig.  28.  East-west  vertical  cross-section  of  simulated  vertical  velocity  contoured  every 
3.0  m  s  ’  starting  at  1.5  m  s  ' for  CONTROL  simulation  at  6480  s.  Solid  (dashed)  con¬ 
tours  show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig.  27. 


98 


CONTROL  EXP 


WESTERLY  WIND 

Y  *  0.00  KM  T  =  6480.  S 


Fig.  29.  East-west  vertical  cross-section  of  simulated  u  wind  component  contoured  ev¬ 
ery  3.0  m  s  '  starting  at  1.5  m  s  '  for  CONTROL  simulation  at  6480  s.  Solid  (dashed^ 
contours  show  westerly  (easterly)  winds.  Axes  are  labeled  as  in  Fig.  27. 
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od,  the  convergent  area  remained  at  approximately  x  =  -2.5  km.  At  the  same  time,  a 
shallow  layer  of  easterly  flow  developed  beneath  the  westerly  flow,  creating  a  second 
low-level  convergent  area  atx  =  -5.0  km  (Fig.  30).  The  updraft  associated  with  CCl 
at  6660  s  separated  into  two  distinct  updrafts  after  6840  s.  The  upper  portion  of  the 
original  updraft  drifted  westward  with  the  condensate  while  the  lower  portion  of  the 
updraft  remained  anchored  at  approximately  x  =  -2.5  km.  CC2’s  original  updraft 
developed  at  about  6840  s  and  moved  slightly  westward  before  diminishing  after  7200 
s.  At  about  this  same  time,  the  condensate  associated  with  CCl  moved  over  the  an¬ 
chored  updraft  at.r  =  -2.5  km.  By  7200  s  both  major  updrafts  exceeded  4.5  m  s  '. 
Figures  3 1  and  32  show  analyses  of  total  condensate  mixing  ratio  and  vertical  veloc¬ 
ity,  respectively,  for  a  time  representative  of  this  period-6840  s. 

Between  7200  s  and  7560  s,  CC2's  condensate  drifted  westward  and  merged  with 
the  lower  remnants  of  CCl's  updraft.  During  this  period,  both  cells’  updrafts  re¬ 
mained  separate  and  drifted  westward.  CC2’s  updraft  increased  to  6.7  m  s  ',  while 
CCl’s  weakened  to  about  1.5  m  s'  by  7560  s.  At  7200  s,  a  new  updraft  developed 
west  of  CCl  atx  =  -5.0  km  in  response  to  the  secondary  convergent  area  at  that  loca¬ 
tion.  This  updraft  began  to  strengthen  as  CCl  moved  into  its  vicinity.  By  7560  s, 
both  CCl’s  original  updraft  and  the  new  updraft  were  approximately  the  same 
strength.  A  relatively  strong  easterly  jet  developed  west  of  CCl  at  z  =  3.75  km  near 
the  level  of  the  capping  inversion.  Several  small  cloud  masses  also  formed  over  the 
ridge  between  x  =  0.0  km  and  x  =  2.0  km  and  remained  very  weak  with  no  asso¬ 
ciated  updrafts.  The  first  significant  development  of  rain  water,  pristine  ice  crystals, 
and  graupel  occurred  during  this  period.  Prior  to  these  developments,  the  condensate 
consisted  almost  entirely  of  cloud  water.  Rain  water  first  developed  at  7100  s  in  CCl 
near  z  =  3.75  km.  By  7380  s,  most  of  the  small  area  of  rain  had  evaporated  as  it  fell 
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Fig.  30.  East-west  vertical  cross-section  of  simulated  u  wind  component  as  in  Fig.  29 
except  for  CONTROL  simulation  at  7200  s. 
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Fig.  31.  Analogous  to  Fig.  27,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  CONTROL  simulation  at 
6H40  s.  Contour  label  "I "  represents  O.l  g  kg  '. 
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Fig.  32.  Analogous  to  Fig.  28,  east-west  vertical  cross-section  of  simulated  vertical  ve¬ 
locity  for  CONTROL  simulation  at  6840  s. 
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out  of  CCl,  and  CCl  drifted  away.  Both  pristine  ice  crystals  and  graupel  developed 
in  CCl  near  z  =  4.0  km  between  71 10  and  7290  s.  Rain  water  developed  at  7560  s 
near  the  center  of  CC2  between  z  =  3.25  and  3.75  km.  Figures  33  and  34  show 
analyses  of  total  condensate  mixing  ratio  and  vertical  velocity,  respectively,  for 
7560  s. 

By  7740  s,  CC2  had  become  the  dominant  cell  in  the  merged  storm  with  maxi¬ 
mum  total  condensate  mixing  ratios  of  3.4  g  kg  '.  A  new  cell  (denoted  hereafter  as 
CC3)  developed  at  x  =  0.0  km  and  z  =  1.75  km  with  a  weak  updraft.  CC3's  devel¬ 
opment  appears  to  have  resulted  primarily  from  upslope  flow  induced  by  the  inflow 
into  CCl  and  CC2.  CCl's  original  updraft  dissipated  by  this  time  and  was  replaced 
by  the  newer  updraft  at  .t  =  -5.0  km  which  had  strengthened  to  4.5  m  s  '.  The  sec¬ 
ondary  area  of  low-level  convergence  continued  to  support  CCl’s  new  updraft  and 
also  provided  support  to  CC2's  updraft.  Graupel  first  developed  in  CC2  at  this  time 
between  z  =  3.0  and  4.5  km.  Ice  crystals  continued  in  CCl  and  also  began  to  devel¬ 
op  in  CC2  at  this  time. 

At  7920  s  the  first  precipitation  developed  as  a  rain  shaft  from  CC2  between  x  = 
-3.0  and  -4.0  km.  At  this  time  also,  CCl  had  regenerated  significantly  while  CC2 
weakened  slightly.  Maximum  total  condensate  mixing  ratios  in  both  CCl  and  CC2 
ranged  from  2.0  to  3.0  g  kg  '.  The  rapid  regeneration  of  CCl  can  be  attributed  to  the 
rapid  strengthening  of  its  new  updraft  to  over  10.0  m  s  '.  At  the  same  time,  another 
new  updraft  also  developed  just  to  the  west  of  CCl  in  response  to  the  low-level  con¬ 
vergence  which  continued  near  x  =  -6.0  km.  Graupel  and  rain  water  redeveloped  in 
CCl  between  z  =  2.25  and  3.75  km.  CCl's  total  condensate  maximum  consisted  en¬ 
tirely  of  cloud  water,  while  CC2's  maximum  consisted  of  a  mixture  of  ice  crystals. 
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Fig.  33.  As  in  Fig.  27,  east-west  vertical  cross-section  of  simulated  total  condensate 
mixing  ratios  overlaid  with  simulated  wind  vectors  for  CONTROL  simulation  at  7560  s. 
Contour  label  "I "  represents  0.1  g  kg  '. 
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Fig  34.  .,4a’  in  Fig.  28,  east-west  vertical  cross-section  of  simulated  vertical  velocity  for 
CONTROL  simulation  at  7560  s. 
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cloud  water,  and  graupel.  Figures  35  and  36  show  analyses  of  total  condensate  mix¬ 
ing  ratio  and  vertical  velocity,  respectively,  for  7920  s. 

CCl  had  become  the  dominant  cell  at  8100  s  with  a  maximum  total  condensate 
mixing  ratio  exceeding  4.0  g  kg By  this  time  both  updrafts  in  the  vicinity  of  CCl 
had  merged  to  produce  a  double  maxima  updraft.  The  upper  portion  of  the  updraft 
had  strengthened  to  13.0  m  s  '  while  the  lower  portion  exceeded  4.5  m  s  '.  Both  por¬ 
tions  of  the  updraft  continued  to  be  suf^rted  by  the  low-level  convergence  near  x  = 
-6.0  km.  CC2's  updraft  remained  distinct  with  a  maximum  strength  of  over  4.5 
ms'.  A  significant  downdraft  of  5.5  m  s'  developed  at  this  time  between  the  up¬ 
drafts  associated  with  CCl  and  CC2  and  at  a  height  of  z  =  3.75  km,  CC3  remained 
very  weak  with  several  weak  updrafts  in  its  vicinity.  A  new  area  of  rain  water  devel¬ 
oped  near  CCl 's  maximum  total  condensate  at  z  =  4.25  km.  CCl's  total  condensate 
maximum  at  this  time  consisted  of  a  mixture  of  ice  crystals,  cloud  water,  rain  water, 
and  graupel.  Figures  37  and  38  show  analyses  of  total  condensate  mixing  ratio  and 
vertical  velocity,  respectively,  for  8100  s. 

By  8280  s,  CCl  was  clearly  the  dominant  cell  with  double  maxima  in  both  the 
total  condensate  mixing  ratio  and  the  vertical  velocity  fields.  Maximum  total  conden¬ 
sate  mixing  ratios  exceeded  4.0  g  kg  '  in  the  upper  portion  of  the  cell  and  2.0  g  kg  ' 
in  the  lower  portion.  The  cell’s  condensate  reached  a  maximum  height  of  about  6.0 
km.  The  upper  portion  of  the  updraft  remained  strtmg  at  11.0  m  s  '  while  the  lower 
portion  strengthened  to  over  7.5  m  s'*.  The  low-level  convergence  continued  beneath 
CCl's  updraft  while  a  divergent  flow  developed  at  z  =  5.25  km  near  the  top  of  the 
updraft.  As  rain  began  to  fall  from  CCl,  the  storm's  rain  shaft  increased  in  size  to 
approximately  2.5  km  wide.  CC2  remained  evident  in  the  total  condensate  analysis  as 
a  turret  of  condensate  on  the  eastern  side  of  the  storm.  CC2's  updraft  diminished 


IKH) 


107 


CONTROL  EXP 


X  (KM) 

TOTAL  CONDENSATE 

Y  =  0.00  KM  T  »  7920.  S 


WINDS 

IN  VECTORS 


IMXtIMI 


Fig.  35.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  av  in  Fig.  27.  except  for  CONTROL  simulation  at 
7920  s.  Contour  label  "I "  represents  0.1  g  kg'. 
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Fig.  36.  East-west  vertical  cross-section  of  simulated  vertical  velocity,  as  in  Fig.  28, 
except  for  CONTROL  simulation  at  7920  s. 
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Fig.  37.  Analogous  to  Fig.  27.  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  CONTROL  simulation  at 
H 1 00  s.  Contour  label  "I "  represents  0.1  g  kg'. 
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Fig.  38.  Analogous  to  Fig.  28,  east- west  vertical  cross-section  of  simulated  vertical  ve¬ 
locity  for  CONTROL  simulation  at  8100  s. 
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considerably  in  size  from  the  previous  analysis.  CC3  had  increased  in  size  (primarily 
by  building  westward),  but  still  contained  no  significant  updrafts.  Aggregated  snow 
flakes  first  developed  at  this  time  in  CC2.  Both  CCl  and  CC2  contained  well-defined 
graupel  shafts.  CC3‘s  condensate  at  this  time  consisted  entirely  of  cloud  water. 
CC2's  area  of  maximum  total  condensate  consisted  primarily  of  pristine  ice  crystals 
but  also  contained  a  mixture  of  graupel,  cloud  water,  and  aggregated  snow  flakes. 
CCl's  maxima  on  the  other  hand,  consisted  primarily  of  cloud  water,  with  smaller 
amounts  of  rain  water,  graupel,  and  pristine  ice  crystals. 

At  8460  s,  CCl  continued  to  strengthen  with  maximum  total  condensate  mixing 
ratios  exceeding  6.0  g  kg  ‘  at  approximately  6.0  km  near  the  top  of  the  cell.  CCl's 
updraft  again  contained  a  double  maxima  with  vertical  velocities  exceeding  10.5  m  s  ' 
and  14.5  m  s  '  in  the  upper  and  lower  maxima,  respectively.  By  this  time,  the  area  of 
low-level  convergence  supporting  CCl's  updraft  had  moved  slightly  westward  to  v  = 
-7.5  km.  This  movement  resulted  from  the  westward  expansion  of  the  cool  easterly 
outflow  associated  with  the  precipitation  shaft.  Meanwhile  the  divergent  flow  contin¬ 
ued  near  the  top  of  CCl's  updraft  at  z  =  5.75  km.  CCl's  condensate  reached  a 
maximum  height  of  about  7.5  km  at  this  time.  CC2  remained  distinct  as  a  turret  of 
condensate  on  the  eastern  edge  of  the  storm,  however,  its  updraft  had  weakened  to 
about  1.5  m  s  '.  At  the  same  time,  CC3  had  strengthened  significantly  and  moved 
slightly  westward.  CC3's  maximum  total  condensate  mixing  ratios  had  increased  to 
2.0  g  kg  ',  and  its  updraft  developed  to  a  maximum  of  4.5  m  s  '.  Both  CCl  and  CC2 
were  dominated  by  the  ice  phase  at  this  time.  The  areas  of  highest  total  condensate 
mixing  ratios  in  both  cells  consisted  primarily  of  ice  crystals,  while  very  little  cloud 
water  existed  anywhere  in  the  cells.  Aggregated  snow  flakes  had  also  developed  by 
this  time  in  CCl.  As  expected  from  the  location  of  the  pristine  ice  crystals,  the  ag- 
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gregates  were  located  near  the  total  condensate  maxima  in  both  cells.  Graupel  also 
continued  in  both  CCl  and  CC2.  Rain  water  developed  for  the  first  time  in  CC3. 

CCl  continued  to  dominate  the  simulation  at  8640  s.  While  CCl's  maximum  to¬ 
tal  condensate  mixing  ratios  decreased  slightly  to  S.9  g  kg'‘,  its  updraft  continued  to 
exhibit  a  double  maxima  with  vertical  velocities  remaining  about  the  same  strength  as 
at  8460  s.  The  low-level  convergence  at  the  outflow  boundary  continued  to  support 
CCl's  updraft  at  this  time,  while  the  divergence  near  the  top  of  the  updraft  continued 
to  increase  in  height  to  approximately  6.5  km.  The  downdraft  located  just  to  the  east 
of  CCl's  updraft  strengthened  again  to  4.5  m  s  ‘.  CCl's  condensate  reached  a  maxi¬ 
mum  height  of  approximately  9.0  km.  A  new  updraft  with  maximum  vertical  veloci¬ 
ties  of  4.5  m  s  '  developed  at  this  time  between  the  condensate  masses  associated  with 
CCl  and  the  remnants  of  CC2,  CC3  continued  to  strengthen  and  move  slowly  west¬ 
ward  to  about  X  =  -1.0  km.  While  CC3's  total  condensate  mixing  ratios  remained 
virtually  unchanged  from  8460  s,  its  updraft  increased  to  a  maximum  of  7.5  ms'.  A 
very  weak  westerly  flow,  resulting  from  the  precipitation-induced  outflow,  developed 
at  this  time  nearx  =  -4.5  km  and  z  =  1.5  km.  The  interaction  of  this  westerly  flow 
with  the  prevailing  easterly  flow  at  that  location  produced  another  area  of  conver¬ 
gence  which  helped  strengthen  CC3's  updraft.  The  ice  phase  continued  to  dominate 
CCl  at  this  time.  The  bulk  of  CCl 's  condensate,  particularly  near  its  maximum,  con¬ 
sisted  of  pristine  ice  crystals.  At  the  same  time,  a  deep  shaft  of  graupel  had  devel¬ 
oped  near  the  center  the  cell.  The  areas  of  aggregates  continued  near  the  ice 
crystal  maxima  in  both  CCl  and  CC2.  The  bulk  of  CC3's  condensate  remained  cloud 
water  with  smaller  amounts  of  graupel  and  rain  water.  Figure  39  shows  analyses  of 
total  condensate,  cloud  water,  rain  water,  pristine  ice  crystals,  graupel,  and  £^re- 
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Fig.  39a.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  .simulated  wind  vectors,  as  in  Fig.  27,  except  for  CONTROL  simulation  at 
H64()  s.  Contour  label  "I "  represents  O.l  g  kg  '. 
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Fig.  39h.  East-west  vertical  cross-section  of  simulated  cloud  water  mixing  ratios  con¬ 
toured  every  2.0  g  kg  starting  at  0.1  g  kg '  for  CONTROL  simulation  at  8640  s.  Axes 
labeled  as  in  Fig.  27 .  Contour  label  ”/"  represents  0. 1  g  kg  '. 
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Fig.  39c.  Analogous  to  Fig.  39b,  east-west  vertical  cross-section  of  simulated  pain  wa¬ 
ter  mixing  ratios  for  CONTROL  simulation  at  8640  .v.  Contour  label  "I  ”  represents  0.  / 
gkg'. 
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Fig.  39d.  Analogous  to  Fig.  39b,  east-west  vertical  cross-section  of  simulated  pristine 
ice  crystal  mixing  ratios  for  CONTROL  .simulation  at  8640  s.  Contour  label  repre¬ 
sents  O.I  g  kg'. 
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Fig.  39e.  Analogous  to  Fig.  39b,  east-west  vertical  cross-section  of  simulated  graupel 
mixing  ratios  for  CONTROL  simulation  at  8640  s.  Contour  label  "I "  represents  0. 1 
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Fig.  39f.  Analogous  to  Fig.  39b,  east-west  vertical  cross-section  of  simulated  aggre¬ 
gated  snowflake  mixing  ratios  for  CONTROL  simulation  at  8640  s.  Conhnir  label  "I " 
represents  0. 1  g  kg  '. 
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gated  snow  flake  mixing  ratios,  reflectively,  for  8640  s.  Figure  40  shows  an  analysis 
of  vertical  velocity  for  8640  s. 

By  8820  s  CC2  had  lost  its  identity  as  CC3  moved  westward  and  CC3's  conden¬ 
sate  merged  with  CCl.  The  maximum  total  condensate  mixing  ratio  excee<kd  S.9 
g  kg  '  in  both  CCl  and  CC3,  with  the  highest  mixing  ratios  located  above  4.0  km  in 
both  cells.  CCl's  condensate  reached  a  maximum  height  of  10.25  km,  and  precipita¬ 
tion  first  developed  from  CC3  at  this  time.  CCl's  updraft  contained  several  local 
vertical  velocity  maxima  with  the  primary  maximum  of  16.6  m  s  ‘  located  at  8.0  km 
near  the  top  of  the  cell.  The  precipitation-induced  outflow  below  CCl  continued  to 
create  areas  of  convergence  both  to  the  west  and  east  of  the  main  precipitation  shaft. 
The  westerly  flow  portion  of  the  outflow  first  noted  at  8640  s,  became  better  defined 
by  this  time.  The  divergent  flow  near  CCl’s  updraft  extended  to  nearly  z  =  8.0  km. 
The  downdraft  on  the  eastern  edge  of  CCl  weakened  considerably,  while  a  new 
downdraft  of  about  6.2  m  s  ‘  formed  near  the  center  of  the  precipitation  shaft.  CC3’s 
updraft  continued  at  its  previous  strength  of  7.5  m  s  *,  but  the  updraft's  structure  be¬ 
came  much  less  organized.  A  weak  divergent  flow  developed  at  4.75  km  near  the  top 
of  CC3's  updraft.  The  new  updraft  which  had  developed  between  CCl  and  CC2  at 
8640  s  weakened  dramatically  to  1 .5  m  s’*  as  it  came  under  the  influence  of  the  down- 
drafts  from  CCl  and  CC3  as  they  merged.  By  this  time,  the  bulk  of  CCl’s  conden¬ 
sate  consisted  of  either  pristine  ice  crystals  or  aggregated  snow  flakes.  Aggregate 
mixing  ratios  exceeded  2.0  g  kg'*.  The  graupel  shaft  continued  in  CCl,  and  graupel 
reached  the  surface  for  the  first  time.  CC3’s  condensate  still  consisted  primarily  of 
cloud  water,  but  pristine  ice  crystals  had  developed  for  the  first  time.  The  area  of 
graupel  expanded,  and  rain  reached  the  surface  for  the  first  time  out  of  CC3.  Figures 
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Fig.  40.  East-west  vertical  cross-section  of  simulated  vertical  velocity  as  in  Fig.  28 
except  for  CONTROL  simulation  at  8640  s. 
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41  and  42  show  analyses  of  total  condensate  mixing  ratio  ami  vertical  velocity,  re¬ 
spectively,  for  8820  s. 

At  9000  s,  CCl  still  dominated  the  merged  storm.  The  maximum  total  crxidensate 
mixing  ratio  exceeded  6.0  g  kg'*  in  CCl,  while  barely  reaching  2.0  g  kg'*  in  CC3. 
CCl's  condensate  reached  a  maximum  height  of  11.0  km.  The  updraft  structure  be¬ 
came  very  complicated  as  CCl  and  CC3  began  to  interact  with  each  otlwr.  CCl's 
main  updraft  contained  a  maximum  vertical  velocity  of  over  lO.S  m  s'*.  Meanwhile  a 
secondary  updraft  core  had  developed  slightly  to  the  east  of  the  original.  This  core 
contained  a  maximum  vertical  velocity  of  7.5  m  s  '.  At  the  same  time,  CC3's  updraft 
had  become  more  organized  and  strengthened  to  10.5  m  s  '.  By  this  time  the  low- 
level  outflow  convergence  areas  were  very  well  defined.  The  convergent  flow  on  the 
western  side  of  the  outflow  spawned  the  development  of  several  updrafts  on  the  west¬ 
ern  edge  of  CCl,  while  the  strong  low-level  convergence  on  the  eastern  side  of  the 
outflow  supported  the  strengthening  of  CC3's  updraft.  The  uj^r-level  divergence 
near  CCl's  updraft  continued,  but  became  less  organized  than  at  previous  times.  The 
same  basic  microphysical  structure  existed  in  the  cells  at  this  time  as  reported  for  the 
previous  time. 

6.2  Comparison  With  Radar  Observations 

This  section  presents  a  comparison  of  the  control  simulation  results  with  the  Dop¬ 
pler  radar  storm  observations  from  the  31  July  1984  case.  The  two-dimensional  na¬ 
ture  of  the  numerical  simulations  conducted  for  this  study  naturally  limit  direct 
comparisons  with  storm  observations,  however,  qualitative  and  some  quantitative 
comparisons  can  yield  a  realistic  assessment  of  the  simulation's  accuracy.  Establish¬ 
ing  the  fidelity  of  the  control  simulation  represents  a  crucial  first  step  of  this  study 
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Fig.  41.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  as  in  Fig.  27,  except  for  CONTROL  simulation  at 
HH20  s.  Contour  label  "I "  represents  0. 1  g  kg  '. 
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Fig.  42.  East-west  vertical  cross-section  of  simulated  vertical  velocity,  as  in  Fig.  27. 
except  for  CONTROL  simulation  at  8820  s. 
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since  the  numerical  simulation  experiments  discussed  in  later  chapters  will  be  com¬ 
pared  to  the  control. 

Given  the  limitations  imposed  by  the  two-dimensionality  of  the  numerical  simu¬ 
lation,  the  simulation  compared  very  favorably  with  the  Doppler  radar  observations 
collected  on  31  July  1984.  The  simulaticm  captured  most  of  the  important  cell  and 
storm  features  and  characteristics  noted  in  the  description  of  the  radar  observations  in 
Chapter  3. 

6.2.1  General  Comparisons 

The  control  simulation  captured  the  multicellular  nature  of  the  storms  as  observed 
by  the  radars.  Three  dominant  cells  developed  in  the  simulation  while  seven  cells 
were  observed  by  the  radars.  Some  of  this  difference  can  be  attributed  to  the  two- 
dimensional  nature  of  the  simulation  which  naturally  decreases  the  amount  of  space 
available  for  cell  development.  Second,  as  described  later,  the  convective  develop¬ 
ment  in  the  simulation  lags  the  observations  by  between  15  and  40  min.  The  final 
four  observed  cells  developed  in  the  last  18  min  of  the  analysis.  It  is  likely  that  the 
simulation  lag  resulted  in  the  simulation  ending  before  further  cell  development  oc¬ 
curred. 

The  development  of  low-level  convergence  areas  in  the  simulation  compares  rea¬ 
sonably  well  with  that  seen  in  the  observtuions.  Since  the  simulation  was  conducted 
in  two  dimensions,  it  could  not  capture  the  terrain-induced  channelling  convergence 
found  in  the  observations.  The  simulation  did  show  the  development  of  low-level 
convergence  areas  in  response  to  the  mountain-valley  flow  as  observed  in  the  synthe¬ 
sized  wind  fields.  In  both  the  simulation  and  in  the  observations,  the  low-level  con¬ 
vergence  areas  were  crucial  to  convective  initiation.  These  convergence  areas 
primarily  determined  the  location  of  initiation.  The  simulation  also  reasonably  repro- 
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duced  the  development  of  upper-level  divergence  areas  above  the  strongest  cells.  For 
the  most  part,  the  simulation's  upper-level  divergence  areas  were  weaker  than  those  in 
the  synthesized  wind  fields. 

In  both  the  observations  and  in  the  simulation,  the  convective  initiation  occurred 
over  or  very  near  the  ridges.  The  simulated  initiation  locations  tended  to  be  slightly 
to  the  west  of  the  ridge  lines  while  in  the  observations,  the  initiation  locations 
generally  were  either  directly  over  or  slightly  to  the  east  of  the  ridge  lines.  The  mag¬ 
nitude  of  the  distance  differences  between  observed  and  simulation  initiation  locations 
was  always  less  than  5.0  km.  The  two-dimensional  nature  of  the  simulation  again  ex¬ 
plains  in  part  the  difference  in  the  pattern  of  initiation  locations.  As  described  in 
Chapter  3,  the  initiation  locations  were  strongly  influenced  by  the  channelling  conver¬ 
gence  induced  by  the  narrow  valleys.  These  features  could  not  be  adequately  repre¬ 
sented  in  a  two-dimensional  simulation.  As  a  result,  the  simulation's  initiation 
locations  are  more  strongly  determined  by  the  thermally-driven  mountain-valley 
flows. 

As  mentioned  previously,  the  simulation's  convective  development  lagged  that  ob¬ 
served  by  the  radars  by  a  significant  amount  of  time.  The  two-dimensional  nature  of 
the  simulation  makes  it  difficult  to  detemune  the  exact  timing  difference.-  As  ex¬ 
pected,  investigation  of  individual  features  yields  different  timing  differences.  For 
example,  a  comparison  of  the  simulated  maximum  updrafts  and  downdrafts  with  those 
synthesized  from  the  radar  observations  (Fig.  43)  shows  good  agreement  between  the 
simulation  and  the  observations  with  a  simulation  lag  time  of  about  14  min.  The 
simulation's  condensate  first  developed  at  1115  (simulation  time)  and  the  condensate 
first  extended  down  to  1.5  km  (3.3  km  MSL)  at  1151  (simulation  time).  Comparing 
this  time  with  the  observed  radar  reflectivities  of  30  dBZ  at  3.3  km  MSL  at  1116 
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SYNTHESIZED  AND  SMULATED  VERTICAL  VELOCITIES 
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Fig.  43.  Time  evolution  of  synthesized  (heavy  curves  with  ticks)  and  CONTROL  .simu¬ 
lation  simulated  (light  curves  without  ticks)  maximum  and  minimum  vertical  velocities. 
Solid  curves  show  maximum  upward  vertical  velocities  (positive  w).  Dashed  cunes 
show  maximum  downward  vertical  velocities  (negative  w).  Vertical  axis  is  veiiical  ve¬ 
locity  in  m  .s  ',  bottom  horizontal  axis  is  simulation  time  in  seconds,  and  top  horizontal 
axis  Ls  observation  time  in  hours  and  minutes  (MST). 
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MST  yields  a  lag  time  of  at  least  35  min.  It  is  important,  therefore,  to  exercise  cau¬ 
tion  when  attempting  to  make  one-to-one  temporal  comparisons  between  the  simulated 
storms  and  the  observed  storms. 

The  simulated  updrafts  qualitatively  compared  quite  well  with  those  analyzed  from 
synthesized  wind  fields.  The  updrafts  in  both  cases  agreed  quite  well  in  both  horizon¬ 
tal  and  vertical  extent.  The  width  of  the  simulated  updrafts  averaged  only  slightly  less 
than  the  2.0  km  widths  seen  in  the  observations.  The  simulated  updrafts  were  about 
8.0  km  deep  compared  to  the  7.0  to  8.0  km  depth  of  the  observed  updrafts.  The 
comparison  of  the  maximum  updrafts  (again  using  Figure  43)  shows  that,  taking  into 
consideration  the  timing  lag  discussed  above,  the  maximum  strength  of  the  simulated 
updrafts  agrees  very  well  with  the  observed  maximum  updrafts.  In  general,  the  de¬ 
velopment  of  the  simulated  updrafts  also  follows  quite  well  that  seen  in  the  synthe¬ 
sized  wind  fields.  In  both  cases  cells  frequently  contained  single  updrafts  with 
vertical  velocities  of  12.0  m  s  '  or  less.  Also  in  both  cases,  early  in  the  convective 
development,  the  updrafts  are  strongest  in  the  low-  to  mid-levels,  and  both  show  a 
tendency  for  the  updraft  maxima  to  move  to  the  mid-  to  upper-levels  as  the  cells  ma¬ 
ture.  The  simulation  showed  that  some  cells  developed  either  multiple  updrafts  or 
updrafts  with  two  or  more  local  maxima  in  the  vertical  direction.  The  observed 
storms  also  showed  a  similar  phenomenon.  The  simulated  updrafts  for  the  most  part 
are  nearly  vertical  with  very  little  tilt.  This  contrasts  markedly  with  the  observed  up¬ 
drafts  which  frequently  tilted  backward  (away  ft'om  the  inflow).  As  observed  in  the 
synthesized  wind  fields,  the  simulated  updrafts  developed  most  often  near  regions  of 
low-level  convergence.  The  stronger  updrafts  were  associated  also  with  regions  of  at 
least  weak  upper-level  divergence. 


128 


The  simulated  downdrafts  do  not  compare  nearly  as  well  with  observations.  A 
comparison  of  the  horizontal  and  vertical  extent  of  the  simulated  downdrafts  with  the 
observed  shows  only  a  poor  correspondence.  The  simulated  downdrafts  for  the  most 
part  are  much  narrower  and  shallower  than  those  observed.  A  comparison  of  the 
maximum  downdrafts  (again  using  Figure  43)  shows  agreement  between  both  the 
simulation  and  the  observations  after  taking  into  consideration  the  simulation's  timing 
lag.  The  simulation  fails  to  develop  the  strong  downdrafts  of  over  9.0  m  s  '  seen  in 
the  observations  between  1158  and  1201  MST.  The  simulation  also  fails  to  capture 
the  broad  areas  of  strong  downdrafts  as  observed.  Instead  the  simulated  downdrafts 
were  usually  confined  to  very  narrow  regions.  The  simulated  downdrafts  were  lo¬ 
cated  primarily  in  low-levels  as  was  predominantly  the  case  with  the  observations. 
However,  the  simulated  downdrafts  could  also  be  found  throughout  nearly  the  entire 
vertical  extent  of  the  model  domain. 

The  analysis  of  the  validity  of  the  simulation's  precipitation  development  is  some¬ 
what  hindered  by  the  fact  that  the  lowest  radar  analyses  were  performed  at  3.3  km 
MSL.  Thus  it  is  not  possible  to  determine  precisely  where  or  when  the  precipitation 
was  associated  with  radar  returns  at  that  level.  If  it  is  assumed  that  the  precipitation 
shaft  is  approximately  the  size  of  the  radar  maximum  at  the  3.3  km  MSL  analysis  lev¬ 
el,  then  the  simulation's  precipitation  shafts  compared  very  favorably  with  those  ob¬ 
served.  In  both  cases,  the  precipitation  shafts  early  in  the  precipitating  stage  ranged 
from  1 .0  to  1 .5  km  in  diameter  and  increased  to  2.5  to  3.0  km  in  the  latter  part  of  the 
mature  stage.  Also  in  both  cases,  several  precipitation  shafts  developed.  The  ob¬ 
servations  contained  more  precipitation  shafts,  and  they  were  farther  apart  than  in  the 
simulation;  a  discrepancy  most  likely  resulting  from  the  limitations  imposed  by  the 
two-dimensional  nature  of  the  simulation.  In  both  the  simulation  and  in  the  observa- 
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tions,  the  precipitation  shafts  were  located  near  or  directly  beneath  major  updrafts. 
The  separation  between  precipitation  and  updrafts  was  generally  greater  in  the  ob¬ 
servations  than  in  the  simulation. 

The  simulation's  ability  to  capture  the  vertical  development  of  cells  can  be  as¬ 
sessed  by  comparing  the  height  of  simulated  condensate  with  the  height  of  radar  re¬ 
turns.  A  comparison  of  the  maximum  heights  of  the  simulated  condensate  with  the 
maximum  height  of  radar  returns  (Fig.  44)  shows  a  lag  of  approximately  23  min  be¬ 
tween  the  simulated  development  and  the  actual  development  similar  to  that  discussed 
earlier  in  this  section.  Taking  this  lag  into  consideration,  the  simulation  captured  the 
rapid  increase  in  maximum  radar  return  height  observed  between  1 146  MST  and  1 152 
MST.  Through  most  of  the  simulation,  the  maximum  condensate  heights  remained 
2.0  to  3.0  km  lower  than  the  observed  radar  return  heights.  The  simulation  indicated 
a  slow  steady  growth  in  height  to  about  5.2S  km  (7.0  km  MSL)  followed  by  a  single 
rapid  surge  in  height  to  near  the  absolute  maximum  seen  in  the  radar  data.  The  rapid 
development  which  followed  the  penetration  of  the  5.25  km  (7.0  km  MSL)  level  ap¬ 
pears  to  result  from  the  favorable  moisture  and  temperature  profiles  above  that  level. 
The  relatively  strong  inversion  and  the  very  dry  layer  which  exM,ted  below  5.25  km 
(7.0  km  MSL)  capped  the  convective  development  until  late  in  the  simulation.  A 
similar  capping  phenomenon  can  be  seen  in  the  observed  early  convective  develop¬ 
ment  where  no  echoes  penetrated  above  6.5  km  MSL  until  after  11  ;.3  MST.  The 
simulation  failed  to  capture  the  observed  double  peak  in  maximum  heights  seen  in 
Figure  44. 

The  lifetime  of  the  simulated  cells  compared  very  closely  with  those  observed  in 
the  radar  observations.  The  longest-lived  simulated  cell,  CCl,  existed  for  42  min 
prior  to  the  simulation's  termination.  This  cell  likely  would  have  existed  for  up  to 
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OBSERVED  AND  SIMULATED  CLOUD  TOP  HEIGHTS 

OBSERVATION  TIME  (MST) 


SIMULATION  TIME  (SEC) 


Fig.  44.  Time-height  cross  section  showing  maximum  simulated  total  condemate 
height  for  the  CONTROL  simulation  (plain  curve)  and  maximum  height  of  observed  re¬ 
flectivity  (solid  curve  with  ticks).  Plain  curve  shows  maximum  height  of  0.1  g  kg  '  total 
condensate  mixing  ratios.  Ticked  curve  shows  maximum  height  of  10  dBZ  radar  reflec¬ 
tivities.  Right  vertical  axis  is  labeled  in  kilometers  above  the  minimum  terrain  elevation 
in  the  simulation  grid  (in  this  case  1.82  km)  as  described  in  the  text,  and  l^  vertical 
axis  is  height  relative  to  mean  sea  level  in  km.  Bottom  horizontal  axis  is  simulation 
time  in  seconds  and  top  horizontal  axis  is  observation  time  in  hours  and  minutes  MST. 
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another  10  tnin  had  the  simulation  been  extended.  CC2  existed  for  36  min  before  be¬ 
ing  absorbed  in  CCl.  Finally,  CC3  existed  only  21  min  before  the  simulation  ended. 
This  cell  was  in  the  early  stages  of  development  and  likely  would  have  developed  for 
some  time  had  the  simulation  been  extended.  These  lifetimes  compare  well  with  the 
20  to  60  min  lifetimes  seen  in  the  observations.  It  is  likely  that  had  the  simulation 
been  extended,  the  storm  would  have  grown  and  eventually  decayed  with  the  result 
that  the  overall  storm  lifetime  also  would  have  compared  favorably  with  the  1.0  to  1.5 
hr  storm  lifetime  observed  by  the  radars. 

Only  one  of  the  simulated  cells,  CCl,  showed  the  tendency  for  periodic  reintensi¬ 
fication  seen  in  the  observed  cells.  In  the  simulated  cell,  the  cycle  of  reintensification 
occurred  early  in  the  cell's  development  between  7200  and  7740  s  simulation  time. 
The  simulated  cell's  reintensification  period  of  9  min  agreed  well  with  the  observed 
reintensification  periods  of  6  to  12  min. 

Finally,  the  simulation  captured  reasonably  well  the  intercellular  relationships  seen 
in  the  radar  data  and  the  synthesized  wind  fields.  The  simulation  accurately  repro¬ 
duced  the  alternating  intensification  between  cells  described  in  Chapter  3.  CC2  ini¬ 
tially  developed  as  the  strongest  cell  in  the  simulation.  Then  by  8100  s,  CCl  had 
become  the  dominant  cell.  Later,  as  CCl  began  to  mature  and  weaken,  CC3  began  to 
intensify.  Additionally,  the  simulation  reproduced  the  observed  tendency  for  intensi¬ 
fying  cells  to  entirely  absorb  the  condensate  and  circulation  of  surrounding  cells.  For 
example,  CC2  was  absorbed  into  CCl  as  CCl  and  CC3  merged.  The  simulation  to  a 
lesser  extent  captured  the  strong  observed  relationships  between  the  circulations  of 
neighboring  cells.  Only  in  a  few  cases  did  the  divergent  outflow  at  the  top  of  a  cell 
influence  or  feed  the  inflow  of  a  neighboring  cell.  The  best  example  of  this  occurred 
between  CC3  and  CCl  late  in  the  simulation. 
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6.2.2  Detailed  Comparisons 

Comparisons  of  simulation  results  at  selected  times  throughout  the  storm  develop¬ 
ment  with  radar  observations  and  synthesized  wind  fields  provide  more  detailed  in¬ 
sight  into  the  accuracy  of  the  simulation.  The  comparisons  were  made  using  the 
vertical  cross  sections  of  the  observed  storm  cell  referred  to  as  C2  in  Chapter  3.  This 
cell  was  chosen  because  of  its  proximity  to  the  location  of  the  simulation  domain. 
The  comparison  times  were  chosen  to  account  for  the  simulation's  apparent  lag  time 
from  the  observed  time  as  discussed  in  the  previous  section.  The  most  rapid  cell  de¬ 
velopment  occurred  in  the  simulation  immediately  after  penetration  of  the  capping  in¬ 
version.  Unfortunately,  the  corresponding  penetration  in  nature  occurred  during  the 
radar  data  gap  between  1119  and  1134  MST.  As  a  result,  no  comparisons  between 
the  simulation  and  observations  could  be  made  for  this  crucial  period  in  the  storm's 
development. 

The  first  comparison  was  made  using  the  simulation  results  at  6840  s  (Figs.  31 
and  32)  and  the  ob  •^rvations  from  1119  MST  (Fig.  8).  This  comparison  examines  the 
very  early  stages  of  cell  development.  Both  the  simulation  and  the  observations  show 
only  very  weak  development  occurring  at  this  time.  No  vertical  cross-sections  were 
made  from  the  observations  at  this  time,  however,  the  horizontal  cross-sections  show 
only  very  weak  echoes  existed  above  5.3  km  MSL  in  the  vicinity  of  C2.  The  simu¬ 
lated  cells  extended  to  only  3.75  km  (4.5  km  MSL).  The  observed  echo  and  the 
simulated  cell  were  both  approximately  1.5  km  wide.  The  observations  at  1119  MST 
did  not  support  the  synthesis  of  vertical  velocities  in  the  vicinity  of  the  echo.  The 
simulated  vertical  velocities  were  only  1.5  m  s  '. 

The  next  comparison  was  made  using  simulation  time  8460  s  and  observations 
from  1134  MST.  This  comparison  examines  the  early  part  of  the  cell's  mature  stage. 
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Figure  45  shows  the  simulated  total  condensate  and  wind  fields  for  8460  s,  and  Fig¬ 
ure  46  shows  the  radar  reflectivity  and  synthesized  wind  fields  for  1134  MST.  The 
two  figures  show  many  similar  features  present  in  both  the  simulation  and  in  the  ob¬ 
servations.  In  both  cases,  the  cell  extends  to  approximately  the  same  MSL  height. 
Both  cells  contain  a  relatively  strong  core  of  high  reflectivity  (observed)  or  total  con¬ 
densate  mixing  ratio  (simulated)  located  just  to  the  west  of  the  cell's  main  updraft. 
Both  cases  also  show  low-  to  mid-level  inflow  into  the  eastern  side  of  the  cell  with 
relatively  strong  mid-  to  upper-level  outflow  on  the  western  side.  Additionally,  both 
cases  exhibit  a  developing  weak  westerly  inflow  into  the  cell  at  very  low  levels. 

The  comparison  also  shows  several  marked  differences  between  the  simulation  and 
the  observations.  First,  the  simulated  cell  developed  several  kilometers  west  of  the 
observed  cell.  The  cause  of  this  difference  is  not  clear,  but  may  be  attributed  to  the 
lack  of  three-dimensional  forcing  effects  in  the  simulation  or  perhaps  an  unrealistical¬ 
ly  strong  initial  wind  profile.  Second,  the  simulation  contains  a  significant  total  con¬ 
densate  maximum  near  the  top  of  the  cell  which  does  not  appear  as  a  maximum  in  the 
observed  reflectivity  field.  Further  inspection  shows  the  simulated  maximum  con¬ 
sisted  primarily  of  pristine  ice  crystals.  The  lack  of  a  similar  upper-level  maximum  in 
the  observed  reflectivity  analysis  may  be  attributed  to  the  low  radar  reflectivity  of 
such  ice  crystals  (Battan,  1973).  This  phenomena  was  noted  in  almost  all  of  the  simu¬ 
lations  conducted  for  this  study.  Third,  the  simulated  cell's  updraft  was  significantly 
stronger  than  that  observed  in  C2  (14.6  m  s  '  vs.  9.0  m  s  ').  Fourth,  the  simulated 
cell  measured  only  about  3.5  km  across,  while  the  observed  cell  measured  nearly  7.0 
km  across. 

The  final  comparison  was  made  using  simulation  time  9000  s  and  observations 
from  1 152  MST.  This  comparison  examines  the  latter  part  of  the  cell's  mature  stage. 
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Fig.  45.  Analogous  to  Fig.  27,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  CONTROL  simulation  at 
A460  s.  Contour  label  "I "  represents  O.l  g  kg  '. 
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Fig.  46a.  Horizontal  cross-section  of  observed  radar  reflectivities  as  in  Fig.  lOa,  ex¬ 
cept  at  1134  MST  and  4.3  km  MSL.  Line  AB  shows  location  of  vertical  cross-section  in 
Fig.  46h.  Line  CD  was  not  used. 
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Fig.  46b.  As  in  Fig.  I  Oh,  vertical  cross-section  of  observed  radar  rtflectivities  and  vec¬ 
tor  velocities  at  1134  MST.  Cross-section  location  AB  is  shown  in  Fig.  46a. 
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Figure  47  shows  the  simulated  total  condensate  and  wind  fields  for  9000  s,  and  Fig¬ 
ure  48  shows  the  radar  reflectivity  and  synthesized  wind  fields  for  1152  MST.  It  is 
important  to  note  that  in  the  previous  two  comparisons,  the  single  simulated  cell 
(denoted  CCl  in  earlier  discussions)  was  compared  with  the  observed  cell  C2,  while 
in  this  comparison,  CCl  corresponds  more  closely  to  the  observed  cell  Cl.  While 
this  change  is  significant,  the  qualitative  nature  of  these  comparisons  still  allows  an 
assessment  of  the  simulation's  accuracy.  The  two  figures  again  show  many  similar 
features  between  the  simulation  and  the  observations.  In  both  cases,  two  cells  exist  in 
nearly  the  same  configuration.  The  western  cell  in  both  cases  extends  to  nearly  the 
same  MSL  height.  Both  western  cells  contain  strong  downdrafts  associated  with  the 
cells'  precipitation  shafts.  In  both  cases,  low-level  flow  enters  the  eastern  side  of  the 
eastern  cell,  rises  through  the  eastern  cell’s  updraft,  exits  the  western  side  of  the  east¬ 
ern  cell,  then  enters  the  western  cell’s  updraft  to  finally  exit  in  the  western  cell's  out¬ 
flow  region.  The  horizontal  dimensions  of  both  simulated  cells  compared  very 
favorably  with  those  observed.  Also  in  both  cases,  there  are  two  precipitation  shafts 
with  an  area  between  the  shaft  which  contained  little  or  no  precipitation. 

In  addition  to  the  change  in  cell  correspondence  described  in  the  previous  para¬ 
graph,  this  comparison  showed  several  other  major  differences  between  the  simulation 
and  the  observations.  Several  of  the  differences  noted  in  the  previous  comparison 
were  also  noted  in  this  comparison.  First,  the  simulated  cells  were  again  located  sev¬ 
eral  kilometers  west  of  the  observed  locations.  Second,  the  western  simulated  cell 
again  contained  a  tot^i  condensate  maximum  near  the  top  of  the  storm  with  no  cc^re- 
sponding  reflectivity  maximum  in  the  observations.  Third,  the  maximum  simulated 
vertical  velocity  exceeded  that  observed  in  CCl  by  over  4.0  m  s  ' .  Finally,  the  simu- 
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Fig.  47.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  as  in  Fig.  27,  except  for  CONTROL  simulation  at 
9000  s.  Contour  label  "I "  represents  O.I  g  kg  '. 
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Fig.  4Sa.  Analogous  to  Fig.  lOa,  horizontal  cross-section  of  observed  radar  reflectivi¬ 
ties  except  at  1152  MST  and  4.3  km  MSL.  Line  AB  shows  location  of  vertical  cross- 
section  in  Fig.  4Hh.  Line  CD  was  not  used. 
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Fig.  4Ht  Analogous  to  Fig.  I  Ob.  vertical  cross-section  of  observed  radar  reflectivities 
and  vector  velocities  at  1152  MST.  Cross-section  location  AB  is  .shown  in  Fig.  4Ha. 
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lated  eastern  cell  did  not  extend  nearly  as  high  as  the  observed  eastern  cell  (S.7  km 
[7.5  km  MSLJ  10.0  km  MSL). 

6.3  Comparison  With  Micro|rfiysical  Retrievals 

The  accuracy  of  the  simulated  storm's  microphysical  structure  and  morphology 
can  be  assessed  by  comparing  it  with  the  microphysical  structure  and  morphology  re¬ 
trieved  by  Lang  (1991)  from  the  Doppler  radar  observations  from  the  3  August  1984 
case.  This  section  presents  a  brief  overview  of  Lang's  (1991)  microphysical  retrieval 
results  and  a  comparison  of  those  results  with  the  simulated  microphysical  results. 

Lang's  (1991)  retrieval  showed  that  the  cell  development  could  be  divided  into 
two  distinct  stages:  the  first  phase,  in  which  the  cell  consisted  mostly  of  liquid  cloud 
and  rain  water;  and  the  second  stage  which  occurred  after  penetration  of  the  cell  into 
colder  air,  in  which  the  ice  phases  dominated  the  cell.  The  maximum  retrieved  cloud 
water  mixing  ratios  ranged  between  2.5  and  2.7  g  kg  '  throughout  the  analysis  period. 
The  cloud  water  existed  between  4.3  and  9.3  km  MSL.  The  maximum  retrieved  rain 
mixing  ratio  reached  5.0  g  kg  '  near  the  middle  of  the  first  phase  of  development. 
Lang  (1991)  reported  that  the  rain  mixing  ratio  maxima  in  this  stage  of  development 
existed  near  the  top  of  the  updrafts.  Later  maxima  which  occurred  in  the  second  stage 
of  development  occurred  much  lower  in  the  analysis  domain.  The  lower  elevation 
rain  maxima  in  the  second  stage  were  attributed  to  the  melting  of  hail/giaupel  falling 
from  higher  levels.  Rain  water  existed  between  3.3  and  7.8  km  MSL  through  most  of 
the  analysis  period  and  between  3.3  and  5.3  km  MSL  late  in  the  analysis.  The  maxi¬ 
mum  retrieved  cloud  ice  crystal  mixing  ratio  exceeded  3.0  g  kg '  about  midway 
through  the  second  phase  of  development.  Cloud  ice  crystals  existed  between  7.8  and 
11.8  km  MSL  with  the  maximum  located  near  10.0  km  MSL.  The  maximum  re¬ 
trieved  hail  mixing  ratio  reached  3.5  g  kg  '  early  in  the  second  phase  of  development. 
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The  hail  existed  between  4.3  and  8.3  km  MSL  during  the  entire  second  phase.  A  dis¬ 
tinct  double  maxima  existed  in  the  hail  mixing  ratio  analysis.  One  maximum  was  lo¬ 
cated  at  7.3  km  MSL  while  the  second  was  located  at  5.3  km  MSL.  Lang  (1991) 
noted  the  correspondence  of  these  maxima  with  the  top  of  the  updraft  and  with  the 
melting  level,  respectively.  Finally,  the  maximum  retrieved  snow  mixing  ratio 
reached  0.6  g  kg '  about  mid- way  through  the  second  phase.  The  snow  existed  be¬ 
tween  4.3  and  10.3  km  MSL  during  most  of  the  second  phase. 

In  most  respects,  the  simulated  microphysical  evolution  closely  resembled  that 
found  by  Lang's  (1991)  microphysical  retrievals  for  the  3  August  1984  case.  Ex¬ 
amination  of  the  simulated  microphysical  evolution  (Fig.  49)  showed  that  early  in  the 
cells'  development,  the  condensate  consisted  primarily  of  cloud  and  rain  water.  This 
agrees  well  with  Lang's  (1991)  observation  that  warm  rain  processes  dominated  the 
early  cell  development.  As  in  Lang's  (1991)  retrievals,  as  the  simulated  cells  devel- 
oj)ed  and  penetrated  into  the  cold  air,  the  ice  species  rapidly  began  to  dominate  the 
cells'  development. 

The  maximum  simulated  cloud  water  mixing  ratios  generally  ranged  between  2.5 
and  3.0  g  kg  '  with  a  few  maxima  near  4.0  g  kg  '.  The  cloud  water  was  located  be¬ 
tween  2.0  and  6. 1  km  (3.8  and  7.9  km  MSL).  These  values  and  locations  agree  quite 
well  with  those  found  in  the  retrieval.  The  maximum  simulated  rain  water  mixing 
ratio  reached  3. 1  g  kg  '  at  8730  s.  The  rain  water  existed  between  the  surface  and  4.0 
km  (5.8  km  MSL)  early  in  the  early  precipitating  stages  of  cell  development  and  be¬ 
tween  the  surface  and  3.2  km  (5.0  km  MSL)  in  the  later  stages  of  development.  The 
simulated  mixing  ratios  are  about  2.0  g  kg  '  less  than  those  retrieved.  Additionally, 
the  simulated  maximum  rain  water  mixing  ratio  occurred  much  later  in  the  develop¬ 
ment  than  in  the  retrieval.  In  fact,  in  the  simulation  the  rain  maximum  occurred  after 
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Fig.  49.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the  CON¬ 
TROL  .simulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash  cun’e 
shows  cloud  water  mixing  ratios,  dashed  curve  shows  rain  water  mixing  ratios,  double 
dot-dash  curve  shows  pristine  ice  crystal  mixing  ratios,  .shaded  cune  shows  graupel 
mixing  ratios,  and  dotted  curve  shows  aggregated  snow  flake  mixing  ratios.  Vertical 
axis  is  mixing  ratio  in  grams  per  kilogram.  Horizontal  axis  is  .simulation  time  in  .sec¬ 
onds  (bottom)  and  in  hours  and  seconds  MST  (top). 
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the  graupel  maximum,  while  in  the  retrieval  the  rain  maximum  occurred  before  the 
graupel  maximum.  The  simulation  did  show  a  tendency,  similar  to  that  seen  in  the 
retrieval,  for  the  rain  water  maximum  to  be  located  near  the  top  of  the  updraft  early 
in  the  development,  and  descend  to  lower  elevations  beneath  the  graupel/hail  shaft 
later  in  the  development. 

The  maximum  simulated  pristine  ice  crystal  mixing  ratio  reached  a  maximum  of 
6.8  g  kg '  with  a  typical  late  simulation  value  of  about  6.0  g  kg  '.  The  ice  crystals 
existed  between  3.7  and  11.2  km  (5.5  and  13.0  km  MSL).  The  simulated  mixing  ra¬ 
tios  were  about  twice  the  magnitude  of  those  found  in  the  retrieval.  The  ice  crystals' 
location  agreed  fairly  well  with  those  found  by  Lang  (1991).  The  maximum  simu¬ 
lated  graupel  mixing  ratios  reached  2.3  g  kg  '  at  8550  s  and  near  2.2  km  (4.0  km 
MSL).  The  graupel/hail  existed  between  the  surface  and  6.0  km  (7.8  km  MSL).  The 
simulated  mixing  ratios  agree  quite  well  with  those  retrieved.  The  simulation  fails  to 
depict  the  double  hail  maxima  exhibited  in  the  retrieval.  In  fact,  the  single  simulated 
graupel/hail  maximum  is  located  much  lower  than  the  retrieved  lower  elevation  maxi¬ 
mum.  Additionally,  the  simulated  graupel/hail  field  extends  to  the  surface  while  the 
retrieved  hail  shafts  extend  down  to  only  4.3  km  MSL.  The  maximum  simulated  ag¬ 
gregated  snowflake  mixing  ratio  reached  1.7  g  kg  '  very  late  in  the  simulation.  The 
aggregates  existed  between  3.2  and  10.7  km  (5.0  and  12.5  km  MSL).  The  maximum 
simulated  mixing  ratios  were  about  three  times  the  magnitude  of  those  in  the  retrieval. 
The  location  of  the  aggregates  agreed  quite  well  with  the  retrieval. 

The  differences  between  the  simulation  and  the  retrieval  may  be  attributed  to  sev¬ 
eral  factors.  First,  although  the  31  July  and  3  August  cases  were  similar,  there  were 
certainly  differences  between  them.  For  example  the  thermodynamic  and  wind  pro¬ 
files  were  certainly  different  between  the  two  days.  Second,  the  difference  in  the  mi- 
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crophysical  parameterizations  between  RAMS  and  the  microphysical  retrieval  model 
certainly  was  an  important  factor  in  the  differing  results. 


Chapter  7 

Sensitivity  To  Wind  Flow 

The  first  two  experiments  conducted  for  this  study  investigated  the  sensitivity  of 
the  simulated  storms  to  the  initial  base-state  wind  profiles.  In  Experiment  1  (EXPl), 
the  simulation  was  initialized  with  calm  winds  throughout  the  depth  of  the  atmo¬ 
sphere.  The  next  section  provides  a  brief  description  of  the  evolution  of  EXPl's 
simulated  storm.  Section  7,2  compares  EXPl’s  results  with  those  of  the  control.  In 
Experiment  2  (EXP2),  the  simulation  was  initialized  with  a  wind  profile  with  wind 
speeds  50%  of  those  in  the  control's  profile.  This  simulation  developed  serious 
boundary  condition  instabilities  which  rendered  the  simulation  useless  for  the  purposes 
of  this  study.  Section  7.3  provides  a  brief  description  of  the  problems  encountered 
with  this  simulation. 

7.1  Experiment  1  -  No  Wind  Case:  Analysis  Of  Simulated  Storm  Evolution 

Two  cells  developed  directly  over  the  mountain  peaks  early  in  the  simulation. 
The  cells  are  denoted  hereafter  as  ElCl  and  E1C2  for  the  western  and  eastern  cells, 
respectively.  Until  6300  s,  E1C2  dominated  with  maximum  total  condensate  mixing 
ratios  of  over  4.0  g  kg'^  and  maximum  vertical  velocities  of  1 1.0  m  s  '.  Precipitation 
first  developed  as  a  rain  shaft  out  of  E1C2  at  6300  s.  By  6300  s,  both  cells  already 
contained  some  pristine  ice  crystals.  Figures  SO  and  51  show  analyses  of  total  con¬ 
densate  mixing  ratio  and  vertical  velocity,  respectively,  for  6300  s. 

El  Cl  strengthened  slightly  and  E1C2  weakened  between  6300  and  6660  s.  Dur¬ 
ing  this  period,  rain  first  began  reaching  the  surface  from  ElCl  and  graupel  reached 
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Fig.  50.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors  for  EXPl  simulation  at  6300  s.  Contours  are  ev¬ 
ery  2.0  g  kg'  starting  at  0. 1  g  kg'.  Wind  vector  length  shown  in  lower  right-hand  cor¬ 
ner  of  figure  equals  50  m  s'.  Vertical  axis  is  labeled  in  kilometers  above  the  minimum 
terrain  elevation  in  the  simulation  grid  (in  this  case  1.82  km)  as  described  in  the  text. 
Horizontal  axis  is  horizontal  distance  labeled  in  kilometers.  The  cross-section  shown  is 
located  at  y  =  0.0  km  and  runs  east-west,  with  the  origin  at  Langmuir  Laboratory. 
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Fig.  51.  East-west  vertical  cross-section  of  simulated  vertical  velocity  contoured  every 
3.0  in  s'^  starting  at  1.5  m  s  '  for  EXPl  simulation  at  6300  s.  Solid  (dashed)  contours 
show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig.  50. 
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the  surface  tom  E1C2  at  6480  s.  Maximum  total  condensate  mixing  ratios  ranged 
between  3.0  and  3.7  g  kg  ',  while  maximum  updrafts  decreased  to  3.1  m  s  '  before 
increasing  again  in  E1C2  to  5.5  m  s  '.  The  cells  achieved  a  maximum  height  of  6.0 
km  during  this  period,  and  as  a  result,  the  ice  i^iase  dominated  the  development  as  the 
cells  pushed  into  the  colder  air  at  these  higher  levels. 

Between  6660  and  7560  s,  both  ElCl  and  E1C2  continued  to  develop.  The  cells 
exhibited  a  degree  of  the  oscillatory  development  individually  and  between  each  other 
as  noted  in  the  observations  described  in  Chapter  3  and  described  by  Ziegler  er  al. 
(1986).  During  this  period,  the  maximum  total  condensate  mixing  ratios  reached 
peaks  of  3.6  g  kg  '  in  E1C2  at  6840  s  and  4.2  g  kg  '  in  ElCl  at  7290  s.  The  maxi¬ 
mum  updrafts  during  this  period  reached  peaks  of  14.0  m  s  '  in  E1C2  at  6840s  and 
1 1.5  m  s  '  in  ElCl  at  7290  s.  During  most  of  the  period,  precipitation  only  fell  tom 
E1C2,  but  by  7380  s,  precipitation  was  falling  tom  both  cells.  Graupel  reached  the 
surface  again  tom  E1C2  at  7200  s.  SigniEcant  pristine  ice  crystal  development  also 
occurred  during  this  period.  Figure  52  shows  an  analysis  of  total  condensate  mixing 
ratio  for  7380  s. 

Between  7560  and  8190  s  further  development  occurred  in  both  ElCl  and  E1C2, 
and  a  new  cell  (hereafter  denoted  as  E1C3)  developed  just  to  the  east  of  E1C2.  The 
new  cell  developed  slowly  at  Erst,  but  by  8190  s  contained  total  condensate  mixing 
ratios  of  over  3.0  g  kg  '  and  updrafts  of  over  10.0  m  s  '.  PrecipitaEon  Erst  reached 
the  surface  tom  E1C3  between  7920  and  8100  s.  Meanwhile  ElCl  and  E1C2  devel¬ 
oped  to  a  maximum  height  of  8.0  km  and  contained  maximum  total  condensate  mix¬ 
ing  ratios  of  over  4.0  g  kg  '  and  maximum  updraEs  of  up  to  1 1.0  m  s  '.  During  this 
period,  ElCl  also  began  to  build  westward  slightly.  Precipitation  fell  throughout  the 
period  with  graupel  reaching  the  surface  tom  ElCl  at  7560  s,  tom  E1C2  at  7740  s. 
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Fig.  52.  Analogous  to  Fig.  50,  east-west  vertical  cross-sectUm  of  simulated  total  con- 
den.sate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP  I  .simulation  at 
73H0  s.  Contour  label  "I "  represents  O.I  g  kg  '. 
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and  from  E1C2  and  E1C3  between  7920  and  8100  s.  The  ice  phases  dominated  the 
development  during  this  time  with  pristine  ice  crystals  and  aggregated  snow  flakes 
showing  maximum  mixing  ratios  of  4.1  and  2.3  g  kg'‘,  respectively.  Figures  S3  and 
54  show  analyses  of  total  condensate  mixing  ratio  and  vertical  velocity,  respectively, 
for  8100  s. 

By  8280  s  E1C2  had  essentially  merged  with  ElCl.  At  the  same  time  a  new  cell 
(hereafter  denoted  as  E1C4)  developed  near  x  =  -4.0  km.  Between  8280  and  9(XX)  s, 
E1C4  rapidly  developed  from  a  small  cell  to  one  with  a  peak  maximum  total  conden¬ 
sate  mixing  ratio  in  excess  of  4.0  g  kg '  and  a  peak  maximum  updraft  in  excess  of 
15.0  m  s  '.  During  this  period,  E1C4  consisted  primarily  of  liquid  water  phases,  with 
ice  phases  developing  near  the  end  of  the  simulation.  The  cell  began  producing  rain 
at  the  surface  by  8820  s.  During  the  same  period,  El  Cl  continued  to  develop,  pri¬ 
marily  at  the  upper  levels,  where  the  ice  phases  were  dominant.  Maximum  total  con¬ 
densate  mixing  ratios  exceeded  4.0  g  kg  '  through  most  of  this  period  while  the 
maximum  updrafts  remained  between  7.5  and  10.5  m  s'.  The  cell  reached  a  maxi¬ 
mum  height  of  10.25  km  at  9(X)0  s  and  produced  rain  at  the  surface  between  8640  and 
9(XX)  s  and  graupel  between  8820  and  9(X)0  s.  E1C3  continued  its  rapid  development 
through  the  period  with  maximum  total  condensate  mixing  ratios  reaching  over  6.0 
g  kg  '  and  maximum  updrafts  in  excess  of  15.0  m  s  '  with  a  peak  of  19. 1  m  s  '  at 
8370  s.  This  cell  consisted  of  both  ice  and  the  liquid  water  phases.  The  cell  pro¬ 
duced  rain  at  the  surface  throughout  the  period  and  graupel  after  8460  s.  E1C3  devel¬ 
oped  to  reach  a  maximum  height  of  just  under  10.0  km  by  9(XX)  s.  Over  the  entire 
storm,  the  ice  phases  continued  to  dominate  during  this  period  with  pristine  ice  crystal 
mixing  ratios  reaching  a  peak  of  5.0  g  kg ',  aggregated  snow  flakes  a  peak  of  2.5 
g  kg ',  and  graupel  2.4  g  kg  '.  At  the  same  time,  though,  rain  and  cloud  water  mix- 
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Fig.  53.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors,  as  in  Fig.  50,  except  for  EXP  I  simulation  at  8100 
s.  Contour  label  "I "  represents  O.l  g  kg  '. 
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Fig.  54.  Analogous  to  Fig.  51.  east-west  vertical  cross-section  of  simulated  vertical 
locity  for  EXP  1  simulation  at  8100  s. 
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ing  ratios  reached  their  simulation  peaks  of  2.S  and  4.2  g  kg  respectively.  A  pro¬ 
nounced  low-level  inflow  developed  into  the  storm  over  this  time  period  also.  The 
low-level  flow  created  convergence  areas  near  the  location  where  E1C4  developed 
and  where  E1C3  continued  to  develt^.  This  low-level  convergence  continued  to  sup¬ 
port  the  development  of  these  cells  throughout  the  period.  Compensating  divergence 
near  the  top  of  the  storm  also  aided  continued  development.  Figures  55  and  56  show 
analyses  of  total  condensate  mixing  ratio  and  vertical  velocity,  respectively,  for  a  time 
representative  of  this  period-8820  s. 

Figures  57  and  58  show  the  evolution  of  the  simulation's  maximum  condensate 
mixing  ratios  and  upward  and  downward  vertical  velocities,  respectively. 

7.2  Experiment  1  -  No  Wind  Case:  Comparison  With  Control 

EXPl's  simulated  storm  showed  several  similarities  to  the  control  simulation's 
storm,  but  also  exhibited  many  differences.  Qualitatively,  the  storms  appear  similar. 
Both  storms  contain  several  cells- four  in  EXPl  compared  with  three  in  the  control. 
Both  attained  nearly  the  same  maximum  heights- 10. 25  km  in  EXPl  compared  with 
1  l.O  km  in  the  control.  Finally,  both  storms  generally  look  similar  in  shape. 

The  storms  exhibit  some  major  differences  also.  First,  EXPl's  storm  is  about  3.0 
km  wider  than  the  control.  Second,  EXPl's  cells  appear  to  have  had  a  more  equal 
chance  for  development  than  those  in  the  control.  The  two  major  cells  in  EXPl  grew 
to  about  the  same  height,  while  in  the  control,  CCl  dominated  the  vertical  develop¬ 
ment.  These  two  differences  may  be  attributed  to  the  absence  of  an  initial  environ¬ 
mental  wind  flow  in  EXPl.  Compared  to  the  control,  EXPl's  cells  likely 
experienced  less  entrainment  of  drier,  less  buoyant  environmental  air.  This  decreased 
entrainment  allowed  the  overall  storm  to  grow  more  in  horizontal  extent  than  the  con¬ 
trol,  and  also  allowed  EXPl's  eastern  cell  to  grow  stronger  than  the  eastern  cell  in  the 
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Fig.  55.  As  in  Fig.  50,  east-west  vertical  cross-section  of  simulated  total  conden.sate 
mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP  I  simulation  at  8?i20  s. 
Contour  label  "I "  represents  0. 1  g  kg  '. 
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Fig.  56.  A.S  in  Fig.  51,  east- west  vertical  cross-section  of  simulated  vertical  velocity  Ji. 
EXP  I  simulation  at  HH20  s. 
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Fig.  57.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the  EXP  I 
simulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash  curve  shows 
cloud  water  mixing  ratios,  da.shed  curve  shows  rain  water  mixing  ratios,  double  dot- 
dash  cwre  .shows  pristine  ice  crystal  mixing  ratios,  shaded  ctirs’e  .shows  graupel  mixing 
ratios,  and  dotted  curve  shows  aggregated  snow  flake  mixing  ratios.  Vertical  axis  is 
mixing  ratio  in  grams  per  kilogram.  Horizontal  axis  is  .simulation  time  in  seconds 
(bottom)  and  in  hours  and  seconds  MST  (top). 
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Fig.  58.  Time  evolution  of  maximum  and  minimum  simulated  vertical  velocities  for  the 
EXP  I  simulation.  Solid  curve  shows  maximum  upward  vertical  velocities  (positive  w). 
Dashed  curve  shows  maximum  downward  vertical  velocities  (negative  w}.  Vertical  axis 
is  vertical  velocity’  in  m  s' .  Horizontal  axis  is  .simulation  lime  in  seconds  (bottom)  and 
in  hours  and  seconds  MST  (top). 
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control.  Third,  the  convective  initiation  occurred  much  earlier  in  EXPl  than  in  the 
control"4230  s  in  EXPl  compared  to  5670  s  in  the  control.  The  initial  calm  winds 
limited  the  turbulent  transport  of  heat  away  from  the  surface,  allowing  maximum  so¬ 
lar  heating  and  rapid  convective  development.  Fourth,  the  cells  developed  between 
3.0  and  5.0  km  farther  east,  more  directly  over  the  mountain  tops,  in  EXPl  than  in 
the  control.  Several  studies  discussed  in  Chapter  2  (Banta,  1984;  Braham  and  Dragi- 
nis,  1960;  Glass  and  Carlson,  1963;  Orville,  1965b;  Barker  and  B^ta,  1985;  Banta 
and  Schaaf,  1987;  Schaaf  et  al.,  1988)  described  the  tendency  for  leeside  convective 
development  and  the  mechanisms  for  such  development.  In  the  absence  of  environ¬ 
mental  wind,  there  is  less  translation  of  convective  updrafts  toward  the  lee  of  the 
mountains  (Orville,  1965b).  Thus,  under  calm  wind  conditions  such  as  EXPl,  it  is 
natural  to  expect  the  storms  to  develop  closer  to  the  mountains  than  under  windier 
conditions  such  as  the  control.  Finally,  EXPl’s  storm  developed  more  steadily  and 
more  slowly  than  the  storm  in  the  control.  The  reason  for  this  difference  in  develop¬ 
ment  is  not  entirely  clear.  The  control’s  more  rapid  development  late  in  the  simu¬ 
lation  may  be  attributed  to  the  relatively  strong  flow  above  6.0  to  8.0  km  which 
served  as  an  "exhaust  mechanism"  for  the  developing  cell.  This  flow  helped  remove 
heat  from  the  top  of  the  cell  and  the  surrounding  environment,  keeping  the  environ¬ 
mental  lapse  rate  favorable  for  rapid  development.  The  absence  of  such  a  flow  in 
EXPl  may  have  contributed  to  its  slower,  more  steady  development. 

Quantitatively,  the  storms  are  markedly  different  in  many  respects.  The  maxi¬ 
mum  total  condensate  mixing  ratio  evolution  shown  in  Figure  57  shows  a  rapid  rise  in 
EXPl  between  5850  and  6210  s  followed  by  a  steady  increase  through  the  rest  of  the 
simulation.  This  development  is  much  different  than  the  extended  period  of  rapid  in¬ 
creases  in  the  maximum  total  condensate  mixing  ratio  seen  in  the  control.  Second, 
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the  peak  total  condensate  mixing  ratio  in  EXPl  reached  6.4  g  kg  '  compared  to  7.2 
g  kg  '  in  the  control.  EXPl's  pristine  ice  crystal  maximum  mixing  ratio  evolution 
differed  dramatically  from  that  seen  in  the  control.  In  EXPl,  the  maximum  mixing 
ratio  barely  reached  S.O  g  kg  ',  while  in  the  control,  the  maximum  pristine  ice  crystal 
mixing  ratio  never  went  below  S.O  g  kg  '  after  8460  s.  The  maximum  aggregated 
snow  flake  and  graupel  mixing  ratios  increase  earlier  aiid  persist  longer  in  EXPl  than 
in  the  control.  All  of  the  maximum  mixing  ratio  traces  in  Figure  57  show  a  large 
amount  of  variability  not  seen  in  the  control.  The  evolution  of  maximum  updrafts  and 
downdrafts  shown  in  Figure  58  show  a  similar  amount  of  variability.  The  maximum 
updrafts  exceeded  5.0  m  s  '  much  earlier  in  EXPl  than  in  the  control.  EXPl's  peak 
updraft  exceeded  the  control's  by  nearly  3.0  m  s  ',  while  EXPl's  peak  downdraft  ex¬ 
ceeded  the  control's  by  about  1.0  m  s'.  Again,  the  maximum  updrafts  and  down- 
drafts  indicate  a  much  slower,  steadier  development  in  EXPl  than  in  the  control. 

The  marked  differences  between  EXPl’s  simulation  and  the  control’s  precluded 
any  detailed  comparison  between  the  two  simulations  at  similar  analysis  times. 

7.3  Experiment  2  -  50%  Wind  Case 

EXP2  was  initialized  with  a  wind  profile  with  wind  speeds  50%  of  those  in  the 
control's  wind  profile.  This  simulation  developed  serious  boundary  condition  insta¬ 
bilities  which  ultimately  rendered  the  simulation  useless  for  the  purposes  of  this  study. 

The  simulation  started  very  much  like  the  control.  As  early  as  5400  s,  though,  the 
first  signs  of  the  boundary  problems  began  to  appear.  By  this  time  a  4.5  m  s  '  easter¬ 
ly  flow  developed  near  the  eastern  model  boundary  at  about  2.5  km.  A  westerly  flow 
of  3.0  m  s  '  had  developed  at  the  same  boundary  at  0.5  km.  This  circulation  was  ac¬ 
companied  by  a  relatively  large  negative  pressure  perturbation  (from  the  base  state)  in 
the  low  levels  near  the  same  boundary. 
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As  the  simulation  progressed,  the  pressure  perturbation  and  the  wind  maxima  con¬ 
tinued  to  increase  in  magnitude.  Additionally,  strong  updrafts  and  downdrafts  devel¬ 
oped  and  produced  extremely  large  areas  of  condensate.  The  spurious  circulations 
and  cells  continued  to  strengthen  throughout  the  simulation.  Eventually,  the  effect  of 
these  disturbances  destroyed  the  credibility  of  the  simulation  in  the  rest  of  the  model 
domain. 

The  cause  of  the  boundary  condition  problems  observed  in  this  simulation  are  un¬ 
clear.  The  problem  most  likely  resulted  from  some  reflection  or  excitation  of  a  cer¬ 
tain  wavelength  or  wavelengths  in  the  model.  It  is  not  obvious  why  the  50%  wind 
profile  produced  such  disturbances  while  the  stronger  winds  in  the  control  profile  did 
not.  The  problem  appears  to  be  related  to  the  topography  since  it  occurs  in  the  lower 
levels.  However,  the  topography  should  have  little  effect  since  the  initial  wind  pro¬ 
files  for  both  EXP2  and  the  control  had  calm  winds  below  the  mountain  top. 


Chapter  8 

Sensitivity  To  Diabatic  Heating  EffSects 

The  second  two  experiments  conducted  for  this  study  investigated  the  sensitivity  of 
the  simulated  storms  to  diabatic  heating  effects.  Section  8. 1  describes  the  results  of 
Experiment  3's  (EXP3)  simulation.  Section  8.2  provides  a  brief  description  of  the 
evolution  of  Experiment  4’s  (EXP4)  simulated  storm  and  a  comparison  of  the  results 
with  those  of  the  control. 

8.1  Experiment  3  -  No  Short  Wave  Radiation  Case 

In  EXP3,  the  simulation  was  conducted  with  the  RAMS  short-wave  radiation 
package  deactivated,  effectively  "turning  off  the  sun".  While  no  storm  developed  in 
this  case,  the  results  provide  crucial  evidence  of  the  importance  of  the  wind  flow  in 
focusing  the  convective  initiation. 

8.1.1  Analysis  Of  Simulated  Cloud  Evolution 

No  convective  development  occurred  during  the  simulation,  however,  a  small 
cloud  did  develop  and  persist  between  6300  s  and  9000  s  (Fig.  59).  The  cloud  con¬ 
sisted  entirely  of  cloud  water  and  the  maximum  mixing  ratio  never  exceeded  0.6 
g  kg  '.  A  very  weak,  but  steady,  updraft  of  between  2.5  and  3.4  m  s  '  supported  the 
small  cloud.  The  updraft  developed  near  a  quasi-stationary  area  of  low-level  conver¬ 
gence  located  between  x  =  -2.0  and  -3.0  km.  The  convergence  resulted  from  the 
interaction  of  the  initial  easterly  flow  with  a  weak  westerly  flow  below  2.0  km.  The 
westerly  flow  appears  to  have  developed  as  result  of  a  dynamically  induced  mountain- 
wave  circulation,  perhaps  even  a  rotor  circulation,  on  the  lee  side  of  the  mountain. 
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Fig.  59.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors  for  EXP3  simulation  at  8460  s.  ConUmrs  are  ev- 
ety  2.0  g  kg  '  starting  at  0. 1  g  kg  '.  Wind  vector  length  shown  in  lower  right-hand  cor¬ 
ner  nj  figure  equals  50  m  s''.  Vertical  axis  is  labeled  in  kilometers  above  the  minimum 
terrain  elevation  in  the  simulation  grid  (in  this  case  1.82  km)  as  described  in  the  text. 
Horizontal  axis  is  horizontal  distance  labeled  in  kilometers.  The  cross-section  shown  is 
located  at  y  =  0.0  km  and  runs  east-west,  with  the  origin  at  Langmuir  Laboratory. 
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Figure  60  shows  an  example  from  7740  s  of  the  convei^ence  described  above.  A 
similar  pattern  existed  well  into  the  simulation  to  8460  s. 

Both  the  cloud  and  the  updraft  developed  between  x  =  -2.0  and  -3.0  km,  and  re¬ 
mained  there  through  the  rest  of  the  simulation.  The  cloud  was  located  at  approxi¬ 
mately  z  =  2.75  km.  A  very  small  secondary  cloud  developed  just  to  the  east  of  the 
first  near  the  end  of  the  simulation. 

8.1.2  Comparison  With  Control 

Comparison  of  the  results  of  EXP3  with  those  of  the  control  shows  several  notable 
similarities.  In  both  cases,  the  first  clouds  formed  at  nearly  the  same  location.  The 
cloud  developed  slightly  higher  in  EXP3  than  in  the  control.  Also  in  both  cases,  the 
cloud  and  supporting  updraft  remained  nearly  stationary  for  an  extended  period.  In 
the  control,  the  stationary  cloud  and  updraft  remained  in  place  until  the  convective 
development  began  to  dominate  the  circulation  in  the  r  -ea.  Finally,  both  EXP3  and 
the  control  show  the  development  of  the  quasi-stationary  area  of  low-level  conver¬ 
gence  induced  by  the  mountain-wave  circulation  discussed  above.  This  comparison 
points  out  the  crucial  importance  of  environmental  wind  in  focusing  the  convective 
initiation. 

8.2  Experiment  4  -  Warm  Rain  Microirfiysics  Case 

EXP4  was  conducted  with  all  the  ice-phase  microphysical  processes  deactivated 
and  only  the  warm  rain  parameterizations  activated. 

8.2.1  Analysis  Of  Simulated  Storm  Evolution 

Prior  to  7920  s,  EXP4's  development  followed  exactly  or  nearly  exactly  that  de¬ 
scribed  in  Chapter  6  for  the  control  and  will  not  be  discussed  again  here.  By  7920  s, 
three  cells  had  developed,  E4C1  near  x  =  -5.0  km,  E4C2  slightly  to  the  east,  and 
E4C3  near  x  =  0.0  km.  E4C1  dominated  the  simulated  storm  with  total  condensate 
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Fig.  60.  East-west  vertical  cross-section  of  simulated  u  wind  component  contoured  e\'- 
ery  3.0  m  s'  starting  at  1.5  m  s' for  EXP3  simulation  at  7740  s.  Solid  (dashed)  con¬ 
tours  show  westerly  (easterly)  winds.  Axes  are  labeled  as  in  Fig.  59. 
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mixing  ratios  of  up  to  3.2  g  kg '  and  updrafts  of  10.0  m  s  '.  E4C2  contained  total 
condensate  mixing  ratios  in  excess  of  2.0  g  kg '  and  updrafts  in  excess  of  4.S  m  s  '. 
E4C3  contained  low  total  condensate  mixing  ratios  and  weak  updrafts  of  about  I..5 
m  s  '.  The  horizontal  flow  featured  an  area  of  low-level  convergence  near  x  =  -S.5 
km.  A  rain  shaft  associated  with  E4C1  and  E4C2  developed  between  y  =  -3.5  and 
-4.0  km.  The  peak  rain  mixing  ratio  in  the  storm  reached  0.4  g  kg  '.  The  maximum 
cloud  top  at  this  time  reached  5.25  km.  Figures  61  and  62  show  analyses  of  total 
condensate  mixing  ratio  and  vertical  velocity,  respectively,  for  7920  s. 

At  8100  s,  E4C1  continued  to  dominate  the  storm  with  maximum  total  condensate 
mixing  ratios  of  4.8  g  kg  '  and  maximum  updrafts  of  13.7  m  s  '.  E4C2  still  contained 
total  condensate  mixing  ratios  in  excess  of  2.0  g  kg  '  and  updrafts  of  over  7.5  m  s  '. 
A  significant  downdraft  of  5.1  ms'  developed  at  this  time  between  the  updrafts  asso¬ 
ciated  with  E4C1  and  E4C2.  E4C3  remained  weak  with  updrafts  near  1.5  m  s  '  and 
low  total  condensate  mixing  ratios.  The  rain  shaft  widened  to  about  1.0  km  with  a 
peak  storm  rain  mixing  ratio  of  1.6  g  kg  '.  The  maximum  cloud  top  reached  5.5  km 
by  this  time. 

By  8280  s,  both  E4C1  and  E4C2  had  weakened  slightly.  E4C1  continued  to 
dominate  with  total  condensate  mixing  ratios  of  4.8  g  kg '  and  a  maximum  updraft  of 
10.1  m  s  '.  E4Crs  updraft  showed  a  double  maxima  structure  with  the  upper  maxi¬ 
mum  located  near  z  =  4.5  km  and  the  lower  maximum  located  near  z  =  2.5  km. 
E4C2's  total  condensate  mixing  ratio  remained  slightly  above  2.0  g  kg  '  while  its  up¬ 
draft  weakened  to  about  4.5  m  s'.  Meanwhile  E4C3  grew  slightly  in  size,  but  re¬ 
mained  at  about  the  same  strength  as  at  previous  analysis  times.  The  rain  shaft 
associated  with  E4C1  and  E4C2  widened  again,  this  time  to  nearly  2.0  km.  The 
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Fig.  61.  Analogous  to  Fig.  59,  east-west  vertical  cross-section  of  simulated  total  con- 
densate  mixing  ratios  overlaid  with  .simulated  wind  vectors  for  EXP4  simidation  at 
7920  s. 
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Fig.  62.  East-west  vertical  cross-section  of  simulated  vertical  velocity  contoured  eveiy 
3.0  m  s'  starting  at  1.5  m  s'  for  EXP4  simulation  at  7920  s.  Solid  (dashed)  contours 
show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig.  59. 
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maximum  storm  rain  mixing  ratio  reached  1.9  g  kg '  in  E4C1,  and  the  maximum 
cloud  top  reached  6.25  km. 

At  8460  s  E4C1  clearly  dominated  the  storm.  The  cell  contained  a  double  maxi¬ 
ma  in  the  total  condensate  mixing  ratio.  The  upper  maximum  was  located  at  about 

5.5  km  and  contained  mixing  ratios  in  excess  of  4.0  g  kg  ^  The  lower  maximum  was 
located  near  4.0  km  and  contained  mixing  ratios  of  5.5  g  kg The  updraft  associated 
with  E4C1  increased  to  11.8  m  s  '.  Meanwhile,  E4C2  continued  to  weaken  with 
maximum  total  condensate  mixing  ratios  dropping  below  2.0  g  kg '  and  maximum 
updrafts  only  reaching  1.5  m  s  '.  At  the  same  time,  E4C3  developed  significantly 
over  previous  analyses  with  total  condensate  mixing  ratios  exceeding  2.0  g  kg '  and 
maximum  updrafts  exceeding  4.5  m  s  '.  The  maximum  rain  mixing  ratio  in  the  storm 
reached  3.5  g  kg  '  in  E4C1.  The  rain  shaft  widened  to  2.5  km,  and  the  maximum 
cloud  height  reached  6.5  km.  Figures  63  and  64  show  analyses  of  total  condensate 
mixing  ratio  and  vertical  velocity,  respectively,  for  8460  s. 

The  simulation  at  8640  s  showed  that  E4C1  continued  to  dominate  the  storm.  Its 
maximum  total  condensate  mixing  ratio  remained  nearly  steady  at  5.4  g  kg ',  and  its 
maximum  updraft  strengthened  to  13.0  m  s  '.  E4C2  remained  relatively  weak  with 
total  condensate  mixing  ratios  barely  exceeding  2.0  g  kg '  and  maximum  updrafts  of 

4.5  m  s  '.  E4C3  remained  at  about  the  same  strength  as  at  the  previous  analysis  time, 
but  grew  in  size.  At  this  time,  a  strong  mid-level  flow  out  of  E4C3,  through  E4C2, 
and  into  E4C1  matured.  The  rain  shaft  continued  to  widen  at  this  time  to  3.0  km. 
The  shaft  began  to  show  signs  of  a  split  between  E4C1  and  E4C2.  The  maximum 
rain  mixing  ratio  at  this  time  reached  4.4  g  kg '  in  E4C1,  while  the  maximum  cloud 
height  remained  at  6.5  km. 
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Fig.  63.  Analogous  to  Fig.  59,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP4  simulation  at 
8460  s.  Contour  label  "1 "  represents  0.1  g  kg  '. 
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Fig.  64.  Analogous  to  Fig.  62,  east-west  vertical  cross-section  of  simulated  vertical 
locity  for  EXP4  simulation  at  S460  s. 
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At  8820  s,  the  condensate  from  all  three  cells  joined  to  form  a  single  storm  domi¬ 
nated  by  E4C1.  The  cell's  updraft  contained  double  maxima  in  both  the  total  conden¬ 
sate  mixing  ratio  field  and  the  updraft.  The  maximum  total  condensate  mixing  ratio 
exceeded  4.0  g  kg  ',  and  the  maximum  updraft  reached  9.2  m  s  '.  E4C2  showed  little 
change  in  strength  from  the  previous  analysis  time.  A  strong  downdraft  of  5.7  m  s  ' 
developed  between  the  updrafts  associated  with  E4C1  and  E4C2.  E4C3  remained  at 
about  the  same  strength  as  in  the  previous  analysis.  The  rain  shaft  associated  with 
E4C1  and  E4C2  widened  to  3.5  km  with  a  split  between  E4C1  and  E4C2.  The  maxi¬ 
mum  rain  mixing  ratio  decreased  to  3.7  g  kg  *.  The  maximum  cloud  height  increased 
at  this  time  to  7.0  km. 

By  9000  s,  E4C1  still  dominated  the  storm  with  a  maximum  total  condensate  mix¬ 
ing  ratio  of  4.9  g  kg  '  and  a  maximum  updraft  of  10.2  m  s  '.  Both  of  these  maxima 
were  located  very  near  the  top  of  the  storm  at  elevations  above  6.0  km.  E4C2  be¬ 
came  nearly  absorbed  in  E4C1  and  E4C3  at  this  time.  E4C3  contained  two  updrafts 
with  the  strongest  exceeding  4.5  m  s  '.  The  maximum  total  condensate  mixing  ratio 
in  E4C3  exceeded  4.0  g  kg  '.  Rain  fell  out  of  all  three  cells  at  this  time,  with  a  maxi¬ 
mum  rain  mixing  ratio  of  3.1  g  1^  '.  The  maximum  cloud  height  reached  7.5  km. 
Figures  65  and  66  show  analyses  of  total  condensate  mixing  ratio  and  vertical  veloc¬ 
ity,  respectively,  for  9000  s. 

Figures  67  and  68  show  the  evolution  of  the  simulation's  maximum  condensate 
mixing  ratios  and  upward  and  downward  vertical  velocities,  respectively. 

8.2.2  Comparison  With  Control 

In  general,  the  results  of  EXP4's  simulation  were  similar  to  those  of  the  control- 
especially  given  the  marked  difference  in  the  microphysical  parameterizations  used. 
Prior  to  7920  s,  EXP4's  simulation  nearly  reproduced  exactly  that  seen  in  the  control. 
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Fig.  65.  Analogous  to  Fig.  59,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP4  .Emulation  at 
9000  s.  Contour  label  "1 "  represents  O.I  g  kg  '. 
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Fig.  66.  Analogous  to  Fig.  62,  east-west  vertical  cross-section  of  simulated  vertical 
locitv  for  EXP4  simulation  at  9000  s. 
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EXP4  CONDENSATE  MIXING  RATIOS 


Fig.  67.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the  EXP4 
.Emulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash  curve  shows 
cloud  water  mixing  ratios,  dashed  curve  shows  rain  water  mixing  ratios,  double  dot- 
dash  curve  shows  pristine  ice  crystal  mixing  ratios,  .shaded  ciore  shows  graupel  mixing 
ratios,  and  dotted  curve  shows  aggregated  snow  flake  mixing  ratios.  Vertical  axis  is 
mixing  ratio  in  grams  per  kilogram.  Horizontal  axis  is  .simulation  time  in  seconds 
(bottom)  and  in  hours  and  seconds  MST  (top). 
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EXP4  VERTICAL  VELOCITY 
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Fig.  68.  Time  evolution  of  maximum  and  minimum  simulated  vertical  velocities  for  the 
EXP4  simulation.  Solid  curve  shows  maximum  upward  vertical  velocities  (positive  w). 
Dashed  curve  shows  maximum  downward  vertical  velocities  (negative  w).  Vertical  axis 
is  vertical  velocity  in  m  s  '.  Horizontal  axis  is  simulation  time  in  seconds  (bottom}  and 
in  hours  and  seconds  MST  (top). 
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After  that  time,  the  two  simulations  began  to  differ  substantially,  but  they  still  gener¬ 
ally  corresponded  quite  well.  A  comparison  of  EXP4's  condensate  mixing  ratio 
evolution  (Fig.  67)  with  the  control's  evolution  (Fig.  49)  shows  nearly  perfect  agree¬ 
ment  in  the  timing  of  total  condensate  mixing  ratio  maxima  and  minima.  Naturally, 
the  values  do  not  agree  quantitatively,  but  the  degree  of  qualitative  agreement  is  re¬ 
markable.  A  similar  comparison  of  the  vertical  velocity  evolution  (Fig.  68  for  EXP4 
and  Fig.  43  for  the  control)  shows  good  agreement  prior  to  about  8280  s.  After  that 
time,  the  correspondence  between  the  two  simulations  deteriorates  with  the  control's 
maximum  updrafts  and  downdrafts  exceeding  the  strength  seen  in  EXP4. 

A  more  detailed  comparison  of  the  simulations  shows  similar  results.  Prior  to 
7920  s,  very  little,  if  any,  distinction  could  be  made  between  the  two  simulations. 
After  that  time,  several  of  EXP4's  storm  characteristics  were  weaker  or  less  signifi¬ 
cant  than  the  control's.  First,  EXP4's  maximum  cloud  heights  evolved  from  nearly 
the  same  height  as  the  control's  at  7920  s,  to  1.0  km  lower  than  the  control’s  at  8460 
s,  to  2.5  km  lower  at  8640  s,  to  3.5  km  lower  at  8820  s,  and  finally  to  4.5  km  lower 
than  the  control's  at  9000  s.  During  this  time  period,  the  ice  phases  dominated  the 
development  of  the  storm  in  the  control.  In  EXP4,  the  ice  phases  were  not  permitted 
and  thus  the  development  was  retarded  significantly.  Second,  the  maximum  updrafts 
in  EXP4  were  sometimes  5.0  to  7.0  m  s'  weaker  than  in  the  control.  Similarly, 
EXP4's  total  condensate  mixing  ratios  in  many  cases  were  up  to  2.0  g  kg  '  less  than 
in  the  control.  Again,  most  of  these  differences  can  be  attributed  to  the  lack  of  the 
ice  phases  in  the  simulation.  The  latent  heat  release  associated  with  the  ice  phases 
accounted  for  the  most  rapid  development  of  the  cells  in  the  control  during  this  time 
period.  Throughout  the  precipitating  phase  of  EXP4’s  simulation,  the  rain  shafts 
were  always  narrower  and  weaker  than  those  observed  in  the  control.  This,  too,  can 
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be  attributed  to  the  lack  of  the  ice  phases.  The  rain  shafts  in  the  control  were  very 
strongly  influenced  by  the  presence  of  graupel  shafts  aloft  above  the  rain  shafts. 
Without  the  graupel  in  EXP4,  the  rain  production  process  did  not  proceed  as  efficient¬ 
ly  as  it  did  in  the  control.  EXP4*s  simulated  cells  did  not  exhibit  the  strong  total  con¬ 
densate  mixing  ratio  maxima  found  near  the  top  of  the  control  simulation’s  cells.  The 
maxima  in  the  control  were  a  direct  result  of  maxima  in  the  ice  phases'  mixing  ratios, 
and  thus  were  not  a  factor  in  EXP4*s  simulation.  Finally,  both  E4C2  and  E4C3  were 
noticeably  less  significant  in  EXP4  than  their  counterparts  in  the  control  simulation. 


Chapter  9 

Senativity  To  Initial  Moisture  Profile 

The  final  experiment  conducted  for  this  study  investigated  the  sensitivity  of  the 
simulated  storm  to  the  initial  base-state  moisture  profile.  Section  9. 1  provides  a  brief 
description  of  the  modified  moisture  profile  used  to  initialize  EXP5,  Section  9.2 
briefly  describes  the  evolution  of  EXPS's  simulated  storm,  and  Section  9.3  compares 
EXP5's  results  with  those  of  the  control. 

9.1  Experiment  5  -  Moistened  Lower  Levels  Case 

EXP5  was  designed  to  test  the  sensitivity  of  the  simulated  storm  development  to 
the  initial  atmospheric  moisture  profile.  As  briefly  described  in  Chapter  5,  EXPS's 
initial  moisture  profile  was  constructed  to  produce  a  profile  with  a  higher  moisture 
content  in  the  lower  atmosphere.  For  this  case  the  initial  mixing  ratio  was  increased 
in  the  lowest  model  levels  by  0.5  g  kg  '  over  that  in  the  control's  base-state  moisture 
profile.  The  amount  of  increase  over  the  control  was  gradually  reduced  to  zero  at  720 
mb.  Above  720  mb,  EXPS's  initial  moisture  profile  and  the  control's  were  identical. 
Figure  69  shows  EXPS's  moisture  profile  superimposed  over  the  control's. 

9.2  Analysis  Of  Simulated  Storm  Evohitiou 

By  6840  s,  two  cells,  ESCl  near  x  —  -3.0  km  and  E5C2  near  x  =  -1.0  km,  had 
developed.  ESCl  contained  a  maximum  updraft  of  5.4  m  s  '  and  a  maximum  total 
condensate  mixing  ratio  of  2.2  g  kg  '.  The  cells'  development  prior  to  this  time  fol¬ 
lowed  that  of  the  control  with  EXPS's  development  occurring  slightly  more  rapidly 
than  the  control's.  For  example,  by  6840  s,  the  control's  maximum  updraft  and  total 
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Fig.  69.  Skew  T,  log  P  plot  of  modified  base-state  moisture  profile  used  to  initialize  the 
numerical  model  for  the  EXP5  simulation.  Heavy  solid  and  dashed  lines  represent  sen¬ 
sible  and  dewpoint  temperature  profiles,  respectively.  Light  dashed  line  shows  original 
base-state  moisture  profile  for  reference.  Skew  T.  log  P  diagram  labeled  as  in  Fig.  2. 


181 


condensate  mixing  ratio  reached  only  3.2  m  s  '  and  1.2  g  kg  ',  respectively,  compared 
with  EXPS’s  5.4  m  s  '  and  2.2  g  kg  '.  Figures  70  and  71  show  analyses  of  total  con¬ 
densate  mixing  ratio  and  vertical  velocity,  respectively,  for  6840  s. 

Between  6840  and  7380  s,  the  storm  matured  significantly  with  E5C1  dominating 
the  storm  until  7200  s  and  E5C2  dominating  late  in  the  period.  The  cells  exhibited  a 
degree  of  the  oscillatory  development  individually  and  between  each  other  as  noted  in 
the  observations  described  in  Chapter  3  and  described  by  Ziegler  et  ul.  (1986).  E5C1 
dominated  the  storm  until  7200  s  with  maximum  updrafts  reaching  6.2  m  s  '  at  7020  s 
and  maximum  total  condensate  mixing  ratios  reaching  2.6  g  kg  '  at  the  same  time. 
Tliese  values  diminished  quickly  to  about  1.5  m  s  '  and  less  tlian  2.0  g  kg  '  by  7380  s. 
A  new  updraft  of  about  1.5  m  s  '  developed  in  association  with  E5Ci  after  7200  s 
near  .v  =  -5.0  km.  By  7380  s,  cloud  water,  pristine  ice  crystals,  rain,  and  graupel 
had  developed  in  E5CI,  As  E5C1  was  weakening  late  in  the  period,  E5C2  strength¬ 
ened  significantly  to  contain  peak  updrafts  of  5.5  m  s  '  and  total  condensate  mixing 
ratios  of  2.4  g  kg '.  At  7380  s,  a  new  cell  (denoted  hereafter  as  E5C3)  developed 
near  .v  =  1.0  km.  Figures  72  and  73  show  analyses  of  total  condensate  mixing  ratio 
and  vertical  velocity,  respectively,  for  7380  s.  During  this  period,  two  areas  of  low- 
level  convergence  developed:  one  near  x  =  -3.0  km  and  the  other  near  x  =  -5.0  km 
(Fig.  74).  The  convergence  resulted  from  the  interaction  of  wind  flow  patterns  simi¬ 
lar  to  those  described  in  the  analysis  of  the  control  simulation's  evolution. 

E5C3  had  developed  quite  rapidly  by  7560  s.  The  cell  contained  maximum  up¬ 
drafts  of  10.4  m  s  '  and  a  peak  total  condensate  mixing  ratio  of  over  2.0  g  kg  '.  By 
this  time,  cloud  water,  pristine  ice  crystals,  graupel,  and  rain  water  had  developed  in 
the  cell.  E5C2  remained  relatively  unchanged  from  the  previous  analysis  time  with 
maximum  updrafts  of  over  4.5  m  s  '  and  maximum  total  condensate  mixing  ratios  of 
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Fig.  70.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  simulated  wind  vectors  for  EXP5  simulation  at  6840  s.  Contours  are  ev¬ 
ery  2.0  g  kg  '  .starting  at  0. 1  g  kg  '.  Wind  vector  length  shown  in  lower  right-hand  cor¬ 
ner  of  figure  equals  50  m  Vertical  axis  is  labeled  in  kilometers  above  the  minimum 
terrain  elevation  in  the  simulation  grid  (in  this  case  1.82  km)  as  described  in  the  text. 
Horizontal  axis  is  horizontal  distance  labeled  in  kilometers.  The  cross-section  .shown  is 
located  at y  =  0.0  km  and  runs  east-west,  with  the  origin  at  Langmuir  Laboratory. 
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Fig.  71.  East-west  vertical  cross-section  of  simulated  vertical  velocity  contoured  every 
3.0  m  s  '  starting  at  1.5  m  s  '  for  EXP 5  simulation  at  6840  s.  Solid  (dashed)  contours 
show  positive  (negative)  vertical  velocities.  Axes  labeled  as  in  Fig.  70. 
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Fli,.  72.  Analogous  to  Fig.  70,  east-west  vertical  cross-section  of  simulated  total  con- 
cten.mte  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP5  simulation  at 
73H0  s.  Contour  label  "1 "  represents  0.1  g  kg  '. 
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Fig.  73.  Analogous  to  Fig.  71,  east-west  vertical  cross-section  of  simulated  vertical  ve¬ 
locity  for  EXP5  .simulation  at  7380  s. 
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Fig.  74.  East-west  vertical  cross-section  of  simulated  u  wind  component  contoured  ev¬ 
ery  3.0  m  s  '  starting  at  1.5  m  s  ' for  EXP5  simulation  at  7380  s.  Solid  (dashed)  con¬ 
tours  show  westerly  (easterly)  winds.  Axes  are  labeled  as  in  Fig.  70. 
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3.0  g  kg  '.  E5CI  remained  weak  with  only  very  weak  updrafts  and  maximum  UMal 
condensate  mixing  ratios  of  less  than  2.0  g  kg  '. 

By  7740  s,  E5C2  clearly  dominated  the  storm  with  maximum  updrafts  of  11.8 
m  s  '  and  peak  total  condensate  mixing  ratios  of  4.3  g  kg  '.  Graupel  developed  for 
the  first  time  in  E5C2  at  this  time.  E5C2  also  contained  a  relatively  strong  downdraft 
of  4.5  m  s  '.  At  the  same  time,  ESC3  weakened  slightly  with  a  maximum  updraft  of 
7.5  m  s  '  and  maximum  total  condensate  mixing  ratios  of  over  2.0  g  kg  '.  This  cell 
also  contained  a  downdraft  in  excess  of  4.5  m  s  '.  E5C1  strengthened  slightly  to  in¬ 
clude  maximum  total  condensate  mixing  ratios  of  over  2.0  g  kg ' . 

At  7920  s,  E5C3  dominated  the  storm  with  a  maximum  total  condensate  mixing 
ratio  of  4.4  g  kg '.  The  cell's  updraft  contained  a  double  maxima  with  the  upper 
maximum  (located  near  4.0  km)  reaching  1 1.2  m  s  '  and  the  lower  maximum  (located 
near  3.0  km)  reaching  over  4.5  m  s  '.  The  cell  also  contained  a  downdraft  of  4.7 
m  s  '.  The  first  precipitation  reached  the  surface  at  this  time  as  a  rain  shaft  from 
E5C3.  E5C2  weakened  from  the  previous  analysis  time  as  its  updraft  split  into  two 
maxima.  The  upper  maximum,  located  near  4.0  km,  exceeded  7.5  m  s  ',  while  the 
lower  maximum,  located  near  3.0  km,  exceeded  4.5  m  s  '.  The  cell's  maximum  total 
condensate  still  exceeded  2.0  g  kg  '.  E5C1  remained  very  weak  with  updrafts  still 
near  1.5  m  s  '  and  maximum  total  condensate  mixing  ratios  dn^ping  below  2.0 
g  kg  '.  Figures  75  and  76  show  analyses  of  total  condensate  mixing  ratio  and  vertical 
velocity,  respectively,  for  7920  s. 

By  8100  s,  E5C1  showed  signs  of  leintensification  while  E5C3  weakened  dramati¬ 
cally.  E5C3  contained  maximum  updrafts  of  over  4.5  m  s  '  and  maximum  total  con¬ 
densate  mixing  ratios  of  3.9  g  kg  '.  Rain  continued  to  fall  out  of  this  cell.  E5C2 
maintained  its  earlier  strength  with  a  maximum  updraft  in  excess  of  4.5  m  s  '  and 
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Fig.  75.  East-west  vertical  cross-section  of  simulated  total  condensate  mixing  ratios 
overlaid  with  .simulated  wind  vectors,  as  in  Fig.  70,  except  for  EXP5  simulation  at  7920 
s.  Contour  label  "/  "  represents  O.l  g  kg  '. 
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maximum  total  condensate  mixing  ratios  of  over  2.0  g  kg  ‘.  The  first  precipitation 
developed  out  of  E5C2  at  this  time  as  a  rain  shaft.  Meanwhile,  ESCl’s  maximum 
updraft  increased  to  7.5  m  s  '  while  its  maximum  total  condensate  mixing  ratio  re¬ 
mained  below  2.0  g  kg  ' .  The  maximum  cloud  top  reached  6.0  km  at  this  time.  The 
ice  phases  began  to  dominate  at  this  time  as  ESC2  pushed  into  the  colder  air  at  these 
higher  levels.  ESCl's  intensification  at  this  time  resulted  from  the  intensification  of 
the  low-level  convergence  near  x  =  -5.5  km.  This  intensification  was  accompanied 
by  a  gradual  weakening  of  the  convergent  area  nearx  =  -3.0  km  (Fig.  77). 

E5C1  continued  to  develop  rapidly  at  8280  s  with  maximum  updrafts  of  12.2  m  s  ' 
and  maximum  total  condensate  mixing  ratios  exceeding  2.0  g  kg '.  The  first  preci¬ 
pitation  from  E5C1  developed  as  a  rain  shaft  at  this  time.  Also,  the  first  graupel  de¬ 
veloped  in  the  cell  at  this  time.  All  of  E5Crs  condensate  became  very  closely 
merged  with  that  of  E5C2  at  this  time.  E5C2's  strength  remained  nearly  steady  with 
peak  updrafts  in  excess  of  4.5  m  s  '  and  a  maximum  total  condensate  mixing  ratio  of 
4.4  g  kg '.  The  condensate  maximum  occurred  near  the  top  of  the  cell  and  consisted 
primarily  of  pristine  ice  crystals.  The  first  aggregated  snow  flakes  of  the  simulation 
developed  near  this  maximum  at  this  time.  Rain  continued  to  fall  from  E5C3,  but  its 
maximum  updraft  decreased  to  just  over  1.5  m  s  '  while  its  maximum  total  condensate 
mixing  ratio  dropped  below  2.0  g  kg  '.  At  the  same  time,  graupel  from  the  cell  near¬ 
ly  reached  the  surface  and  some  weak  development  of  aggregated  snow  flakes  oc¬ 
curred  near  the  top  of  the  cell. 

By  8460  s,  E5Crs  updraft  dominated  the  storm  with  a  maximum  value  of  14.8 
m  s  '.  The  cell's  maximum  total  condensate  mixing  ratio  continued  to  remain  above 
2.0  g  kg'.  E5C2  contained  a  double  maxima  updraft  with  the  upper  maximum  ex¬ 
ceeding  7.5  m  s  '  and  the  lower  maximum  exceeding  4.5  m  s  '.  The  cell's  maximum 
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Fig.  77.  East-west  vertical  cross-section  of  simulated  u  wind  component  as  in  Fig.  74 
except  for  EXP5  simulation  at  8100  s. 
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total  condensate  mixing  ratio  reached  S.3  g  kg  ‘.  At  this  time,  the  cell  contained 
some  rain  water  which  did  not  reach  the  surface.  A  relatively  strong  downdraft  of 
4.5  m  s  '  developed  between  the  updrafts  associated  with  E5C1  and  E5C2.  E5C3 
continued  to  dissipate  with  only  very  weak  updrafts  aloft  and  a  small  amount  of 
condensate— mainly  pristine  ice  crystals.  A  new  cell  developed  at  this  time  (denoted 
hereafter  as  E5C4)  near  x  =  2.5  km.  The  new  cell  contained  a  maximum  updraft  of 
over  4.5  m  s  ',  but  maximum  total  condensate  mixing  ratios  did  not  exceed  2.0  g  kg  ' . 

E5C1  clearly  dominated  the  storm  by  8640  s  with  a  maximum  updraft  of  16.3 
m  s  '  and  peak  total  condensate  mixing  ratios  of  4.7  g  kg  '.  E5C1  developed  its  own 
distinct  rain  shaft  at  this  time.  E5C2  continued  to  exhibit  a  double  maxima  updraft 
with  the  lower  maximum  reaching  over  7.5  m  s  '  while  the  upper  maximum  exceeded 
4.5  m  s  '.  The  cell's  maximum  total  condensate  mixing  ratio  exceeded  4.0  g  kg '. 
For  the  purposes  of  this  discussion  E5C3  had  dissipated  by  this  time.  E5C4  remained 
at  about  the  same  strength  as  at  the  previous  analysis  time  with  maximum  updrafts  of 
over  4.5  m  s  '  and  maximum  total  condensate  mixing  ratios  of  less  than  2.0  g  kg  '. 
Rain  water  developed  in  the  cell  at  this  time,  but  did  not  reach  the  surface.  Figures 
78  and  79  show  analyses  of  total  condensate  mixing  ratio  and  vertical  velocity,  re¬ 
spectively,  for  8640  s.  The  horizontal  wind  analysis  showed  an  area  of  strong  upper- 
level  divergence  over  the  area  of  low-level  convergence  associated  with  E5C1 . 

E5C1  continued  to  dominate  the  simulation  at  8820  s  with  a  large  updraft  contain¬ 
ing  several  maxima.  The  absolute  maximum  had  decreased  slightly  to  12.6  m  s  '. 
The  cell  contained  a  double  maxima  in  the  total  condensate  mixing  ratio  with  the  up¬ 
per  maximum  reaching  6.8  g  kg  '  and  die  lower  maximum  reaching  over  4.0  g  kg  '. 
Aggregated  snow  flakes  developed  and  the  graupel  shaft  reached  the  surface  for  the 
first  time  in  E5C1  at  this  dme.  E5C2  remained  at  about  the  same  strength  as  the  pre- 
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Fi^.  78.  Analogous  to  Fig.  70,  east-west  vertical  cross-section  of  simulated  total  con¬ 
densate  mixing  ratios  overlaid  with  simulated  wind  vectors  for  EXP5  simulation  at 
8640  s.  Contour  label  "I "  represents  O.J  g  kg  '. 
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Fig.  79.  Analogous  to  Fig.  71.  east-west  vertical  cross-section  of  simulated  vertical  ve¬ 
locity  for  EXP5  simulation  at  8640  s. 
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vious  analysis  time  with  peak  updrafts  of  over  7.5  m  s  ’.  The  cell  contained  several 
total  condensate  mixing  ratio  maxima  with  most  of  them  over  4.0  g  kg  E5C4  con¬ 
tinued  to  intensify  with  peak  updrafts  reaching  over  7.5  m  s  '  and  maximum  total  con¬ 
densate  mixing  ratios  of  over  2.0  g  kg  ’.  Pristine  ice  crystals  and  graupel  developed 
in  E5C4  at  this  time. 

By  the  end  of  the  simulation  at  9000  s,  E5C1  continued  to  dominate  the  storm. 
The  cell  contained  several  updraft  maxima  with  the  strongest  peaking  at  14.6  m  s  '. 
The  maximum  total  condensate  mixing  ratio  reached  7. 1  g  kg  '  near  the  top  of  the 
cell.  The  aggregated  snow  flake  mixing  ratio  increased  to  over  2.0  g  kg  '  for  the  first 
time.  E5C2  weakened  slightly  with  a  maximum  updraft  continuing  to  exceed  7.5 
m  s  ',  but  with  the  maximum  total  condensate  mixing  ratio  dropping  to  just  over  2.0 
g  kg  '.  E5C4  remained  at  about  the  same  strength  as  before  with  updrafts  in  excess  of 
4.5  m  s  '  and  maximum  total  condensate  mixing  latios  of  over  2.0  g  kg  '. 

Figures  80  and  81  show  the  evolution  of  the  simulation’s  maximum  condensate 
mixing  ratios  and  upward  and  downward  vertical  velocities,  respectively. 

9.3  Comparison  With  Control 

Given  the  relatively  small  change  in  the  initial  moisture  profile,  the  results  of 
EXP5's  simulation  differed  dramatically  from  those  of  the  control.  EXP5’s  simu¬ 
lation  produced  four  cells  compared  to  three  cells  in  the  control.  The  first  three  cells 
to  develop  in  EXP5  developed  in  generally  the  same  locations  as  their  counterparts  in 
the  control.  The  development  of  these  cells  differed  dramatically  from  those  in  the 
control.  While  the  western-most  cell  (E5C1  in  EXP5  and  CCl  in  the  control)  ulti¬ 
mately  dominated  the  simulation  in  both  cases,  in  EXP5,  this  domination  was  preced¬ 
ed  by  strong  development  of  E5C2  and  E5C3.  In  fact,  the  rapid  development  of 
E5C3  seemed  to  retard  for  a  time  the  development  of  the  two  cells  to  its  west.  The 
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EXP5  CONDENSATE  MIXING  RATIOS 


SIMULATION  TIME  (SEC) 


Fig.  80.  Time  evolution  of  maximum  simulated  condensate  mixing  ratios  for  the  EXP5 
simulation.  Solid  curve  shows  total  condensate  mixing  ratios,  dot-dash  curve  shows 
cloud  water  mixing  ratios,  dashed  curx’e  shows  rain  water  mixing  ratios,  double  dot- 
dash  curve  shows  pristine  ice  crystal  mixing  ratios,  shaded  cuire  .shows  graupel  mixing 
ratios,  and  dotted  curve  shows  aggregated  snow  flake  mixing  ratios.  Vertical  axis  is 
mixing  ratio  in  grams  per  kilogram.  Horizontal  axis  is  simulation  time  in  .seconds 
(bottom}  and  in  hours  and. seconds  MST (top). 
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Fig.  81.  Time  evolution  of  maximum  and  minimum  simulated  vertical  velocities  for  the 
EXP5  simulation.  Solid  curve  shows  maximum  upward  vertical  velocities  (positive  w). 
Dashed  cur\’e  shows  maximum  downward  vertical  velocities  (negative  w).  Vertical  axis 
is  vertical  velocity  in  m  s  '.  Horizontal  axis  is  .simulation  time  in  seconds  ( bottom t  and 
in  hours  and  seconds  MST  (top). 
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circulation  and  cool  downdrafts  associated  with  E5C3  inhibited  the  flow  of  warm, 
moist  air  into  the  cells  to  its  west,  thus  slowing  their  development.  E5C2  maintained 
its  strength  much  longer  into  the  simulation  than  its  counterpart  cell  in  the  control. 

A  comparison  of  EXP5's  maximum  condensate  ratio  evolution  (Fig.  80)  with  the 
control's  evolution  (Fig.  49)  shows  several  differences.  First,  EXPS's  development 
appears  steadier,  and  did  not  reach  its  absolute  maximum  total  condensate  mixing  ra¬ 
tio  until  the  end  of  the  simulation.  The  early  development  (before  8280  s)  produced 
higher  total  condensate  mixing  ratios  in  EXP5  than  in  the  control.  Second,  EXPS's 
rain  water  mixing  ratios  were  slightly  higher  early  in  the  simulation  and  then  lagged 
those  of  the  control  by  up  to  1.0  g  kg  '  later  in  the  simulation.  EXPS's  peak  graupel 
mixing  ratios  were  higher  much  earlier  than  in  the  control.  The  pristine  ice  crystal 
mixing  ratios  in  both  cases  followed  about  the  same  development  pattern  with  a  rapid 
increases  between  7920  s  and  8280  s. 

A  similar  comparison  of  the  vertical  velocity  evoludon  (Fig.  81  for  EXPS  and 
Fig.  43  for  the  control)  shows  general  agreement  in  the  upward  trend  of  updraft 
strength.  The  comparison  shows  little  correspondence  between  the  timing  of  the 
maxima  and  minima  in  the  updrafts  of  both  simulations.  The  absolute  maximum  for 
EXPS  was  about  4.0  m  s  '  greater  than  in  the  control.  A  comparison  of  the  downdraft 
time  series  shows  generally  good  agreement  between  both  simulations. 

A  more  detailed  comparison  of  the  simulations  shows  similar  results.  Prior  to 
7S60  s,  the  simulations  in  both  EXPS  and  in  the  control  were  very  similar.  The  dif¬ 
ferences  were  minor  with  EXPS's  updrafts  slightly  stronger  than  the  control's  and  the 
maximum  total  condensate  mixing  ratio  in  EXPS  reaching  peaks  about  O.S  g  kg  ' 
greater  than  the  control's  between  6840  and  7S60  s.  Between  7S60  and  7920  s,  ESC3 
developed  very  rapidly  and  even  began  producing  precipitation.  The  corresponding 
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cell  in  the  control,  CC3,  showed  only  very  weak  development  during  this  period. 
EXPS's  two  western  cells,  E5C1  and  E5C2,  still  were  not  producing  precipitation  by 
7920  s,  while  the  merged  storm  in  the  control  was  producing  precipitation  by  this 
time.  Between  8100  and  8460  s,  as  ESC3  diminished,  the  pattern  of  the  merged 
storm  in  EXPS  appeared  very  similar  to  that  in  the  control.  The  major  difference 
during  this  period  was  that  in  EXP5,  the  dominant  ceil  was  E5C2,  the  eastern  cell, 
while  in  the  control,  the  dominant  cell  was  CCl,  the  western  cell.  The  updrafts  and 
maximum  total  condensate  mixing  ratios  for  the  two  simulations  agreed  reasonably 
well  during  this  period.  During  the  period  from  8640  to  9000  s,  the  maximum  cloud 
tops  in  EXPS  were  consistently  lower  than  those  in  the  control.  The  difference 
ranged  from  0.75  km  at  8640  s  to  1.75  km  at  8820  s.  Also  during  this  time,  both 
E5C1  and  E5C2  attained  approximately  the  same  height  while  in  the  control,  CCl 
grew  to  much  higher  levels  than  did  CC2. 
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Chapter  10 
Conchisioiis 

This  research  examined  the  initiation,  development,  and  maintenance  of  a  small 
mountain  thunderstorm  using  numerical  simulations.  The  storm  considered  in  this 
study  occurred  on  31  July  19S4  and  is  typical  of  the  many  storms  which  form  each 
summer  in  and  around  the  Magdalena  Mountains  of  central  New  Mexico.  As  a  result 
of  their  frequency  and  relatively  fixed  location,  the  storms  have  been  studied  exten¬ 
sively  for  at  least  the  past  15  years.  The  current  study  examines  the  31  July  1984 
storm  using  the  two-dimensional,  non-hydrostatic  cloud  model  configuration  of  the 
Colorado  State  University  Regional  Atmospheric  Modeling  System  (Tripoli  and  Cot¬ 
ton,  1982,  1986;  Cotton  et  ai,  1982,  1986;  Tripoli,  1986).  The  study  also  employs 
the  multiple  Doppler  analyses  performed  by  Lang  (personal  communication,  1992) 
following  the  technique  of  Ray  etal.  (1981)  and  described  in  detail  by  Lang  (1991)  as 
well  as  the  three-dimensional  microphysical  retrieval  performed  by  Lang  (1991)  for  a 
similar  storm  that  occurred  on  3  August  1984. 

The  numerical  model  was  used  to  isolate  the  processes  and  initial  conditions  most 
important  in  the  initiation,  development,  and  maintenance  of  these  types  of  storms.  A 
total  of  six  simulations  were  conducted— one  control  and  five  experiments.  The  re¬ 
sults  of  the  control  siniaLrr.on  were  compared  with  Doppler  analyses  and  with  Lang's 
(1991)  retrieval  to  assess  the  validity  of  the  simulation  and  its  usefulness  in  studying 
the  storm.  The  experiments  were  conducted  to  test  the  sensitivity  of  the  simulated 
storms  to  initial  wind  profile,  the  lack  of  solar  heating,  restriction  to  warm  rain 
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processes,  and  the  initial  moisture  profile.  The  results  of  these  simulations  were  com¬ 
pared  with  the  control's  results. 

Analyses  of  the  Doppler  radar  observations  and  synthesized  wind  fields  provided 
the  basis  for  evaluating  the  accuracy  of  the  control  numerical  simulation.  The  analy¬ 
ses  showed  that  early  cell  development  occurred  directly  over  or  slightly  east  of  the 
ridge  lines.  Convective  initiation  was  favored  in  these  locations  by  upslope  flow 
through  convergent  valleys  located  on  the  eastern  and  southeastern  slopes  of  the 
mountain.  As  the  storm  developed,  circulations  from  existing  cells  affected  the  initia¬ 
tion  locations  of  new  cells.  For  example,  several  short-lived  cells  developed  along 
the  outflow  boundaries  from  older  cells,  and  in  some  cases  the  low-level  flow  induced 
by  the  existing  cells  enhanced  or  created  upslope  and/or  convergent  flows  which  in 
turn  helped  focus  new  cell  development. 

The  observed  storms  generally  had  lifetimes  of  1.0  to  1.5  hr  with  individual  cell 
lifetimes  ranging  from  20  to  60  min.  Many  of  the  cells  exhibited  periodic  intensifica¬ 
tion  and  weakening  with  a  period  of  6  to  12  min.  Radar  returns  from  the  strongest 
cells  exceeded  50  dBZ  several  times,  but  maximum  returns  from  most  cells  ranged 
from  30  to  40  dBZ.  Most  cells  never  extended  higher  than  10.0  km  MSL,  however, 
the  strongest  cells  frequently  extended  to  over  12.0  km  MSL.  Most  cells  contained 
single  updrafts  approximately  2.0  km  wide  with  maximum  vertical  velocities  of  12.0 
m  s  '  or  less,  however,  the  strongest  updrafts  reached  18.0  m  s'.  At  least  one  cell 
contained  multiple  low-level  updrafts  which  merged  into  a  single  updraft  at  mid¬ 
levels.  Mature  cells  contained  downdrafts  with  large  areal  coverages  and  maximum 
vertical  velocities  frequently  exceeding  6.0  m  s  '  with  an  absolute  maximum  of  12.0 
m  s  '.  During  one  cell’s  mature  stage,  the  interaction  between  the  updraft  and  the 
downdraft  appeared  to  be  symbiotic  in  a  manner  similar  to  that  often  observed  in 
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long-lived  severe  storms.  In  this  case,  the  precipitation  and  its  associated  downdraft 
were  located  behind  and  below  the  updraft  (relative  to  the  low-level  inflow)  and  thus 
did  not  interfere  with  the  inflow  or  the  updraft. 

A  strong  interrelationship  existed  between  all  the  cells  in  the  analysis  domain.  In 
the  mature  stage,  the  outflow  from  downdrafts  appeared  in  some  cases  to  enhance 
low-level  convergence  that  either  focused  the  development  of  new  cells  or  helped  sus¬ 
tain  existing  cells.  In  several  instances  the  updrafts  of  two  or  more  cells  merged  into 
one  updraft  at  upper-levels.  Often  the  updraft  of  a  newer  cell  developed  in  the  low- 
to  mid-levels  and  fed  into  a  higher-level  updraft  associated  with  an  older  cell.  Tire 
periodic  intensification  noted  above  appeared  to  alternate  among  the  cells.  In  other 
words,  when  one  cell  intensified,  the  others  weakened  slightly.  Finally,  the  observa¬ 
tions  showed  several  cases  where  older  or  weaker  cells  were  overtaken  and  absorbed 
by  newer  and/or  stronger  cells. 

Comparison  of  the  evolution  of  the  control’s  simulated  storm  with  observations 
and  with  Lang's  (1991)  microphysical  retrieval  show  that  the  simulation  accurately 
captured  the  major  features  observed  by  the  radars  and  accurately  depicted  the  micro¬ 
physical  evolution  of  the  storm  cells.  In  fact,  given  the  limitations  imposed  by  the 
use  of  a  single  initialization  sounding  and  the  two-dimensional  nature  of  the  simu¬ 
lation,  the  simulation  performed  above  expectations.  The  simulation  accurately  de¬ 
picted  the  development  of  a  multi-cellular  storm  with  horizontal  and  vertical 
dimensions  very  similar  to  those  observed.  The  simulation  captured  low-Ievel  inflow 
and  convergence  regions  as  well  as  upper-level  outflow  regions.  In  general,  the 
strength  of  the  simulated  cells  was  close  to  that  observed  by  the  radars.  The  imposed 
limitations  did,  however,  introduce  several  discrepancies  between  the  simulation  and 
the  observations.  These  included  an  initial  delay  in  convective  development  of  about 
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20-30  min  from  that  observed;  more  rapid  development  than  observed;  slightly  stron¬ 
ger  updrafts  and  weaker,  less-organized  downdrafts  than  observed;  a  tendency  for  de¬ 
velopment  farther  west  than  observed;  and  the  development  of  an  upper-level  total 
condensate  mixing  ratio  maximum  with  no  corresponding  maximum  in  the  radar  re¬ 
flectivity  fields. 

The  simulated  microphysical  evolution  compared  favorably  with  that  described  by 
Lang  (1991).  In  particular,  the  simulation  accurately  depicted  the  importance  of 
warm  rain  processes  early  in  the  cell  development  and  the  subsequent  rapid  increase 
in  the  importance  of  the  ice  phase  processes.  The  simulation  also  differed  from  the 
retrieval  results  in  several  areas.  The  simulated  mixing  ratios  of  ice  crystals  and  ag¬ 
gregated  snow  flakes  were  much  higher  than  those  retrieved,  while  the  simulated  rain 
mixing  ratios  were  much  lower  than  those  retrieved.  Additionally,  the  maximum 
simulated  rain  mixing  ratio  occurred  after  the  graupel  maximum,  while  in  the  retriev¬ 
al,  the  rain  maximum  occurred  first. 

The  first  two  experiments  tested  the  sensitivity  of  the  simulated  storm  to  the  initial 
wind  profile.  The  first  experiment,  EXPl,  used  an  initial  wind  profile  consisting  of 
calm  winds  throughout  the  depth  of  the  atmosphere.  Comparison  of  the  results  of  this 
simulation  with  those  of  the  control  shows  that  the  initial  wind  profile  plays  an  impor¬ 
tant  role  in  the  initiation  and  evolution  of  the  simulated  storms.  The  simulated  cells 
in  EXPl  developed  directly  over  the  ridge  lines  instead  of  to  the  west  of  the  ridges  as 
in  the  control.  The  location  of  cell  development  in  EXPl  is  more  consistent  with  the 
observed  cell  development  than  is  the  control.  EXPl's  storm  cells  developed  earlier, 
but  grew  at  a  slower,  steadier  pace  than  the  control’s.  EXPl's  initially  calm  winds 
limited  the  turbulent  transfer  of  heat  away  from  the  surface,  allowing  maximum  solar 
heating  and  rapid  convective  development.  The  control’s  more  rapid  growth  may  be 
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attributed  to  a  relatively  strong  flow  at  mid-  to  upper-levels  which  served  as  an  "ex¬ 
haust  mechanism"  for  the  developing  cells.  EXPl's  storm  contained  two  cells  of 
nearly  equal  strength  and  height  through  most  of  the  development  period.  The  reason 
for  this  difference  from  the  control  was  riot  clear.  Finally,  EXPl's  microphysical  and 
vertical  velocity  evolution  showed  a  much  greater  variability  than  that  seen  in  the  con¬ 
trol. 

The  results  of  EXPl  showed  that  the  initial  wind  profile  strongly  influences  the 
rate  of  development  of  all  cells  and  may  even  inhibit  the  development  of  some  cells. 
The  wind  profile  also  exerts  influence  over  the  microphysical  evolution  and  vertical 
velocity  evolution  of  the  storm.  Finally,  the  initial  wind  profile  plays  a  significant 
role  in  the  convective  initiation  location. 

The  second  experiment,  EXP2,  was  initialized  with  a  wind  profile  containing 
wind  speeds  50%  of  those  in  the  control.  This  simulation  developed  serious  boundary 
condition  instabilities  which  rendered  the  simulation  useless  for  the  purposes  of  this 
study. 

The  third  experiment,  EXP3,  was  conducted  with  the  model's  short-wave  radi¬ 
ation  package  deactivated  to  test  the  effect  on  the  simulated  storm  of  a  "no  sun”  sce¬ 
nario.  In  this  simulation,  very  small,  non-convective  clouds  form  at  nearly  the  same 
locations  as  the  initial  clouds  in  the  control—slightly  west  of  the  ridge.  A  weak,  but 
steady,  updraft  supported  the  clouds.  The  updraft  developed  near  a  nearly  stationary 
area  of  low-level  convergence  which  resulted  from  interaction  of  the  predominantly 
easterly  flow  with  a  weak  westerly  flow  near  the  surface.  The  westerly  flow  devel¬ 
oped  as  a  result  of  a  dynamically  induced  mountain  wave  circulation  on  the  lee  side  of 
the  mountain.  The  results  of  this  simulation  show  that  solar  heating  plays  the  most 
important  role  in  the  development  of  these  mountain  thunderstorms.  The  results  of 
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this  experiment  combined  with  those  of  EXPl  also  show  that  the  environmental  wind 
field  plays  a  crucial  role  in  determining  the  convective  initiation  locations. 

The  fourth  experiment,  EXP4,  tested  the  sensitivity  of  the  simulated  storm  to  a 
restriction  to  warm  rain  microphysical  processes.  In  this  experiment  all  the  model's 
ice  phase  microphysical  parameterizations  were  deactivated.  The  results  from  early  in 
this  simulation  agreed  very  well  with  those  of  the  control.  Later  in  the  simulation 
EXP4's  results  showed  the  simulated  storm  contained  weaker  updrafts  and  lower  total 
condensate  mixing  ratios.  Additionally,  the  cells  did  not  extend  as  high  as  those  in 
the  control.  EXP4's  rain  shafts  were  much  weaker  than  in  the  control,  and  the  total 
condensate  mixing  ratio  maxima  noted  near  the  cell  tops  in  the  control  were  much 
weaker  in  EXP4.  The  good  early  agreement  between  the  results  of  this  experiment 
and  those  of  the  control  confirm  the  importance  of  warm  rain  processes  in  the  storm's 
early  development.  The  growing  differences  between  the  two  simulations  show  that 
as  the  storm  matured,  the  ice  phases  became  increasingly  important  to  the  evolution  of 
the  storm's  development.  Indeed,  the  absence  of  the  ice  phase  processes  and  their 
associated  latent  heat  release  processes  in  EXP4  prevented  the  storm  from  reaching 
the  strengths  achieved  in  the  control.  The  lack  of  the  efficient  precipitation-producing 
mechanisms  associated  with  the  ice  phases  prevented  the  rain  shafts  in  EXP4  from 
reaching  the  strengths  achieved  in  the  control.  This  experiment  demonstrated  the  im¬ 
portance  of  the  ice  phases  in  the  development  of  these  mountain  storms— particularly 
in  their  middle  to  late  stages  of  development. 

The  final  experiment,  EXP5,  tested  the  sensitivity  of  the  simulated  storms  to  the 
initial  base-state  moisture  profile.  The  initial  moisture  profile  in  this  case  was  moist¬ 
ened  by  up  to  0.5  g  kg  '  in  the  lowest  1.0  km.  This  simulated  storm  generally  fol¬ 
lowed  that  of  the  control,  however,  one  of  the  very  minor  cells  in  the  control 
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dominated  EXP5's  simulation  for  a  time.  In  general  EXPS's  simulated  storm  devel¬ 
oped  at  a  steadier  pace  than  the  control  and  eventually  contained  stronger  updrafts  and 
higher  graupel  mixing  ratios  than  the  control.  EXPS's  cells  did  not  reach  the  same 
heights  as  the  control's,  and  as  in  EXPl,  two  cells  developed  to  nearly  the  same 
strength.  This  experiment  clearly  shows  the  importance  of  low-level  moisture  to  the 
evolution  of  these  mountain  storms.  In  particular,  the  moisture  profile  seemed  to  in¬ 
fluence  the  ability  of  some  cells  to  develop  beyond  cumulus  stages.  Additional  low- 
level  moisture  contributes  to  stronger  updrafts  and  slightly  higher  condensate  mixing 
ratios,  but  interestingly  does  not  increase  the  height  of  the  cells.  The  lower  cloud 
heights  may  be  attributed  to  the  increase  in  the  number  of  larger  cells.  The  increased 
competition  among  the  larger  cells  may  prevent  any  of  them  from  reaching  the  heights 
seen  in  the  control. 

In  summary,  the  results  of  the  numerical  simulations  showed  the  following  condi¬ 
tions  and  processes  to  be  most  important  (in  descending  order)  to  the  development  and 
maintenance  of  the  mountain  storms  under  study  here:  solar  heating,  ice  phase  micro¬ 
physical  processes,  environmental  wind  profile,  and  low-level  moisture  profile.  The 
environmental  wind  profile  was  the  single  most  important  factor  in  the  convective  ini¬ 
tiation  location. 

The  most  obvious  limitation  of  the  simulations  conducted  for  this  study  was  their 
two-dimensional  nature.  This  limitation  prevented  a  detailed  quantitative  analysis  of 
the  validity  of  the  control  simulation  compared  with  the  Doppler  analyses.  It  also 
prevented  the  simulation  from  capturing  the  channelling  of  the  wind  flow  up  the  vari¬ 
ous  valleys  in  the  Magdalena  Mountains.  The  observations  demonstrate  the  impor¬ 
tance  of  this  phenomenon  in  focusing  the  convective  initiation  and  in  contributing  to 
the  maintenance  of  the  storm  development. 
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Future  studies  of  this  nature  should  attempt  to  overcome  these  limitations  by  em¬ 
ploying  fully  three-dimensional  simulations.  The  three-dimensional  simulations 
should  also  be  conducted  with  larger  domains  so  as  to  further  minimize  the  boundary 
condition  effects  on  the  areas  of  interest  such  as  those  noted  in  EXP2.  Future  studies 
should  also  attempt  to  employ  higher  resolution  grid  systems  in  both  the  horizontal 
and  vertical  directions.  A  higher  resolution  grid  system  will  naturally  require  a  high¬ 
er  resolution  input  terrain  field.  Such  high  resolution  grids  will  minimize  the  need  for 
extensive  parameterizations  of  processes  so  important  to  the  development  of  these 
small  storms. 

Comparison  of  numerical  simulation  results  with  Doppler  analyses  would  be 
markedly  simplified  with  the  employment  of  a  technique  to  synthesize  radar  reflectivi¬ 
ties  based  on  the  simulated  microphysical  storm  structure.  Finally,  the  accuracy  of 
three-dimensional  simulations  could  be  improved  with  improved  initial  conditions. 
This  study  used  a  single  profile  to  initialize  the  simulation.  A  complex  four¬ 
dimensional  data  assimilation  system  using  all  available  data  including  Doppler  radar 
observations  and  perhaps  retrieved  microphysical  parameters  could  provide  the  model 
with  a  significantly  more  realistic  representation  of  the  base-state  atmosphere.  With 
further  model  developments,  such  high-quality  analyses  could  provide  information 
required  to  "nudge”  the  simulation  as  it  is  progressing  toward  a  more  realistic  solu¬ 
tion. 
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